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Barium (Ba!*134) obtained by dn and d2n reactions of 
spectroscopically pure cesium chloride was separated in 
quantities of <110-* gram without addition of common 
barium. The procedure studied involves adsorption of 
barium on hydrous ferric oxide under conditions at which 
cesium ions remain in solution. The effect of (a) pH of 
solution and (b) varying quantities of ferric hydrous oxide 
on the fraction of Ba adsorbed from solution containing 
<10-* g Ba, keeping other factors constant, showed that 
the fraction of Ba adsorbed, y=k[Z(X —Xo)}, Z being the 
quantity of adsorbent (from 110-5 to 3X10 g atom 
Fe); X, the OH~ concentration in solution (from 3X 10-9 
to3X10-*); Xo, thecritical concentration of OH~ necessary 


to maintain the adsorbent in the solid state. When pH was 
kept constant the value of \ for y=k’Z* was found to be 
0.35, and when the quantity of adsorbent was constant 
the A for y=k’""(X —Xo)* was 0.33..The fact that the two _ 
experimental slopes were so nearly the same in spite of 
certain variations in experimental conditions is taken as 
support for two assumptions made in the derivation of the 
expression viz., (1) the fraction of barium adsorbed from a 
very dilute solution depends on the quantity of adsorptive 
centers in the solid phase and (2) the quantity-of adsorptive 
centers in the solid phase available for barium adsorption 
is proportional to the OH~ concentration in solution, other 
factors being constant. Adsorption procedure is described. 





INTRODUCTION 


N previous work it has been shown that ab- 
sorbents such as hydrous ferric oxide will re- 
move part of such divalent ions as strontium, 
barium, and radium from their solutions at any 
feasible low concentration if experimental condi- 
tions are appropriate.)? — 

Consequently it was suggested that chemisorp- 
tion might be applied first, to isolation, without 
carriers, of minute quantities of radioactive and 
stable species produced by nuclear transmuta- 
tions, from elements of different atomic number 
and, second, to determination of quantities as low 
as 1X10" g atom, that is, quantities less than 
those measurable by the most sensitive previ- 
ously known schemes.® 

1J. D. Kurbatov, J. L. Kulp, and Edward Mack, Jr., J. 
Am. Chem. Soc. 67, 1923 (1945). 

2M. H. Kurbatov and J. D. Kurbatov, J. Am. Chem. 
Soc. 69, 438 (1947). 


5M. H. Kurbatov and J. D. Kurbatov, Phys. Rev. 71, 
466 (1947), 
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The adsorption of divalent ions on hydrous 
ferric oxide is influenced by such variables as pH, 
salt concentration, quantities of adsorbent and 
adsorbate, volume of liquid phase, temperature 
and duration of adsorption. However, the depend- 
ence of adsorption of divalent ions in extreme 
dilution on the above variables is insufficiently 
known especially for purposes of isolation, with- 
out carriers, of species produced by nuclear 
transmutations. 

In the present work a scheme for isolation of 
Ba".134 (produced by activation of cesium with 
10.5-Mev deuterons) in quantities of less than 
10-° g has been studied with reference to its 
dependence on hydrogen ion concentration and 
quantity of adsorbent. 


PROCEDURE 


Ba".134 as the chlorides were separated (as 
described later in this paper) in quantities of less 
than 10-" g atom from the target material; the 




















TABLE I. 
0.1N iN 0.1N 
Fe HCl Water NH.OH NH.OH 
g atom./exp. ml ml ml ml 
6X 10-* 0.0 23.0 1,50 S.o0 
1x10- 16.7 6.3 1.60 5.90 
1X10-5 19.7 3.3 1.65 5.40 








target being spectroscopically pure cesium chlo- 
ride and the holder platinum. The purity of tracer 
(Ba!* as chloride*) was confirmed by (1) evapo- 
ration of the whole tracer in a micro beaker and 
weighing, with the result that no weight could be 
detected to 10-® g (10-* g atom) and (2) by 
adsorption procedure indicating less than 10~° g 
atom Ba. That is, 40 percent of the tracer solu- 
tion was added to 1X 10-5 g atom Fe, the solution 
titrated to pH 8, and the fraction adsorbed found 
_to be 25 percent. Then ten percent of the residual 
solution after adsorption (7.5 percent of the total) 
was added to a fresh portion of iron and volume 
and pH adjusted as before. The second portion 
adsorbed was the same fraction, i.e., 25 percent, 
indicating the original concentration was in the 
region in which the fraction adsorbed is inde- 
pendent of concentration.’ 

The quantity of tracer used in individual ad- 
sorption experiments did not exceed 110+ of 
the total tracer obtained after bombardment, 
that is not more than 10-* g atom Ba. Also the 
bombarded target was spectroscopically pure 
cesium chloride so it contained less than 10-6 g 
barium per 0.1 g of activated material and this in 
turn was diluted to less than 10-° g (or 10-” g 
atom). Considering a maximum of 2000 counts 
used per adsorption experiment with a counter 
efficiency of 9 percent the amount of radioactive 
Ba'® was of the order of 10-4 g atom. 

The purity of the tracer with reference to 
radioactivity was verified by following its dis- 
integration om an electrometer. The curve ob- 
tained corresponded to Ba‘ (TJ = 38.8 hours) and 
was the same as that obtained from a sample 
purified with stable barium carrier added. 

For the adsorption studies on the effect of pH 
and quantity of adsorbent, in addition to the 
purified tracer solution just described, solutions 
of stable barium chloride, ferric chloride, hydro- 


* Ral is stable. 
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chloric acid and ammonia, were required. Spectro- 
scopically pure barium chloride was diluted to 
1<10-" g atom/ml in 0.01N hydrochloric acid. 
Since the spectroscopic analysis indicated no 
impurities present to the extent of 10° g/g 
barium, the diluted barium could not contain 
more than 10~'* g atom/ml of impurities from the 
barium chloride. The water used in preparing all 
solutions was triple distilled and checked by pH 
measurement. The ferric chloride used as source 
of adsorbent was purified by repeated precipita- 
tion of the hydrous oxide at low #H, under which 
conditions divalent impurities are removed in 
solution. 

The ammonia solutions used in titrating the 
samples were prepared by ammonia gas transfer 
from concentrated solutions into triple-distilled 
water in a closed vessel. 

For the individual adsorption determinations 
under a given set of conditions a series of samples 
was made up of stable barium, tracer, and ferric 
chloride. The required volumes of these solutions 
having been transferred into outside ground 
weighting bottles, water or hydrochloric acid or 
both were added to adjust salt concentration and 
volume. Then the samples were rapidly titrated 
with ammonia to the predetermined pH and 
samples closed and set aside to allow adsorbent to 
settle. After a given time interval, portions of 
solution were removed by pipette, these portions 
evaporated in 4-cm Petri dishes, and the activity 
measured on a G-M counter. 

More specifically, in the case of the determina- 
tion of pH effect the factors which were held 


TABLE II. Effect of amount of hydrous ferric oxide on 
adsorption of barium. 


Constant factors: Temp. = 25°+1°; total Ba=1X10-' g 


atom, ~H=8.01+0.01; NH,Cl conc. =5.4XK10- N; Vol. 
=32.3—34.4 ml. 








Total Ba left 





in sample Ba adsorbed 
Ba adsorbed after 10 ml after total 
Total Fe in in 4 hours sol'’n. removed of 23 hours 
g atom g atom X10" gatomX10-2 g atom X10-" 

1x 10-5 2.41 7.71 4.42 
2x10-5 2.88 7.87 5.51 
3x 10-5 3.40 7.96 5.67 
6X 10-5 4.37 8.29 6.84 
1X10 5.92 8.82 7.61 
2x10-* 6.85 9.11 8.40 
3X10 7.95 9.37 8.98 
6X 10-4 9.20 9.76 9.61 
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made up with 1-ml barium chloride solution con- 
taining 1<10-" g atom Ba in 0.01 N HCl, 1 ml 
FeCl; solution containing 1 X 10-* g atom Fe, 1-ml 
tracer in 0.01 N HCl and 5 ml of water. Then each 
was titrated to a different H using both 0.1 N 
and 0.01 N NH,OH, so that the final volumes 
were all within the range of 12 to 13 ml. 

In the case of the experiments using different 
quantities of adsorbent, the samples were made 
up with each containing 1 ml of barium chloride 
solution having 1 X 10-" g atom Ba in 0.01 N HCl, 
1 ml of tracer in 0.01 N HCl, and 1 ml ferric 
chloride solution containing the various quanti- 
ties of iron shown in Table I. Since the varying 
quantities of ferric chloride meant varying 
amounts of acid from hydrolysis, the more dilute 
iron solution had acid added to compensate for 
the greater acidity in the more concentrated iron 
solutions and varying amounts of water were 
added to maintain constant final volume. All 
then were titrated to the same pH. 

Typical samples were as shown in Table I. 


DEPENDENCE OF ADSORPTION OF BARIUM IN 
EXTREME DILUTION ON HYDROGEN ION 
CONCENTRATION AND QUANTITY 
OF ADSORBENT 


The adsorption of divalent ions on a given 
quantity of hydrous ferric oxide may be looked 
upon as competitive with all positive ions present 
in solution. With very dilute solutions of divalent 
ion the change in hydrogen ion concentration in 
solution—resulting from an adsorption—would 
not be measurable. However, as the #H of the 
solution is changed by the addition of an appro- 
priate reagent it would be expected that the 
amount of adsorption would be inversely related 
to the hydrogen ion or directly related to the 
hydroxyl ion concentration in solution. 

If we assume that (1) the fraction of barium 
adsorbed (y) depends only upon the quantity of 
adsorptive centers in the solid phase when other 
factors such as salt concentration, temperature, 
volume, and duration of adsorption are held con- 
stant (in dilution below 10—'° g atom/ml the frac- 
tion adsorbed is known to be experimentally 
independent of total barium) and (2) the avail- 
able quantity of adsorptive centers (S) in the 
solid phase for a given quantity of Fe (Z) is 
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constant are shown in Table II. Each sample was 
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Fic. 1. Effect of amount of hydrous ferric oxide on 
adsorption of barium. 


inversely proportional to hydrogen ion concen- 
tration, or directly proportional to the hydroxy] 
ion concentration in solution, that is, 


S/Z=CTOH-], 


where C is a proportionality constant; then since 
both S and Z refer to the solid phase the effective 
[OH~] is that above a critical value at which the 
solid phase can be maintained. This critical 
[OH~] is designated Xo and the higher [OH~] as 
X so the above expression becomes 


S=KZ(X —X). 


If a fraction y of barium is adsorbed in the 
presence of a quantity S of adsorptive centers in 
the solid phase, then an increase of S by dS 
would give rise, according to (1), to a propor- 
tionate increase in y by dy. Therefore, 

d 
y S 


Integration gives 


Iny =A In[Z(X — Xo) ]+const., 


TABLE III. Effect of | a on adsorption of barium ion by 
hydrous ferric oxide. 


Constant factors: Temp., 25°+1°; total Ba, 1X10-" g 
atom; Fe, 1X10-* g atom; NH,CI conc,, 2.5X107? N; 
Vol., 13.340.3 ml. Duration of adsorption ~3} hours. 








Activity in 5 ml 





pH : of solution Tracer in exp. Ba adsorbed 
Counts/min. Counts/min. g atom X 10" 
5.48 303 878 1.02 
6.00 290 882 1.45 
6.53 259 884 2.32 
6.89 239 858 2.79 
7.48 166.5 853 4.82. 
7.98 114.5 880 6.61 
8.45 26.8 856 9.15 
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and 
y=k[Z(X ~Xo) }. 


From this expression it follows that a log-log 
plot of adsorbed barium against total iron present 
when #H and other factors are constant should 
give a straight line. 


EXPERIMENTAL RESULTS 


The effect of varying the quantity of adsorbent 
while keeping the H, total barium, salt, volume, 
and duration of adsorption constant is shown in 
Table II and Fig. 1. 

The range of iron quantities which could be 
used readily was rather limited. The lower limit 
was an amount which would settle well under the 
conditions used and the upper limit an amount 
which would introduce acid in quantity not too 
much in excess of that used in the other series of 
experiments (to keep salt concentration com- 
parable). As a check on the results obtained for a 
four-hour interval a second portion was removed 
from each sample twenty-three hours after the 
titration. 

Within the range of iron quantities from 
1X10-5 to 3X10~ g atom Fe the log-log plot of 
adsorbed barium vs. iron is an experimentally 
straight line, in agreement with the equation, but 
above this range the curve breaks. Extrapolation 
of the straight line would indicate 100 percent 
adsorption by 5.5X10~ g atom of Fe but 
6X 10-* g atom Fe does not adsorb all the barium 
even at the longer time interval. 

The curve at 23 hours verifies the fact that the 
curve at 4 hours is not a straight line throughout 
and that the range of iron quantities which gives 
a straight line is dependent on duration of adsorp- 
tion in addition to other factors mentioned.** In 
addition, from comparison of the two curves on 
Fig. 1, it can be seen that (with extreme dilutions 
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of barium and limited quantities of hydrous 
ferric oxide) the adsorption occurs sufficiently 
slowly to show its dependence on time. Equi- 
librium is not established in a period of a few 
hours. 

In Table III and Fig. 2 may be seen the result 
of changing the fH of the solution on barium 
adsorption when the total barium, iron, salt, 
volume, and duration of adsorption were kept 
constant. The log-log plot is an experimentally 
straight line over the range of H values used 
(pH 5.5 to 8.5). 

Extrapolation of this curve indicates that at 
pH 3.3 (2 10-"[OH— }), the region in which iron 
precipitates, the amount of adsorption would be 
two percent, however, since the iron does not 
settle well enough at this pH for accurate meas- 
urements, the deviation from a straight line was 
not established. On the other end of the curve 
extrapolation indicates 100 percent adsorption at 
3.510-§ [OH-] or pH 8.56 and this is experi- 
mentally not observed. Anexplanation of this isthe 
time factor. So the application of the equation for 
PH effect must be considered as limited to about 
the range of [OH~] studied. It is to be noted that 
the value of Xo (2X10-") becomes insignificant 
in this range. 

The slope of the straight line part of the curve 
of the “‘iron effect’’ is 0.35 while that of the “pH 
effect’”’ curve is 0.33. Considering the differences 
in volume and some experimental variations in 
salt concentrations and time intervals, the simi- 
larity in slopes is taken as supporting the assump- 
tions made-in the derivation of the equation. 


TARGET CHEMISTRY OF CESIUM ACTIVATED 
WITH DEUTERONS 
The bombarded cesium chloride was allowed to 
stand for several hours for disintegration of the 
radioactive chlorine. 
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** The short half-life of the tracer required study of adsorption over short time intervals. 
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The target material was then washed from the 
platinum holder with triple-distilled water, dis- 
solved in hydrochloric acid (bH~2) and filtered 
through analytical paper. The impurities insolu- 
ble in hydrochloric acid were rejected (see Fig. 3). 
To the acid filtrate was added 1X10-* g atom 
iron as ferric chloride and the solution titrated 
with ammonia to pH 10. The filtrate from this 
was rejected. 

The precipitate containing 91.5 percent of the 
active barium was dissolved in acid and the 
procedure repeated twice. 

The additional loss of 1.8 percent of active 
barium was less than previously lost because of 
lower salt concentration and volume. 

Potassium, activated by deuterons, was purified 
by similar adsorption procedure, and ~300 
counts added to the barium "4 as a tracer for 
monovalent ions (Cs) that might be present in 
the barium solution. 

The adsorption procedure with hydrous ferric 
oxide was carried out again. 

Finally the iron was precipitated at pH 5.4, 
and the filtrate containing 83.1 percent of the 


initial barium activity was freed from salt by 
sublimation, dissolved in 0.01 N HCl, filtered, the 
operation repeated, the sample weighed and 
found to contain no weighable material. 

The active material was dissolved in 0.01 N HCl 
and diluted 10‘ times for use as tracer. It is 
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assumed that the tracer contained minute quanti- 
ties of potassium and cesium. 
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Dimethyltrisulfide (CH3)2S3 has been investigated by the electron diffraction method. The 
molecule was found to have the chain configuration, with the following interatomic distances 


and angles: C~H=1.09A (assumed), C—S=1.78+0.04A, S—S=2.04A, +0.02A, ZH—C-—S 
= 112° (assumed), 2S—S—S=104°+5°, and ZC—S—S=104°+5°. The results also suggest 
that the C—S bonds oscillate with rather large amplitude about a mean position approximately 


106° out of the plane of the three sulfur atoms. 











HE present electron diffraction investiga- 
tion was undertaken mainly to determine 
whether dimethyltrisulfide is a simple trisulfide 
(I), a thiodisulfide (II), or a dithiosulfone (III). 
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Knowledge of this structure is of interest in 
connection with the problems of the constitu- 
tions of the inorganic persulfides and polythio- 
nates,! as well as of the alkyl polysulfides,? and 
the Thiokols.* With the exception of the tri- 
thionate ion, O;S—S—SO3;-, there is no direct 
structural evidence which settles the question 
whether the molecules of these substances have 
the chain configuration or are related to II or 
III. In the case of the disulfides H2S2, SeCl2, and 
(CH3)2Se, the existing evidence from electron 


* The work herein was done under Contract OEMsr-753 
between the National Defense Research Committee and 
the California Institute of Technology. 

** Contribution No. 1108 from the Gates and Crellin 
Laboratories of Chemistry, California Institute of Tech- 


nology. 
iD. M. Yost and H. Russell, Systematic Inorganic 
Chemistry (Prentice Hall, Inc., New York, 1944), p. 361ff. 
2R. Connor, Organic Chemistry (John Wiley and Sons, 
Inc., New York, 1943), edited by H. Gilman, p. 864. 
3S. M. Martin and J. C. Patrick, Ind. Eng. Chem. 28, 
1144 (1936). 
4W. H..Zachariasen, Phys. Rev. 37, 1295 (1931). 









diffraction favors chain structures for all three 
molecules.5~7 


EXPERIMENTAL 


The sample used in this investigation was pro- 
vided by Professor R. C. Fuson. It was pale 
yellow in color and was used without further 
purification. Excellent electron diffraction photo- 
graphs were made in the apparatus described 
by Brockway® with the use of the high tempera- 
ture nozzle.® At the end of each run the boiler 
and remaining liquid were examined; no corro- 
sion of the boiler and only slight darkening of 
the liquid were observed, indicating that very 
little decomposition occurred. 


INTERPRETATION OF THE PHOTOGRAPHS 


The characteristic features of the photographs 
are represented by curve V of Fig. 1, which was 
drawn in accordance with the measured ring 
diameters (except that the unobservable first 
maximum, which is not structure sensitive, was 
added arbitrarily), so as to have the general 
properties of theoretical intensity curves calcu- 
lated with the simplified expression 


ZLi, _{F 
exp(—4;;q?) sin( “ra. 





Ig=X 


“i ij 


5K. J. Palmer, J. Am. Chem. Soc. 60, 2360 (1938). 

6D. P. Stevenson and J. Y. Beach, J. Am. Chem. Soc. 
60, 2872 (1938). 

7P. G. Ackerman and J. E. Mayer, J. Chem. Phys. 4, 
377 (1936). 

8L. O. Brockway, Rev. Mod. Phys. 6, 234 (1936). 
Camera distance, 10.91 cm. Electron wave-length, 0.0610A 
calibrated with respect to zinc oxide smoke (C. S. Lu and 
E. W. Malmberg, Rev. Sci. Inst. 14, 271 (1943)) and 
corrected for film expansion. 

®L. O. Brockway and K. J. Palmer, J. Am. Chem. Soc. 
59, 2181 (1937). 



















STRUCTURE OF DIMETHYLTRISULFIDE 


The radial distribution integral,!° curve RD of 
Fig. 1, was calculated by use of the formula 


100 


X= 1(g) exp(—agq’) sin[ (x/10)qr ], 


q=1, 2, eee 


rD(r) = 


with exp(—a(100)?) =0.1. The integral has max- 
ima at 1.78A, 2.04A, 2.45A, 3.02A (very broad, 
asymmetric peak), and 4.0A (very broad, small 
peak). Of these peaks the first two may be 
interpreted as representing, respectively, the 
bonded C—S and S—S distances, the third the 
non-bonded S...H distance, and the asym- 
metric fourth the non-bonded C...S and S...S 
interrelations. The interpretation of the weak 
maximum at 4.0A is discussed below. 

An extensive investigation of the straight 
chain structure I by the correlation procedure 
was made. For structure I four parameters are 
required to define the structure if the positions 
of the hydrogen atoms are assumed and if the 
long, rotation-dependent, C...S distance is neg- 
lected. These may be chosen as: (1) bonded 
>—S; (2) bonded S—S; (3) non-bonded S...S; 
and (4) non-bonded C...S. 

The first two of these parameters may be 
determined directly from the radial distribution 
integral. As a check on the accuracy of this 
determination, twenty-one intensity curves, some 
of which are shown in Fig. 1, were calculated for 
models with varying ratios between these two 
distances. In these models the S—S distance was 
held at 2.04A and the non-bonded C...S and 
S...S distances were varied in a systematic 
fashion so that their average was in general 
agreement with the broad maximum of the in- 
tegral at 3.0-3.2A. Best agreement between the 
calculated curves and the appearance of the 
photographs was obtained for models with C—S 
=1.78A as was expected, changes of more than 
0.03A giving marked discrepancies regardless of 
the particular values used for the non-bonded 

..S and S...S distances. 

The determination of the two remaining 
parameters is a more difficult problem since the 
terms which are involved in their variation con- 
tribute only about 20 percent of the total 
scattering. Since the term representing the two 

104. R. Spurr and V. Schomaker, J. Am, Chem. Soc. 64, 


2693 (1942); b. P. A. Shaffer, Jr., V. Schomaker, and L: 
Pauling, J. Chem. Phys 14, 659 (1946). 
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Fic. 1. Electron diffraction o— for dimethyltrisulfide. 
Chain models with C—H=1.09A, S—S=2.04A, S...H 
=2.45A, C...S=3.02A,S...S= 3.22A( ZS— —-S- —S=104°), 
and bdc...7=0.00016; the other molecular parameters are 
as follows. 


Model bs...H (zC-—S-—S) 


d 0 102° 45’ 
B d 0.0003 103° 30’ 
Cc . 0.0003 104° 
D ° 0.0003 105° 
E ° 0 106° 30’ 


All other 6;;=0. The long C....S term at 4.00A, was 
included with b=0.0009, for models B, C, and D, and was 
omitted for the other models. 


non-bonded C...S distances and the term for 
the non-bonded S...S distance are of nearly 
equal weight, it is impossible to distinguish 
between them. For example, a curve calculated 
for a model with z2S—S—S= zS—S—C=104° 
is indistinguishable from one calculated for a 
model with zS—S—S=953° and zS—S—C 
=115°. However, it is appropriate to choose 
S...S<C...S so that the sulfur bond angles 
ZS-S-—S and ZS—S-—C have similar values, 
and this is done in the following analysis. 
Twenty-two scattering curves, some of which 
are shown in Fig. 2, were calculated for models 
with varying values of these distances. It was 
found that it is not possible to obtain theoretical 
scattering curves in agreement with the appear- 
ance of the photographs unless (1) a distance at 
4.0A as suggested by the radial distribution in- 
tegral is included in the intensity formula and 
(2) the S...S distance is given a rather severe 
temperature factor. The distance at 4.0A corre- 
sponds to the two longer non-bonded C...S 
distances in a non-planar model with the C—S 
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bonds rotated by 106° out of the plane containing 
the three sulfur atoms. This is just the angle of 
twist between the O—H bonds in hydrogen 
peroxide, as deduced from the crystal structure 
of H2O2-urea." Unless this distance is included in 
the calculation, the character of the first maxi- 
mum and adjacent minima is in definite disagree- 
ment with the appearance of the photographs 
and curve V, the theoretical curves having their 
sharp first minima just at the observed position 
of the weak first maximum, and a weak shoulder 
(which might be interpreted as corresponding to 
the observed first minimum) at the observed 
position of the second minimum. This discrep- 
ancy cannot be regarded as merely a difficulty 
in interpretation, in view of our knowledge of the 
appearance of electron diffraction photographs 
and their interpretation, but must be regarded 
as conclusive. On the curve for the accepted 
model, O, on the other hand, the positions of the 
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Fic. 2. Electron diffraction curves for dimethyltrisulfide. 
Chain models with C—H=1.09A, C—S=1.78A, S—S 
=2.04A, S...H=2.45A, and bc...4=0.00016; all other 
b;; =0 except for models J and M, for which bg,,.5 = 0.0003 ; 
the other parameters are as follows. 


Model C...8 S...S (zS—S-S) (zC-S-S) 
F 3.10 3.10 99° 108° 20’ 
G 2.96 3.10 99° 101° 30’ 
H 2.97 3.14 100° 30’ 101° 
I 3.02 3.17 102° 104° 
J 3.05 3.20 103° 105° 40’ 
K 3.02 3.22 , 102° 
M 3.05 3.25 105° 30’ 105° 40’ 
N 3.00 3.25 105° 30’ 103° 20’ 


The long C....S at 4.00A term was included with 
b=0.0009, for models J and M, and was omitted for the 
other models. 


uC, S. Lu, E. W. Hughes, and P. A. Giguére, J. Am. 
Chem. Soc. 63, 1507 (1941). 
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features corresponding to the measured first and 
second minima and the first maximum are satis- 
factory. The relative depths of the indicated 
first and second minima, also, are in better 
agreement with curve V than on the curves 
calculated without a 4.0A term. The agreement 
is not perfect either in this respect or for the 
magnitude of the first maximum, which on the 
photographs is clearly visible, although very 
weak, but on curve O is indicated only as a region 
of unusually small upward curvature; however, 
these discrepancies do seem to represent diff- 
culties of interpretation. 

Similar differences between the correct and the 
obvious interpretation of a weak inner ring—the 
ring appearing stronger on the photographs than 
would be suggested by the correct theoretical 
curve, or the apparent relative depths of the 
adjacent minima being different from that shown 
by the theoretical curve—have been previously 
observed in the case of numerous other com- 
pounds of which propane” is an example. Curve 
ZZ, for which the customary coefficients ZZ ;/r i; 
were replaced by Z;Z;, was calculated to illus- 
trate the notion that the apparent strength of 
the weak—and also relatively sharp—inner ring 
may be enhanced by a tendency of the eye to 
detect more readily the higher frequency com- 
ponents of the pattern. 

It is evident that, in general, the appearance 
of the innermost rings of electron diffraction 
photographs cannot be so simply interpreted in 
terms of simplified theoretical intensity curves as 
can the rest of the pattern, perhaps in part 
because of the very rapid change of total in- 
tensity with change in scattering angle at small 
scattering angles, and that it is never justifiable 
to make important use of an inner feature—such 
as weak inner ring—without reference to similar 
features on the photographs of other substances 
of known structure. In addition to propane, we 
have here made particular reference to the photo- 
graphs and curves of cyclopropane, ethylene 
imine, and ethylene oxide,!* which show some- 
what similar inner rings to that of dimethyl- 
trisulfide. 


12S, H. Bauer, J. Chem. Phys. 4, 406 (1936); L. Pauling 
and L. O. Brockway, J. Am. Chem. Soc. 59, 1223 (1937). 

18H. Pfeiffer, H. A. Levy, and V. Schomaker (to be 
published, J. Am. Chem. Soc.). 
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STRUCTURE OF DIMETHYLTRISULFIDE 


Curve O was calculated with the following 
values of the temperature factor 6: C—H: 
0.00016; S...H:0.0003; S...S:0.0003; long 
C....S:0.0009; all others zero, except for the 
remaining distances involving hydrogen, which 
were neglected. Curves C and P are identical 
with curve O except that the values zero and 
0.0009, respectively, were used for bg...s. Curve C 
is unsatisfactory in regard to a number of minor 
details in the region 35<qg<60 of which one, 
regarding the relative intensity of the fourth 
ring, can be seen most clearly in Fig. 3. Curve P 
is only slightly less satisfactory than curve O. 
This indication of a relatively large value re- 
quired for dg...s is rather surprising. We know of 
no other molecules in which the temperature 
effect is so large for a distance between heavy 
atoms bonded to the same atom. The large 
temperature factor required for the long C....S 
distance, on the other hand, is quite reasonable, 
and indicates that only a small potential restricts 
internal rotation at the temperature at which 
the photographs were made. 

Scattering curves for several of the models of 
type II and III which were calculated are shown 
in Fig. 3. The only important essential differ- 
ences between these models and chain structure 
models are differences of amplitude of the short 
C...S term and the long C....S term, of which 
the latter is most important for our argument, 
the essential long C....S term having only 
half-strength in type II models and being absent 
in type III models. It is seen that these curves 
are quite unsatisfactory with respect to the 
first ring. (At large scattering angles these curves 
and curve O are closely similar, none having an 
appreciable 4.0A contribution.) Since no varia- 
tion in molecular parameters whatever makes 
the initial portions of these curves correspond to 
the appearance of the photographs, we conclude 
that, to the accuracy of our observations of the 
inner rings, neither II nor III is the structure of 
dimethyltrisulfide. 

Our final values for the parameters of the 
satisfactory straight chain structure are: C—H 
=1.09A (assumed), C—S, 1.78+0.04A, S—S 
=2.04+0.02A, ZH—C—S=112° (assumed), 
ZS—S—S=104°+5°, and zC—S—S=104° 
+5°. (C...S=3.02+0.04A and S...S=3.22 
+0.04A.) The terminal CH;—S bonds are ro- 
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Fic. 3. Electron diffraction curves for dimethyltrisulfide. 
Molecular parameters for the various models are C—H 
=1.,09A, C—S=1.78A, S—S=2.04A, S...H=2.45A, long 
C....S=4.00A, and dc...—47=0.00016; the other param- 
eters are as follows. 

Chain models I 
Model bs...H 


K 0 

Cc 0.0003 

O 0.0003 

P 0.0003 
ZiZi 0.0003 
C...S =3.02A and S...S =3.22A for all chain models 


Dithiosulfone models III 
Model &...8 


Q 3.02 
R 3.12 
Ss 3.02 
T 3.22 


Thiodisulfide model II 
U 3.22 


b long C....S 
term omitted 


(zC-S~—C) 


tated about 106° from the plane of the sulfur 
chain. The quantitative comparison for model O 
is shown in Table I. 


DISCUSSION 


The chain structure found for dimethyltri- 
sulfide is in agreement with the chain structures 
proposed for the three disulfides which have 
been previously investigated by the electron 
diffraction method, H2S2,° S2Cle,® 7 and (CH3)2S2.° 
In all four molecules the bonded S—S distance 
is 0.02-0.03A shorter than the S—S distance 
2.07A in Ss.!4 The C—S distances in (CH3)2Se 


4 C. S. Lu and J. Donohue, J. Am. Chem, Soc. 66, 818 
(1944). 
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4.42 

6.61 

8.52 
13.35 
17.73 
21.59 
25.70 
29.37 
33.85 
37.86 
42.03 
45.93 
50.93 
56.19 
60.92 
65.17 
70.51 
75.78 
81.50 
86.12 
91.33 
96.74 





0.989 
998 
1.000 


0.988 
975 
985 
.987 

1.006 

0.999 
995 
998 

1,008 

1.001 

1.007 

0.996 
-996 
999 
999 


average 0.996 


average deviation 0.006 








and (CHs3)eS; are the same, and are 0.04A 
shorter than the distance in (CH3;)2S.'° These 
shortenings are probably significant although 
they are within the limits of experimental error 
in each case. When sulfur forms single covalent 
bonds with two other atoms, the existing data 
indicate the angle decreases as the electro- 
negativity of the bonded atoms increases: the 
observed values are 107° in (CH3)2Se, 105° in Sg, 
104° in (CH3)2Ss, 103° in SoCle, and 101° in SCls. 

In addition to the direct experimental evi- 
dence of the electron diffraction data, the sulfur- 
sulfur distance also suggests that the molecule 
has the straight chain structure I, inasmuch as 
it is approximately equal to the S—S distance 
in the singly bonded molecules just mentioned 
and is indeed just equal to the distance in 
(CH3)2Se, which is agreed to have the singly 
bonded structure; it is about 0.15A greater than 


%L. O. Brockway and H. O. Jenkins, J. Am.{Chem. 
Soc. 58, 2036 (1936). 


would be expected for an S—S bond of the type 
to be found in structure III and about 0.08A 
greater than the average to be expected for 
structure II, on the basis of the Pauling-Huggins 
covalent radii.'* Although no test of these co- 
valent-radius predictions of sulfur-sulfur multiple 
bonds is available, a rather close analogy can 
be drawn with the phosphorus-sulfur bonds in 
P,S3, in which the bonds are presumed to be 
single bonds, as contrasted to the phosphorus- 
sulfur bonds in phosphorus sulfoxide, thiophos- 
phory] fluoride, and thiophosphory] chloride. The 
P—S distances in P,O.S,, PSF3, and PSCl; are 
1.85A,!7 1.85A,18 and 1.94A,!® respectively (the 
covalent radius sum for a double bond is 1.94A),!° 
while in P,S; the P—S single bond distance is 
about 2.12A”° (radius sum 2.14A). Here the 0.2A 
difference between the lengths of the single and 
double bonds between phosphorus and sulfur is 
found to be valid, even though there might be 
some question, as there is for the special sulfur- 
sulfur bonds in structures II and III, as to 
whether the phosphorus-sulfur bond in the thio- 
phosphoryl compounds should have the same 
length as that predicted in a conventional double 
bond. 

The same conclusion is suggested by the bond 
angles of 104°. These bond angle values are 
reasonable for the singly bonded atoms in the 
chain structure I but they do not seem reasonable 
for the angles involving the doubly bonded 
sulfur atoms in structures II and III. For struc- 
ture III, moreover, the S—S—S and C—S-—S 
bond angles of 104° correspond to a C—S—C 
bond angle of greater than 130°, an unreasonably 
high value by a wide margin. 


16L. Pauling, The Nature of the Chemical Bond (Cornell 
University Press, Ithaca, New York, 1940), second edition, 
Chapter IV. 

17 A, J. Stosick, J. Am. Chem. Soc. 61, 1130 (1939). 

18D. P. Stevenson and H. Russell, J. Am. Chem. Soc. 
61, 3264 (1939). 

19 J. Y. Beach and D. P. Stevenson, J. Chem. Phys. 6, 
75 (1938). 

20Q. Hassel and A. Petterson, Tids. Kjem. Bergvesen 
Met. 1, 57 (1941); Chem. Abs. 35, 4646 (1941). 
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The characteristics of the faint Raman lines appearing in long exposures have been investi- 
gated for the molecules, SnBr4, SnCl4, and CCl,. All the lines observed have been classified as 
overtones and combination tones of the fundamentals. The depolarizations of the second-order 
lines were measured by means of a Wollaston prism placed between the collimator lens and 
the dispersing prism of the spectrograph. The depolarizations were foufhd to be in accordance 
with theory, on the assumption of tetrahedral symmetry for the molecules. The intensities 
of the second-order lines relative to the fundamentals were measured. For the spherically 
symmetrical frequencies it is shown that the observed intensities can be explained by assuming 
that electrical and mechanical anharmonicities are both operative; in some cases counter- 
acting one another and in other cases enhancing the effect of each other. 





INTRODUCTION 


TRONGLY exposed Raman spectrograms of 
many substances show very faint, or second- 
order, Raman lines in addition to the intense, or 
first-order, lines due to the Raman-active funda- 
mental vibrations of the scattering molecule. 
A systematic study of these lines of second-order 
intensity in a related group of molecules should 


aid in the understanding of the physical processes 
involved in the molecular scattering of light, and 
provide further information concerning the struc- 
ture of molecules. 

In general, Raman lines of very low intensity 
should be attributable to transitions that are 
forbidden by the first-order selection rules, but 
which are allowed when the anharmonicities 
present in the molecules are taken into account. 
Thus, overtones and combination tones of the 
fundamentals should appear, and possibly cer- 
tain fundamental tones which are forbidden in 
the first approximation by symmetry. In some 
cases second-order lines occur with unusually 
high intensity which can be associated with 
accidental resonance degeneracies in the mole- 
cule. Second-order lines can also arise from 
multiple Raman scattering. However, the in- 
tensity to be expected for this process is negligible 
compared with that of the second-order lines 
which have been observed. 

The characteristics of second-order lines that 


*The greater part of the experimental work reported 
here was carried out in 1940-1941 while one of the authors 
(G. D. Scott) was holder of a Studentship granted by the 
National Research Council of Canada. 
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are of primary importance are: (a) the frequency 
shifts, (b) the intensities relative to each other 
and to the first-order lines, (c) the degree of 
depolarization, (d) the intensity contours. The 
frequency shift of a line identifies the transition 
giving rise to it. The relative intensities afford 
information on the electrical and mechanical 
anharmonicities of the molecule, although the 
theory is not yet sufficiently advanced to make 
full use of these measurements. The depolariza- 
tion values can be used to test the polarizability 
theory and the symmetry assumed for the 
molecule. Finally, the intensity contours can be 
of considerable significance in cases in which 
accidental resonance degeneracies are present in 
the molecule. 

Although overtones and combination tones of 
fundamental vibrations are of common occur- 
rence in infra-red spectra, only relatively few 
have been observed in Raman spectra. This is 
probably caused by the experimental difficulties 
of observing these faint, second-order Raman 
lines. The anomalously strong overtone of the 
accidentally degenerate vibration, Avy= 1550 cm 
in CCl, has been studied by a number of 
workers.'? In addition to this overtone, two 
second-order lines in the spectrum of CCl, and 
one in the spectrum of SiCl, were observed by 
Ananthakrishnan.* The work on CCl, was ex- 


1L. S. Ornstein and J. J. Went, Proc. Amst. Akad. Sci. 
35, 1024 (1932). 
1932) Carrelli and J. J. Went, Zeits. f. Physik 76, 236 
wa R. Ananthakrishnan, Proc. Ind. Acad. Sci. 2, 425 
1935). 
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tended by Landsberg and Malysev;‘ they also 
investigated SnBr, and recorded seven second- 
order lines. Intensity and polarization measure- 
ments have been made for the overtone Av = 1550 
cm in CCl, by Ornstein,! and by Carrelli and 
Went.? Landsberg and Malysev estimated that 
the intensities of the second-order lines were 
about 3X 10-* of the intensity of the first-order 
lines. No further intensity or polarization meas- 
urements of second-order lines have been pub- 
lished. 
EXPERIMENTAL 


The Spectrograph 


A double prism glass spectrograph was used. 
The refracting components are a constant devia- 
tion flint glass prism and a 60° crown glass prism. 
The collimator lens is an f:10 doublet of 75-cm 
focal length. Either of two camera lenses can be 
fitted, an f:4 triplet of 30-cm focal length or an 
f:10 doublet of 75-cm focal length. Both lenses 
give excellent definition over the Raman region 
of the spectrum. With the f:4 lens the dispersion 
is 43A/mm (225 cm/mm) at 4358A; with the 
f:10 lens, it is 16A/mm (86 cm-!/mm). 


Raman Tubes 


The Raman tubes were of conventional design 


,/ | with a plane window at one end and a blackened 


“horn at the other. The standard type tube had 
a diameter of 1 cm and an effective length of 
10 cm, and was mounted directly in front of the 
slit of the spectrograph. A tube 100 cm in length 
and 1.6 cm in diameter was also used, in this 
case with a condensing lens. With the same 
illumination per unit length of the tube, the 
brightness of the Raman spectrum from the 100- 


4G. S. Landsberg and V. J. Malysev, C. R. (Doklady) 
Acad. Sci. U.S.S.R. 3. 365 (1936). 
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Fic. 1. Raman tube used with 
condensing lens. 
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cm tube was approximately tenfold that from 
the standard tube. Both types of tubes had 
jackets through which a saturated solution of 
sodium nitrite was circulated. The solution served 
as a light filter and a temperature control. Most 
of the exposures were taken with the scattering 
substance at 25°C. 

When a long Raman tube is used, a condensing 
lens placed between the window of the tube and 
the slit of the spectrograph is essential. It can be 
shown that for a given diameter the usable 
length of the tube is a maximum when the image 
of the window of the tube, formed by the con- 
densing lens, is at the collimator lens and equal 
in size to its effective aperture and the image of 
the horn end of the tube is at the slit with a 
diameter equal to the slit length. Figure 1 
illustrates this arrangement for a cylindrical 
Raman tube. The dimensions of the tube under 
these conditions are limited by the relation, 
D*?=Los/n, where D is the diameter of the tube, 
L is its irradiated length, s is the slit length, 
o is the reciprocal of the stop number of the 
collimator lens, and m is the index of refraction 
of the scattering medium. The focal length of the 
condensing lens is, F=D/o(1—s/cf), where f is 
the focal length of the collimator lens. J. R. 
Nielsen® derived equivalent relations in 1930 and 
pointed out the application to Raman spec- 
troscopy. 


Mercury Lamps 


The irradiating sources were low pressure 
mercury vapor lamps operated at 2.5 to 3.0 
amperes from a 220-volt d.c. supply. They were 
constructed of Pyrex tubing of 9-mm bore and 
had an extra tungsten electrode near the anode 
to facilitate starting with an induction coil. In 


5 J. R. Nielsen, J. Opt. Soc. Am. 20, 701 (1930). 
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this way, lamps with a discharge path up to 
120 cm could be started on a 220-volt supply. 

Lamps of different forms were used: straight 
lamps for polarization, and either coil or straight 
for ordinary spectra. In the latter case irradia- 
tion of the tube was always arranged so that 
the relative intensities of lines of different degrees 
of depolarization were not dependent on reflec- 
tion losses at the faces of the dispersing prisms. 

When the 100-cm tube was irradiated with a 
bank of eight lamps, the first-order Raman spec- 
trum of CCl, was recorded in 5 seconds. 


Depolarization Determinations 


The essential features of the polarization tech- 
nique used here are directional illumination of 
the Raman tube and division of the beam 
emerging from the collimator into two separated 
beams polarized in perpendicular directions. The 
mercury lamps are arranged so that the un- 
polarized light illuminating the Raman tube is 
confined approximately to a plane through the 
axis of the tube. The beam from the collimator 
is divided by a large Wollaston double image 
prism of quartz placed between the collimator 
and the dispersing prism. The double image 
prism, which can be rotated in its mount about 
an axis coincident with the axis of the collimator, 
is so oriented that the polarization of one of 
the emergent beams is parallel and that of the 
other perpendicular to the plane of irradiation. 
With this arrangement two adjacent spectra are 
recorded simultaneously, one in each of the two 
polarizations required to evaluate the depolariza- 
tion factor p. Differential reflectién of the two 
polarized beams at the faces of the dispersing 
prisms can be eliminated by setting the plane of 
irradiation at 45° to the plane of incidence at 
the dispersing prisms. 

When the illumination of the Raman tube is 
not unidirectional, the ratio of the intensities 
of the same line in the two spectra, which is 
called herein the intensity ratio of the L to the || 
component, gives only an approximate value 
of p. An accurate value of p can be obtained by 
approximating closely to the ideal arrangement 
in which the incident radiation is unidirectional 
in the scattering medium, but this necessitates 
much longer exposures. 





An alternative method is to calibrate the 
experimental arrangement in which the irradia- 
tion is only approximately unidirectional by 
measuring the intensity ratio for a number of 
Raman lines whose p values have been deter- 
mined accurately. This method does not require 
elimination of differential reflection of the two 
polarized beams because the calibration auto- 
matically corrects for the difference in reflection 
losses. In fact there is some advantage in setting 
the plane of irradiation parallel to the plane of 
incidence at the dispersing prisms. The intensity 
of the _ component is then increased relative to 
that of the || component, and for a highly 
polarized line their ratio can be measured more 
accurately. Nor is it essential, if the arrangement 
is calibrated, to orientate the double image 
prism so that the polarization of one of the 
emergent beams is parallel to the plane of irradi- 
ation. However, since maximum sensitivity is 
obtained when the double image prism is so 
orientated, it is not advisable to depart far from 
this optimum orientation. 


Hg4339 Hg Hg 4358 
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Fic. 2. Raman spectra of SnBr, excited by Hg 4358; 
“spectral” width of slit 3 wave numbers: (a) 1-minute ex- 
posure showing lines due to the four fundamental vibra- 
tions of the molecule; (b) 6-hour exposure showing the 
higher frequency overtones and combination tones. 

















































































TABLE I. Raman lines of SnBry, liquid. 


WELSH, CRAWFORD, AND SCOTT 











Frequency shift 
(cm-1) 
Meas- Calcu- Peak 
ured lated intensity Remarks 
v3— V2 — 20 covered by wings 
Munn — 59 of Rayleigh line 
v2 68 370 merged 
Vs 88 ‘ 560 
vi-V3 133 * 
vi V2 1 * 
ate 151 156 1.5 
23 176 176 0.7 . 
Va v3 193 192 . 6.5 53 
V4—V2 — 213 * 
v1 221 1000 sharp, symmetrical 
v4 280 180 broad, degraded to red 
vitve — - 289 6 
rites 303 309 “= * 
vets —— 348 
vats 367 3=— 368 1.5 * broad 
2 442 442 2.5 sharp, symmetrical 
vite, 502 S501 0.9 degraded to red 
2 565 560 0.1 broad 








* These lines are on the wings of fundamentals, merged with each 
other or merged with the fundamentals from Hg 4348 or Hg 4339; the 
intensities of such lines are therefore very approximate. 


Since highly unidirectional irradiation would 
necessitate extremely long exposures for second- 
order lines, the alternative method was used in 
these experiments. The plane of irradiation was 
inclined 45° to the plane of incidence and the 
double image prism was set 8° from the optimum 
orientation. This orientation of the Wollaston 
prism makes the | and || components of a 
depolarized line equal in intensity and facilitates 
the visual identification of depolarized lines on 
the spectrogram. The arrangement was calibrated 
using the Avy=459(p=.05) and Av=218(p=.83) 
lines of CCly. The measured intensity ratios for 
these lines were 0.33 and 1.0, respectively. These 
values are sufficient to determine the calibration 
curve since previous calibrations made with a 
more complete range of p values showed that the 
relation between the measured intensity ratios 
and p values is almost linear. 


Intensity Measurements 


A comparison of the intensities of first- and 
second-order Raman lines requires a technique 
capable of measuring, at different frequencies, 
intensities that differ by a factor of 2X10‘. 
Stepped slit spectra of a Kipp and Zonen stand- 
ard lamp were used as intensity standards to 
calibrate the photographic plates. The spectral 





distribution of energy emitted by this lamp for a 
wide range of filament currents is known from 
calibration data supplied by the manufacturer. 
The desired range of intensities is readily ob- 
tained with this method since a 20-fold range of 
intensities is produced by the stepped slit and a 
5000-fold change can be produced by changing 
the current through the lamp. 

Ilford Special Lantern plates were used. Pre- 
vious intensity work showed that these plates 
had adequate sensitivity and fineness of grain, 
and that plates of the same emulsion number 
had very uniform sensitivity characteristics. To 
increase their sensitivity for lines of low density 
the plates were prefogged in green light to a 
density of 0.1. The prefogging was, in general, 
advantageous, since the continuous background 
of the mercury lamps was extremely low. The 
backing on the plates as supplied was not suffi- 
ciently effective for very strong exposures and 
was replaced by Eastman Opaque which elimi- 
nated halation rings. 

To eliminate errors due to reciprocity failure, 
the Raman spectrum and the intensity standards 
are usually exposed for equal times. The incon- 
venience of this procedure for long exposure 
times was obviated by photographing the stepped 
slit spectrum on each plate with an exposure 
time of 15 minutes and then making corrections 


TABLE IJ. Raman lines of SnCl, liquid. 








F wae shift 





(cm7 
Meas- Calcu- Peak 
ured lated intensity Remarks 

ve—V2 — 25 covered by wings 
vu" — 35 of Rayleigh line 

v2 106 475 merged 

v3 131 520 

2v2 213 = -- 
Vi-— V3 23 — 
vatns, 241 937 

2v3 — 262 covered by 71 
vi— ve — 262 from Hg 4339 
va 3 — 272 

covered by »; from 
va—v2 — 297 Hg 4348 and by 
from Hg 4339 

v1 368 1000 sharp, symmetrical 

% 403 265 broad, degraded to red 
v 24 ve 475 = 0.9 broad 
vitvs 
pil ing 502 509 0.4 
vatr, 534 534 0.8 broad 

2 737 =—736 0.6 sharp, symmetrical 
m+y 773 771 0.5 broad, degraded to red 

2% 814? 806 0.05? very faint 
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for reciprocity failure by means of a “‘master”’ 
intensity calibration for plates of the same 
emulsion. number. Calibrations at three wave- 
lengths, 4350A, 4500A, and 4750A, were ob- 
tained by determining the characteristic curves 
for exposure times of 1, 5, 15, 60, and 240 
minutes. From these the relative intensity of a 
Raman line of any frequency shift up to 2000 
cm~ photographed with any exposure time from 
30 seconds to 10 hours could be readily deter- 
mined by interpolation or extrapolation. 

A Moll recording microphotometer and later a 
Leeds and Northrup recording microphotometer 
were used to measure the photographic densities. 


RELATIVE INTENSITIES 


Commercial c.p. samples of CCl, and SnCl, 
were used. The SnBr, was prepared from mossy 
tin and bromine. Each substance was care- 
fully purified by preliminary distillations under 
vacuum, and then distilled into the Raman tube. 
In one sample of CCl, one of the combination 
tones of CCl, was obscured by a stronger, highly 
polarized line; it was identified as the Av=667 
cm line of chloroform. No lines which could 
be ascribed to impurities were found in any of 
the other samples of the three substances. 

A number of plates were taken for each sub- 
stance and at least two exposures were made on 
each plate: a short exposure, 1 to 10 minutes, 
to record the fundamental frequencies with suit- 
able densities and one or more long exposures, 
30 minutes to 12 hours, to record the second- 
order lines (see Fig. 2). The “‘spectral’’ width of 
the slit was from 3 cm= to 15 cm. With these 
widths approximate intensity contours could be 
obtained for all lines. 

The frequencies and relative intensities of the 
observed Stokes Raman lines of SnBr4, SnCl,, 
and CCl, are given in Tables I, II, and III, 
respectively. On the assumption of tetrahedral 
symmetry these molecules have four fundamental 
modes of vibration all of which are active in 
Raman effect. The vibrations are designated as 
v, for the symmetrical frequency, v2 for the 
doubly degenerate frequency, v3; and v4 for the 
triply degenerate frequencies. All second-order 
lines observed have been identified as overtones 
and combination tones and are designated as 
such in terms of the fundamentals, e.g., 211, 
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TABLE III. Raman lines of CCl, liquid. 











Frequency shift 
(cm) 
Meas- Calcu- Peak 
ured lated intensity Remarks 
= ae Bs ‘as em } on wing of Rayleigh line 
v2 218 890 symmetrical 
vive _ 241 
va 314 915 symmetrical, with faint wings 
2v2 _ 436 
v1 459 1000 asymmetrical, faint wings 
ve+vs 535 532 0.7 overlaps v4’ —v2 
va —v2 — 544 
va! —v2 580 572 -- broad, faint 
2v3 630 628 0.7 
rt be 673 677 0.4 overlaps v4’ from Hg 4339 
- o a4 = } symmetrical doublet 
21 915 918 A 
vetvd’ 989 980 0.1 broau, diffuse 
votre” 1008 sesait teutied, suse 
, roa oublet, higher 
a i, et ae os frequency component 
piv 3 | brad double high 
t higher 
n+te"d 1221 1221 0.3 | th ay 
equency component 
vive’ 1250 1249 0.3 slightly sharper 


1524 
2 1540 1552 17 broad, asymmetrical 
1580 
vo+2 1749 1758 0.03 very broad 
vat2v 1860 1854 0.03 very broad 
vit2 1995 1999 0.04 very broad 











vit+v2. The measured shifts for the fundamentals 
have been used to obtain the ‘‘calculated”’ fre- 
quency shifts listed for the overtones and com- 
bination tones. The frequency measurements are 
accurate to +1 cm for strong lines and to 
+5 cm for faint or broad lines. The relative 
intensities are given on a scale fixed by taking 
the intensity of »; as 1000. The error in the 
intensity measurements is about +5 percent for 
the fundamentals and +20 percent for the 
second-order lines when there is no overlapping. 
In most cases the second-order lines were found 
to have half-value widths approximately equal 
to the sums of the half-value widths of the corre- 
sponding fundamentals. 

The intensities of the second-order lines rela- 
tive to the first-order lines vary widely as may 
be seen in Tables I, II, and III. On the average, 
however, the ratio in CCl, (excluding 24) and 
in SnCl, is about 1/1000; but in SnBr, it is larger, 
about 3/1000. 

The relative intensities of the first-order lines 
of the three substances were measured in binary 
mixtures of known concentration. For the totally 
symmetrical frequency 1, of CCly, SnCly, and 
SnBr, the peak intensities per molecule were 
found to be in the ratio 1:6.7:30. This ratio 




















































































































































































































































TABLE IV. Depolarizations of overtones and 
combination tones. 











Depolar- Depolar- 
Measured ization ization 
intensity classifi- classifi- 
ratio cation Symmetry type of cation 
L/il (experi- upper vibrational (theo- 
components mental) state retical) 
SnBr, 
vats 1.0 dp A:+E+F4+ Fe pp 
2, 0.3 p Ai p 
vit 1.0 dp F. dp 
2% 0.7 pp Ai tE+F, pp 
CCl, 
V3— V2 0.8 dp Fi+F,2 dp 
vi— v3 0.9 dp F, dp 
votps 1.3 dp Fi+ Fs dp 
2v3 0.8 dp Ai +£E+F2 bp 
vitve 1.1 dp E dp 
2% 0.5 pp Ai+tE+F; pp 








correlates the intensities in the three tables. 
Since the half-value widths of the »; lines for 
CCl,, SnCl4, and SnBr,y are about 12 cm™, 
9 cm=, and 5.5 cm™, respectively, the integrated 
intensities per molecule for this vibration are in 
the ratio 1:5:14. Thus, as anticipated due to the 
increase of polarizability with atomic number, 
the: scattering per molecule is larger when an 
atom in the molecule is replaced by a heavier 
atom of the same chemical group. 

It is interesting to note that instances of triple 
quantum jumps were observed, namely, the lines 
vit2v4, vot2v4, and v3+2v4 of CCly. These lines 
are extremely faint and broad. Their appearance 
and the high intensity of 2, indicate a strong 
anharmonicity in the v4 mode.** 


DEPOLARIZATIONS 


Polarization spectra were obtained for SnBr, 
and CCl. The intensity ratios of the L to the || 
component for the overtones and combination 
tones that could be measured are given in 
Table IV. The lines are classified as completely 
polarized, if this ratio is 0.3 to 0.4, depolarized 
if it is 0.8 to 1.2, and partially polarized if it is 
between 0.4 and 0.8. The tolerances for the first 
two classifications are set mainly by the experi- 
mental error in the intensity ratios and not by 


** The intensity contour of 2», in CCly, especially in the 
polarization spectra, indicates three overlapping com- 
ponents with different degrees of depolarization. A detailed 
discussion of this and other features of the CCl, spectrum 
in relation to the resonance degeneracy will be presented 
elsewhere. 
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the error in the polarization calibration. The 
types of vibration and the depolarizations ex- 
pected on the basis of tetrahedral symmetry® are 
shown in columns four and five. The agreement 
between theory and experiment is complete for 
all the lines measured except for the v3+ 4 line 
of SnBr, and the 27; line of CCl,. Even for these 
two the agreement can be considered satisfactory 
since the experimental classification of a line as 
dp shows only that p has a value near 6/7 and 
the theoretical classification of a line as pp means 
that p can have any value from 0 to 6/7. 

The results verify such predictions as: (1) the 
overtone of the polarized symmetrical frequency 
v1 is polarized ; (2) the overtone of the completely 
depolarized triply degenerate frequency v is 
partially polarized ; (3) the binary combinations 
involving the polarized frequency » with a 
completely depolarized frequency v2, v3, or v4 are 
completely depolarized. 


DISCUSSION OF THE INTENSITIES OF 
SECOND-ORDER LINES 


A detailed interpretation of the intensities 
requires some of the results of the polarizability 
theory of scattering.’ The intensity of the radia- 
tion scattered by a molecule is proportional to 
the square of the matrix elements of the polariza- 
bility tensor, a, that is to the square of (a), 
where v and v’ are the quantum numbers of the 
initial and final states. a can be expanded in the 
usual way in terms of the normal coordinates qj. 
The matrix elements in the derivative terms in 
the expansion of a give rise to Raman scattering. 
The intensity of a Raman line of frequency 
shift v;, corresponding to the ith normal mode 
of vibration of the molecule is proportional to 
the product of {da/dq;(q;)xi:+1°*}? and the appro- 
priate Boltzmann population factor, summed 
over all values of v;. The intensity of the line is 
also proportional to the fourth power of its 
frequency. 

Transitions corresponding to overtones and 
combination tones can occur owing to electrical 
anharmonicities or to mechanical anharmonicities 
or to both. When electrical anharmonicities are 
present, the derivatives in the expansion of the 
polarizability beyond the first are different from 


6 L. Tisza, Zeits. f. Physik 82, 48 (1933). 
7G. Placzek, Hand. d. Radiologie, VI-2 (1934). 
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zero. When mechanical anharmonicities are pres- 
ent, the classical motion of the nuclei is not 
simply harmonic and hence the matrix element 
(qi)»’” is different from zero for transitions other 
than v;—v;+1, that is, for example, for v;-v;+2 
and for vj, 7%.—v;+1, v,+1. Thus, if both types 
of anharmonicity are present the intensity of the 
first overtone of the ith vibration depends on 
the square of matrix elements of the form: 


2 


al Oa ee eae OP (1) 
ae pal ae 
qi sete 2 dq? 


Di+2 


Similarly, the intensity of a combination tone 
involving the jth and kth vibrations will depend 
on the square of matrix elements of the form: 


0a 0a Djs 
(EeEs) 
0q; Oge 7 vj+1, %+1 
07a 
of. 





Vj, Ve 


Of the two terms in each of the expressions, the 
first is the contribution of mechanical anhar- 
monicity, the second the contribution of elec- 
trical anharmonicity. It is possible that under 
certain conditions the terms due to the two 
anharmonicities can be of opposite sign and 
hence combine to give a value of the intensity 
of the second-order line smaller than that ex- 
pected for either anharmonicity alone. 

The relative intensities of overtones and com- 
bination tones cannot, in general, be calculated 
directly with present theory. Placzek’ derived 
general expressions for the intensities of over- 
tones assuming only electrical anharmonicity; 
however, the derivatives of the polarizability in 
his expressions are undetermined. Wolkenstein® 
has recently developed a classical, semi-empirical 
theory based on the additivity of bond polariza- 
bilities for calculating intensities in Raman 
spectra of polyatomic molecules. This theory is 
useful in dealing with the relative intensities of 
fundamental vibrations, but the expressions de- 
rived for the intensities of overtones and combina- 
tion tones do not admit experimental verification. 

In the special case of the spherically sym- 
metrical vibration of an XY, molecule, it is 


*M. W. Wolkenstein, J. Phys. U.S.S.R. 5, 185 (1941). 
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possible to estimate the term due to the mechani- 
cal anharmonicity in the matrix element (Eq. 
(1)). Since only changes in the X — Y bond length 
occur and the vibration is symmetrical, the 
potential energy can be represented by a Morse 
function and the known wave functions for this 
potential can be used in calculating (q),+2”. 

The potential energy takes the form, 


Vir) = 4D(1 a e~alr—re))2 (3) 


in which D is the X— Y bond energy and a isa 
constant. The change in the length of the bonds 
is (r—r.) which can be considered as the normal 
coordinate g. The vibrational energy is expressible 
in the usual form, 


E,=we(0+}) —xwe(v+4)’. (4) 


If the dissociation energy, 4D, of the molecule 
and the vibrational frequency, w., are known, 
the anharmonicity constant, x., can be estimated 
from the formula: 


X-=w,/16D. (5) 


The matrix elements (g),41” and (q).42” can be 
obtained as functions of x, by using the method 
given by Dunham? for the calculation of the 
intensities of overtones in infra-red spectra. For 
the 0—>1 and 0-2 transitions these are: 


(g)1°=1/a(k—3)} 


; : 6 
x| ~«(-—-)- o«(—)} (©) 


TABLE V. Intensity ratio of the overtone to the fundamental 
for the symmetrical vibration in XY, molecules. 














Io 





I 

(Calculated on the Io 

Anharmonicity basis of mechanical —_ 

constant anharmonicity alone Is 
Xe by Eq. (13)) (observed) 
SnBry 0.00060 1/2900 1/400 
SnCl, 0.00085 1/2200 1/1700 
CCl, 0.00111 1/1800 1/10,000 





TABLE VI. a”’/a’X10-§ cm. 








SnBr, 0.78 or —1.55 
SnCk 0.098 or —0.90 
CCl, —0.29 or —0.75 











®J. L. Dunham, Phys. Rev. 34, 438 (1929). 
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and 


on ve) { ssines e(3) 


—2(h—3) loe(- =) 


. )| (7) 
k—34 
in which k=1/x,.*** 


The application of these formulae is consider- 
ably simplified by expanding them in series 
and introducing approximations justifiable when 
k(=1/x,) is >1. With this procedure (6) and (7) 
reduce to: 


(g)1°=1/a(1/k4) =1/a(x,}) =, (8) 
(q)2°=1/a(1/24k) =1/a(x./2*) =b(x./2)*. (9) 


Since for the Morse function a= (82’ucw,x./h)}, 
the constant 6=x,!/a, introduced in the above 
expressions is merely the zero point amplitude 
of the simple harmonic oscillator of the same 
frequency. Thus, for this approximation, the 
value of (q):° is the same as that for the simple 
harmonic oscillator. 

Since the v=1 level of the molecules studied 
is also populated to some extent at a temperature 
of 25°C, the 1—+2 and 1—3 matrix elements are 
also required. A calculation similar to the above 
gives: 














+(k—4) loe( 


(g)2'=5-2?, - 
(q)3' = b(3x,/2)!. 


When the ground vibrational state only is 
considered as the initial state and when the 
frequency factor is neglected, the ratio of the 
intensity of the overtone due to mechanical 
anharmonicity alone to the intensity of the 
fundamental can be expressed very simply in 
terms of the anharmonicity constant: 

Io {a'(g)2°}*? xe (12) 
Ty {a’(g):°}* 2 


where a’ is written for da/dq. 


*** S. Bhagavantam (Ind. Acad. Sci. 10A, 224 (1939)) has 
used these formulae to calculate Raman overtone intensi- 
ties in diatomic molecules. 


When the transitions from the v=1 state are 
also taken into account and the frequency factor 
is included, the ratio becomes: 
In Xe f1+3e-ee!*Ty sy—2u.\4 
a roow Bn 
Tp 2N142e-we/k? V—We 

in which » is the frequency of the exciting light. 

Thus the determination of this intensity ratio 
requires only the value of x, which can be esti- 
mated from the dissociation energy by means 
of (5). 

In Table V are listed the calculated values of 
x, and J)/Z;, and also the experimental value of 
I)/I;. The dissociation energies used to obtain x, 
were taken from the tables given by Samuel." It 
is seen that there is no agreement between calcu- 
lated and experimental results. However, this is 
not unexpected, as it is unlikely that mechanical 
anharmonicity alone is responsible for the appear- 
ance of overtones. From a comparison of the 
observed and calculated intensity ratios it is 
evident that in the case of CCl, the electrical 
and mechanical anharmonicities counteract each 
other, whereas in the case of SnBr, they augment 
each other. These two possibilities are admitted 
by the theory as pointed out in the discussion 
of expression (1). 

The experimental values of the intensity ratios 
and the calculated values of the terms of the 
matrix elements of the polarizability due to 
mechanical anharmonicity can be used to obtain 
values for a’’/a’, the ratio of the second and first 
derivatives of the polarizability in the neighbor- 
hood of the equilibrium position. If for simplicity 
only transitions out of the ground state are 
considered, the intensities of the overtone and 
the fundamental can be given as: 


To=c(v—2w)*{(a’)ose+3(a’)os2}? (14) 





and 


I;=c(v—w)*{ (a’)o41}? (15) 


in which ¢ is a constant factor proportional to 
the intensity of the exciting light; (a’)o.2 is 
written for a’(q)2°, and (a’’)o.2 for a’’(g)2°. Sub- 
stituting from expressions (8) and (9) in (a’)o 
and (a’)o.2, and putting (e’’)o,2.=0%a”’ in the 
harmonic oscillator approximation, the intensity 


10 R, Samuel, Rev. Mod. Phys. 18, 103 (1946). 
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ratio becomes: 








To (—) eae 
iy V—We ba’ 

y—2w- 4 a’)? 
-(—) {(s./2)+-1/20—| | 
Qa 


V— We 


(16) 


Writing the frequency factor as ** and rearranging 
the equation we have: 

“) 

I; 


+ 
Xe\? 
5 (5) 
It is seen that the experimental value of Jo/J; 
and the calculated value of x, lead to two values 
for a’’/a’. These are given in Table VI. If the 
excitation of higher vibrational levels is taken 


into account, these results are only slightly 
altered. 


(17) 
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The absolute values of a’ and a” are not 
obtainable from these measurements alone. How- 
ever, if the ratio of the intensities of the Raman 
line and the corresponding Rayleigh line were 
available, the value of a’/a could be found. 
Then, if a were known, say, from refractive index 
measurements the values of a’ and a” could be 
derived. 

In conclusion, it may be said that the investi- 
gation presented here shows that overtones and 
combination tones in Raman scattering depend 
on both the electrical and mechanical anhar- 
monicities of the molecule. In special cases the 
intensity of the second-order line relative to the 
first-order line can be used to obtain the second 
derivative of the bond polarizability if the anhar- 
monicity constant of the vibration is known, and 
if the first derivative of the polarizability has 
been obtained from the intensity ratio of Raman 
to Rayleigh scattering. Experimental data of this 
kind should be of use when the theory of the 
polarizability of molecules is further developed. 
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Ionization Potentials of Some Hydrocarbon Series 


RicHARD E, HoNniG 
Research and Development Department, Socony-Vacuum Laboratories, Paulsboro, New Jersey 
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A new method is described which permits ionization and 
appearance potentials to be determined reproducibly, 
within +0.02 volt, by electron impact in a conventional 
180° mass spectrometer. The energy spread of the ionizing 
electrons, while not eliminated, is accounted for in an 
analytical expression that describes the ionization efficiency 
curve of argon near the initial upward break. It is shown 
that an ionization or appearance potential can be deter- 
mined graphically by drawing a tangent of critical slope 
to the efficiency curve which has been plotted on a semilog 
scale. Absolute values are obtained by referring the un- 


I. INTRODUCTION 


ONIZATION potentials of organic molecules 
are not only of theoretical interest, but may 
also be applied to practical problems, such as 
the identification of isomers. A large volume of 
work has been published on critical potentials of 
individual molecules by electron impact? and 
' For recent work on hydrocarbons see: D. P. Stevenson 


and J. A. Hipple, J. Am. Chem. Soc. 64, 1588, 2766, 2769 
(1942); 65, 209 (1943). 


known to the ionization potential of a calibrating gas 
(usually krypton). For the monatomic gases studied, 
electron impact and spectroscopic values agree within 
the experimental uncertainty. For some known poly- 
atomic molecules, the ‘‘vertical’’ impact values exceed the 
“adiabatic” spectroscopic potentials by 0.02 to 0.2 volt, 
depending on molecular complexity. Paraffins, 1-olefins, 
and 2-olefins show a systematic reduction in potential 
with chain length. A similar decrease is observed when 
the mono-olefins are arranged according to the number and 
size of substitutions on the basic ethylene structure. 


by Rydberg series in the ultraviolet,* but, with 
one exception,‘ no systematic survey of any 
particular series has been made so far. The 


( 947) H. Dibeler, J. Research Nat. Bur. Stand. 38, 329 
1947). 

’For work by W. C. Price, D. M. Simpson, T. M. 
Sugden, W. T. Tutte, and A. D. Walsh see: Nature 148, 
373 (1941); Proc. Roy. Soc. A174, 207, 220 (1940); 179, 
201 (1941); 185, 176 (1945); Trans. Faraday Soc, 41, 35, 
76, 381, 498 (1945). 

*R. S. Mulliken, Rev. Mod. Phys. 14, 265 (1942); see 
also W. C. Price, Chem. Rev. 41, 257 (1947). 
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IONIZING VOLTAGE (UNCORRECTED) 

Fic. 1. Lower portion of a typical ionization efficiency 
curve (argon). tion a is exponential in character, 
section } is intermediate, and section ¢ is essentially 


straight. Point A denotes the initial upward break, while 
point B is located at the extrapolated straight-line intercept. 


present work is a mass spectrometric study of the 
variation of ionization potential with carbon 
number and molecular structure. 

In any electron impact method employing a 
hot filament, the spread of electron energies 
makes the exact ‘experimental determination of 
potentials difficult. Only for the Hg ionization 
potential have attempts been made to eliminate 
the energy spread by a magnetic filter.* In mass 
spectrometric work where such filtering would 
meet with practical difficulties, various methods 
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1ONIZING VOLTAGE (UNCORRECTED) 


Fic. 2. Semilog plot of a typical ionization efficiency 
curve (argon). C denotes the empirical critical slope. 


5 W. B. Nottingham, Phys. Rev. 55, 203 (1939). 
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have been devised to permit exact determinations 
without removing the cause of the spread. 

In his study of the CH, molecule,’ L. G. 
Smith used the initial ‘upward break”’ of the 
efficiency curve (point A in Fig. 1) to comnpare 
the potentials in question with a reference poten- 
tial ‘‘under conditions of similar sensitivities” by 
adjusting the ordinate scales of all curves until 
their straight portions were parallel. Since differ- 
ent processes have different absolute cross sec- 
tions, there is no theoretical justification for 
adjusting sensitivities in this manner, and con- 
siderable personal factors are involved in deciding 
where the curve first breaks away. However, 
Stevenson and Hipple showed in their comparison 
of argon and neon potentials’ that this method 
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Fic. 3. Comparison of theoretical expression (using m=2 
in Eq. (4)) with experimental values for the ionization 
efficiency curve of argon. Agreement is satisfactory to 
about 0.8 volt above the ionization potential. 





yields values that are within 0.15 volt of the 
spectroscopic difference, whereas the extrapola- 
tion of the straight-line portion (point B in 
Fig. 1) will yield erroneous results. The initial 
onset method has been used almost exclusively 
ever since. One of the few exceptions is found in 
a paper on SiCl, by Vought.* He attempts to 
show that the appearance potential curves of 
molecular fragments have a shape that can be 
explained by the transition probability from the 
ground state to the antibonding branch of the 
ionized state, and that the extrapolated straight- 
line intercept should be used. This reasoning 
‘L. G. Smith, Phys. Rev. 51, 263 (1937). 2 
7D. P. Stevenson and J. A. Hipple, Phys. Rev. 62, 23/ 


(1942). 
8 R. H. Vought, Phys. Rev. 71, 93 (1947). 
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IONIZATION POTENTIALS OF HYDROCARBONS 


does not take into account the formation of 
molecular fragments by the transition to the 
bonding branch of the ionized state. It extrapo- 
lates from the simple two-dimensional case of He 


to the unknown potential surfaces of a complex — 


molecule, such as SiCl4, and it is based on the 
simple but doubtful assumption that the proba- 
bility of transition is a linear function of the 
excess energy above the required minimum. 

This paper will present a new approach to the 
evaluation of ionization and appearance poten- 
tials, check the experimental accuracy obtained 
against the well-known ionization potentials of 
rare gases, and apply the new method to organic 
molecules. 
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_Fic. 4. Schematic diagram of electron and ion gun 
circuits of mass spectrometer. F—filament, S;S:—electron 
accelerating plates, B—back wall, ionization chamber, 
C—catcher, P;P2—pushers, S;S;—ion accelerating plates, 
and E—battery used for elimination of secondaries or 
(in reverse) for retarding potential measurements. 


Il, THEORY* 


If the ionization efficiency curve shown in 
Fig. 1 is redrawn on a semilog scale (Fig. 2), 
the exponential character of the lower portion 
becomes evident. Since the energy distribution 
of electrons emitted from a hot filament can be 
described by an exponential expression, it is of 
interest to investigate how the shape of an 
efficiency curve would be modified by a given 
distribution. An expression will be derived that 
is essentially exponential below the critical energy 
V., but reflects the shape of the original efficiency 
curve above V,. 

Assuming a quasi-Maxwellian energy distri- 
bution and a transmission coefficient of unity 

* The writer is greatly indebted to G. H. Wannier for a 


detailed derivation of the analytical expression presented 
in this section. 
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Fie. 5. Variation of the uncorrected ionization 
potential of argon with cross-field. 


(i.e., all electrons with thermal energies U larger 
than the work function W leave the filament), 
the number of electrons of energy U leaving the 
filament per second, 


4nrmA 


W-—n+U 
N.(U)= | (1) 


si exp aa 


h3 


where m=electronic mass in g, A =surface area 
in cm’, h=Planck’s constant=6.62X10-?’ erg 
sec., U=thermal energy in ergs, W=work func- 
tion in ergs, 7=Fermi limit in ergs, k= Boltz- 
mann constant = 1.38 X 10~"* erg/°K, and T =ab- 
solute temperature in °K. This expression, when 
integrated over all energies, yields Richardson’s 
equation. 

Before reaching the ionizing region, the elec- 
trons of thermal energy U pass through an 
electric accelerating field so that their total 
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Fic. 6. Typical retarding potential plot: electron current 
to catcher (J,) vs. ionizing voltage (V2—). The slope of the 
exponential straight line yields a distribution temperature 
of 2320°K. Magnet current J,,=0.4 amp. 





energy 
E=U+Y, (2) 


where V=energy associated with the acceler- 
ating voltage V/e. The total number of ions 
NV) produced per second with an electron 
accelerating voltage V/e is given by 


N(V)=] N(U)p(E)dU, (3) 


U=0 


where (EZ) = probability that an electron of total 
energy £ produces an ion which reaches the 
collector. Until now, this probability has not 
been determined theoretically near the critical 
energy V.. If it is permissible to use the empirical 
. relation 


P(E)=0 for EX V., (4) 
p(E)=Ci(E-V.)" for E2 V., 


then Eq. (3) can be evaluated by substituting 
Eqs. (1) and (4). To date, a linear probability 
function (7=1) has been generally used.® In the 
present study it was found empirically (see 
below) that the best agreement between experi- 
mental and computed values is obtained for 
n=2, yielding 


NAV) =2CkT*((V.— V)+3kT ] 
W-nt+V.- -| 
kT 





xX exo| _ 


for V< Vz, 
NV) = CT? 6k?T?+4kT(V — V.) 





+(V-V.)"] ~—| 
—. a 


for V2 Ve. 


If now In(J;) is plotted as a function of V, the 
curve approaches an exponential straight line 
below V, with a final slope 


© (nny) 6) 
or ae { 


For the critical energy V= V,, the critical slope 


: InN : 7 
— ary (7) 
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Thus, if all the simplifying assumptions made 
above were valid, the ionization potential or 
critical energy could be determined by applying 
a tangent of critical slope 2/3kT to the semilog 
plot. Modifications of these assumptions will 
change the absolute value of the critical slope, 
but it will always be a certain fixed fraction of 
the final slope. For the most general case (n= n 
in Eq. (4)), the critical slope can be shown to 


occur at 
eee (8) 
—(InN,).=——_ —. 
av n+1 kT 


In order to test Eq. (5) derived above, a 
typical ionization efficiency curve obtained ex- 
perimentally for argon was compared to the 
theoretical curve computed from Eq. (5). An 
electron distribution temperature T=2340°K 
was used (see Section III), and the two curves 
were matched empirically. Figure 3 shows the 
excellent agreement obtained over a range of 
several volts and presents some justification for 
using the relation m=2. It is to be noted that 
n=1 does not lend itself to such a check. 

It thus appears that by plotting ion yield 
curves on a semilog scale and by drawing 
tangents of a given empirical slope to the un- 
known and reference curves, potential differences 
can be determined accurately from the hori- 
zontal separation of tangent points. Since this 
method depends on the location of curves drawn 
through at least 20 points, its reproducibility 
should be considerably better than that attain- 
able by the initial onset method which depends 
on a few points close to the x axis. 


Ill. EXPERIMENTAL PROCEDURE 


All measurements were made on a consolidated 
mass spectrometer? which employs a 180° mo- 
mentum filter. Electrons are emitted by the hot 
tungsten filament F (see Fig. 4), pass through a 
hole in shield S;, enter the ionizing region through 
slit S2, and are collected at the back wall of the 
ionization chamber B and on “catcher” C. 
Secondary electrons can be eliminated by con- 
necting battery E into the catcher circuit, but 
this method introduces undesirable irregularities 
in the potential curves which presumably are 


*H. W. Washburn, H. F, Wiley, S. M. Rock, and C. E. 
Berry, Ind. Eng. Chem. Anal. Ed. 17, 74 (1945). 
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IONIZATION POTENTIALS OF. HYDROCARBONS 


due to ions formed near the catcher. Since this 
investigation was primarily concerned with the 
onset of ionization where the effects of secondaries 
are negligible, chamber and catcher were run at 
the same potential. A small cross-field (about 10 
volts/cm) between “‘pushers”’ P;, P2 and the first 
ion slit plate draws some of the ions formed 
through slit S; into the strong field (about 1000 
volts/cm) between S; and S,. Those ions that 
pass S4 are analyzed magnetically and recorded 
in the usual manner. 

Electron energies are measured on a Leeds & 
Northrup No. 7655 potentiometer which had 
been high voltage insulated to permit safe opera- 
tion as all the electron gun circuits are at high 
potential. In conjunction with a calibrated poten- 
tial divider of proper size, the potentiometer will 
measure absolute voltages to better than +0.1 
percent. The measured ‘ionizing voltage” V2- 
is attained by the electrons when they pass 
through S_ and does not include the energy 
contribution from the cross-field. This contribu- 
tion is variable since it depends on where in the 
ionizing region the ion is formed, and will con- 
tribute to the electron energy spread. Prelimi- 
nary experiments were carried out to determine 
the maximum cross-field contribution. In Fig. 5, 
the uncorrected argon potential is plotted as a 
function of V,+, the pusher potential. The slope 
of the straight line obtained measures directly 
the maximum contribution as 0.49 V,*. Its inter- 
cept on the y axis yields an argon potential 
(15.00 volts) that differs from the spectroscopic 
value (15.77 volts) by a thermal energy contri- 
bution and a certain contact potential difference. 

The energy distribution in the electron beam 
is principally a function of filament temperature. 
Pyrometric measurements reported by H. W. 
Washburn” point to a filament temperature 
range of 1650—-2000°K ; these apparent values are 
not corrected for the unknown emissivity of the 
WC layer which is believed to cover the treated 
filament. Two other factors increasing the energy 
spread are the potential drop along the filament 
(estimated as about 0.2 volt”) and the cross-field, 
as mentioned above. 

The total energy distribution was determined 
empirically by applying a fixed retarding poten- 


'’°H. W. Washburn, two private communications. 
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Fic. 7. Variation of temperature of electron energy 
distribution with magnet current Im. 


tial to the catcher (see Fig. 4) and measuring 
the catcher current J, as a function of ionizing 
voltage V2. A typical semilog plot (Fig. 6) 
yields, by the slope of its straight-line portion, 
a fictitious distribution temperature of about 
2320°K. A number of experiments were made to 
determine which of the many instrumental 
parameters would affect the energy spread. As 
expected, a decrease in filament current, and 
therefore filament temperature, reduces the 
spread. When gas is introduced, the change in 
work function of the filament surface shifts the 
entire curve shown in Fig. 6 horizontally. Since 
this gas effect is not reproducible, it is seen that 
in potential measurements the calibration gas 
must be introduced simultaneously with the 
unknown. Figure 7 shows that the fictitious distri- 
bution temperature increases as the magnet 
current is reduced below 0.3 amp. This results 
from the fact that at low field values the electron 
beam collimation is poor. Thus electrons that 
were emitted from high or low potential sections 
of the filament are able to pass the slit system. 
For this reason, special precautions must be 
taken when potentials of light compounds, such 
as Hz, De, and He, are to be determined. 

All ionization potential measurements were 
carried out under the following, strictly con- 
trolled conditions which have been found to 
yield optimum reproducibility: first electron 
accelerating potential, V;-=40 volts; cross-field 
potential (average), Vi+=4.8 volts; second ion 
accelerating potential, V;+= 600 volts ; current to 
ionization chamber and catcher, Ipi¢=47.0 A; 
magnet current J», adjusted to focus ion in 
question; reservoir gas pressure, up to 250 
microns total for calibrating gas and unknown. 
Measurements were taken by reducing the ion- 
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Fic. 8. Ionization potentials of some members of the 
normal paraffin and 1- and 2-olefin series. 


izing voltage V2~ in the critical range in steps 
of 0.1 volt and scanning each peak electrically 
several times in succession. The slight change in 
V.+ necessary to produce this sweep (about 10 
volts) did not materially alter conditions in the 
source. To test the effect of pressure reduction in 
the reservoir on the curve shape, a number of 
experiments were made by first reducing and then 
increasing V2- and averaging the values. No 
differences between this and the simpler routine 
could be detected for all compounds tested. 

The calibration gases found to be best suited 
for simultaneous introduction with hydrocarbons 
in the C,; to Cs range were argon and krypton, 
while xenon and benzene were used with all 
other unknowns. These four reference gases were 
carefully intercalibrated, and all values quoted 
below are based on the spectroscopic krypton 
potential of 14.01 volts! as a single primary 
standard. 


IV. RESULTS AND DISCUSSION 


The results given in Table I show that ioniza- 
tion potentials can be determined reproducibly 
by the slope method. On a typical semilog plot 
(see Fig. 2), a smooth curve can be drawn which 

11 Recomputed, using the new conversion factor 8067.49 
cm=1 v, from R. F. Bacher and S. Goutsmit, Atomic 


Energy States (McGraw-Hill Book Company, Inc., New 
York, 1932). 
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will fall within +0.01 volt of all experimental 
points, and the point where the critical slope is 
tangent to the curve can be reproduced graphi- 
cally to +0.01 volt. If the reference gas is 
introduced simultaneously with the unknown, 
repeat runs will yield absolute potentials that 
generally agree within 0.03 volt. Where different 
reference gases were used for the same compound, 
the maximum spread is occasionally 0.04 volt, 
while m-decane with its unusually low sensitivity 
is the only example of a 0.05-volt spread. Thus 
the experimental reproducibility of a given poten- 
tial may be expressed as +0.02 volt. 

To test the reliability of absolute values ob- 
tained by the slope method, the potentials of 
the monatomic gases A and Xe were determined 
with Kr serving as reference. It is seen that for 
these gases the impact values agree with the 
spectroscopic values within the reproducibility 
of the measurements. 

For some of the polyatomic molecules shown 
in Table I, spectroscopic ionization potentials 
computed from Rydberg series are available for 


.comparison with the electron impact values 


which are based on krypton as a primary stand- 
ard. In every case, the impact values, which 
represent ‘‘vertical” transitions favored by the 
Franck-Condon principle, are higher than the 
spectroscopic potentials (‘‘adiabatic”’ transitions 
between vibrational levels V’=0Oand V’” =0)! by 
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Fic. 9. Plot of ionization potential vs. ethylene substi- 
tutions. The last point on the methyl curve was recomput 
from Price’s value (see reference_3). 
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{ONIZATION POTENTIALS OF HYDROCARBONS 


0.02 to about 0.2 volt. This difference seems to 
be a function of the molecular complexity as it 
is smallest for diatomic molecules and largest 
for benzene. According to the Franck-Condon 
principle, the difference between the vertical and 
adiabatic transitions depends on how the minima 
of the potential surfaces representing the ground 
and ionized states are located relative to each 
other. Since the geometry of these potential 
surfaces is known exactly only in the simple 
case Of He, it is not possible to predict theoreti- 
cally such differences. However, the problem 
might be avoided by comparing the molecule in 
question to a known molecule with a similar 
ionization probability rather than to a known 
monatomic gas. 

The potentials of some members of the normal 
paraffin and olefin series are plotted as a function 
of carbon number in Fig. 8. Even though all 
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absolute values are probably high by about 0.1 
volt, the smooth curves clearly indicate that all 
potentials within any one series are self-con- 
sistent. The normal paraffin curve shows a 
systematic reduction of potential with chain 
length which is still noticeable in the Cio range. 
The potentials published in earlier electron im- 
pact studies'® for CH., CsHs, and C3Hs agree, 
within their experimental uncertainties, with the 
present values. On the other hand, Stevenson 
and Hipple’s value! for »—C,Hw (10.4 volts, 
based on J(A+) =15.77 volts) lies approximately 
0.4 volt below the value determined in this 
research and cannot be reconciled wit!. the 
paraffin curve. When extrapolated to zero carbon 
number, the paraffin curve yields the H2 poten- 
tial. Since the H—H, C—H, and C—C bonds 
are all believed to be similar’ in that a o electron 
is removed for the ionization of Hz, CH,, and 


TABLE I. Ionization potentials of standards, paraffins, and mono-olefins. The new conversion factor 
8067.49 cm=!=1 volt was used for recomputing the literature values. 








Ionization potential in volts 


Compound Individual values 


€ Literature 
<—Spectroscopic——> ‘Mass Spec.” 
Value Diff. Ref. Value Ref. 





Remarks 





Krypton 


Argon 

Xenon 

Nitrogen 
Carbon dioxide 
Acetylene 
Benzene 


12.16 
15.60 


9.45 9.43 


Methane 
Tetradeuteromethane 
Ethane 
Propane 
n-Butane 
n-Pentane 
n-Hexane 
n-Heptane 
n-Octane 
n-Nonane 
n-Decane 


Ethylene 
Propene 
1-Butene 
1-Pentene 
1-Hexene 
1-Heptene 
1-Octene 
1-Decene 


inininine sé r 
NP ODOD 


cis-2-Butene { Penn 


4ans-2-Butene 
¢is-2-Butene 
frans-2-Butene 
tvans-2-Hexene 
2-Octene 


State 
NBS 


mere RUN And~é 
CUDeKAre OF 

Re NMowNhN BSEaSa 

= 10 = 00 00 Raw wn 


2-M-Propene 
3-M-2-Butene 
2-E-1-Butene 
trans-3-Hexene 


POS COSSSe Yooos~e 


CONG OOOO BoOooy~, 
Moe Bebwhi 
Ne Un eADeK 0 


= bo Go be 


Cewo Peeves yo 
wOnru 


™ Ko 00 be 
-OO, 


14.01 

15.77 

12.14 F 

15.58 J 15.7 (14) 


13.79 
11.23 (12) 


Primary standard 
Secondary standard 
Secondary standard 
11.40 J 
9.24 ; 9.86 (15) 
13.18 (6) 
11.68 1) 
11.22 1) 
10.42 1) 


Secondary standard 


Spectroscopic value uncertain 


+0.02(?) (3) Spectroscopic value uncertain 








2 P, Kusch, A. Hustrulid, and J. T. Tate, Phys. Rev. 52, 843 (1937). 

*% G, Herzberg, Molecular Spectra and Molecular Structure. 1. Diatomic Molecules (Prentice-Hall, Inc., New York, 1935). 
4H. D. Hagstrum and J. T. Tate, Phys. Rev. 59, 354 (1941). 

‘® A. Hustrulid, P. Kusch, and J. T. Tate, Phys. Rev. 54, 1037 (1938). 
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higher paraffins, respectively, it is not too sur- 
prising that Hz fits so well into the paraffin 
' scheme. The substantial potential difference 
(about 0.17 volt) between CD, and CH, will 
be treated in detail in another paper. 

The results shown in Table I for the various 
mono-olefins may be represented graphically in 
several ways. Figure 8 shows that the 1- and 
2-olefins arrange themselves into two smooth, 
parallel curves that lie about 0.43 volt apart; 
the potentials are seen to decrease in regular 
fashion with increasing carbon number. It is 
generally postulated!*” that an olefin is ionized 
by the removal of’a"x-electron from the C=C 
bond. Thus it seems logical to treat olefins as 
substituted forms of the basic ethylene structure, 
and to study how the number, size, and position 
of substitutions affect the ionization potential. 

The top curve of Fig. 9 shows the marked and 
regular decrease of potential with the number of 
methyl] substitutions in the series ethylene, pro- 
pene, trans-2-butene, 3-methy]-2-butene, and 2,3- 
dimethyl-2-butene, which confirms Mulliken’s 
tabulation‘ of W. C. Price’s data. Since the last 
named compound was not available at the time 
of this research, Price’s recomputed value was 
used to complete the curve. It is seen that the 
ethylene potential fits such a sequence far better 
than the 1-olefin series of Fig. 8. For the three 
possible cases of double methyl substitution, 
the potential decreases slightly with increasing 
symmetry in the sequence 2-methyl-propene 
(9.35), cis-2-butene (9.31), trams-2-butene (9.29). 
In a recent paper, Dibeler? quoted the cis- and 
trans-potentials as 9.4:+0.1 and 9.1;+0.1 volts, 
respectively. In view of the symmetry sequence 
given above, Dibeler’s cis-trans difference appears 
too large to be real, but can be brought in line 
if his experimental error is larger than quoted. 
In the present research, two cis-trans pairs 


1’ E. J. Bowen, The Chemical Aspects of Light (Oxford 
University Press, New York, 1946), Chapter IIT. 
17 A, D. Walsh, Trans. Faraday Soc. 42, 779 (1946). 
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(Source: Penn State and NBS) were run a 
sufficient number of times to establish that the 
cis-trans difference measured (0.02+0.01) is real. 

The curves obtained for the ethyl and propy| 
substitution series lie below the methyl curve. 
Unavailability of pure compounds prevented 
their completion. Two cases of double ethyl 
substitution show again that the most symmetric 
structure has the lowest potential (2-ethyl-1- 
butene: 9.21; trans-3-hexene: 9.12). 4-octene was 
unobtainable to serve as an example of double 
substitution in the propyl series. However, an 
admittedly rough estimate of this point may be 
had by using the 2-octene value and applying a 
downward correction that is somewhat larger 
than the difference between 2- and 3-hexene. 
For single substitutions, the systematic lowering 
of potential with radical size has already been 
pointed out by the 1-olefin curve of Fig. 8. 
Similar tendencies for the double substitution 
sequence methyl-ethyl-propyl are apparent from 
Fig. 9. 

Earlier electron impact studies agree, within 
their experimental uncertainty, with the present 
values for propene and 1i-butene,! but their 
ethylene” and 2-methyl-propene! potentials are, 
respectively, too high and too low to fit into the 
general olefin scheme. 

The experimental data presented graphically 
in Figs. 8 and 9 show that all paraffin and olefin 
potentials are self-consistent, although they may 
be as much as 0.1 volt higher than the spectro- 
scopic values. 
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The slow rate of change of surface tension observed by Addison with aqueous solutions of 
fatty alcohols cannot be explained by diffusion alone. It is assumed that the adsorption of non- 
electrolyte solutes in the surface film is opposed by a constant potential energy barrier, and, on 
this basis, a relation between spreading pressure and surface age is developed. With the re- 
sulting equation, the experimental data of Addison are employed to evaluate the energy term. 
Nearly constant values for the assumed barrier are obtained, the actual calculated values, in 
calories per mole at 20°C, being about 1900 for isoamy] alcohol, 1300 for n-hexy] alcohol, and 700 
for n-hepty] alcohol. A second relation between spreading pressure and surface age is derived on 
the assumption that the energy barrier is a linear function of spreading pressure. The use of 
Addison’s data for isoamyl alcohol in this relation indicates that the spreading pressure- 
dependent portion of the barrier is small and of the same order of magnitude as the energy 
required to compress the surface film sufficiently to admit the entering solute molecule. 





HE curious slow attainment of surface ten- 
sion equilibrium of solutions has often been 
observed.!? In solutions of organic electrolytes 
this variation of surface tension with time is 
much slower than with solutions of non-electro- 
lytes. Addison,*-> employing the vibrating jet 
technique, has measured the surface tension of 
fatty alcohol solutions as a function of time, and 
has found that the final steady value is reached 
in times of less than one second. 

Repeated attempts? ** to describe this phe- 
nomenon in terms of rate of ‘diffusion of solute 
molecules to the surface have been unsuccessful. 
Ross? was forced to the untenable conclusion 
that a polymolecular layer was formed on the 
surface of alcohol solutions, while Addison con- 
cluded that the ‘‘migration velocity” of the solute 
molecules increased with increasing molecular 
weight. It appears that even with the non- 
electrolytes the rate of change of surface tension 
is much less than can be accounted for by 
diffusion alone. Doss? made the reasonable sug- 
gestion that an energy barrier may account for 


1N. K. Adam and H. L. Shute, Trans. Faraday Soc. 31, 
205 (1935); ibid. 34, 758 (1938). 

*K. S. G. Doss, Kolloid Zeits. 84, 138 (1938); ibid. 86, 
205 (1939). 

*C. C. Addison, J. Chem. Soc., 535 (1943). 

*C. C. Addison, J. Chem. Soc., 252 (1944). 

®°C. C. Addison, J. Chem. Soc., 98 (1945). 
(1939) N. Bond and H. O. Puls, Phil. Mag. [7] 24, 864 

7). 

7S. Ross, J. Am. Chem. Soc. 67, 990 (1945). 

(1946) F. H. Ward and L. Tordai, J. Chem. Phys. 14, 453 


the phenomenon of slow adsorption in high 
molecular weight electrolyte solutions and at- 
tempted, in the case of these solutions, to de- 
scribe the barrier qualitatively in terms of 
electrical effects at the surface. 

With solutions of non-electrolytes, electrical 
double layer effects would presumably be absent. 
The results of Addison, particularly with iso- 
amyl alcohol, thus present an opportunity to 
explore the nature of the adsorption barrier in 
a relatively simple system. 

In the somewhat similar problem relating to 
the evaporation of water through compressed 
monolayers, Langmuir and Schaefer® were able 
to describe the very low evaporation rates ob- 
served by the concept of an energy barrier of 
the form: 


\=(A+BF)RT, (1) 


where A and B are constants, and F is the sur- 
face film pressure. By analogy, it might be ex- 
pected that a similar form of barrier would 
appear in the formation of a surface film from 


.solutions. In this instance, however, since the 


solute molecules need not go entirely through 
the surface film, it would seem reasonable to 
expect the coefficient B in Eq. (1) to be smaller 
than in the process studied by Langmuir and 
Schaefer. If B is small, it would be difficult, in 
the absence of extremely accurate experimental 


I. Langmuir and V. J. Schaefer, J. Frank. Inst. 235, 
No, 2, 119 (1943). . 
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TABLE I. Values of constants of Eq. (7) for 
various alcohols at 20°C. 











J D 
Alcohol (cm? X 1016) L (cm? sec.~! X108) 
Isoamyl 16.7 0.827 7.9 
n-Hexyl 19.7 0.710 7.1 
n-Heptyl 18.4 0.600 6.4 








data, to differentiate between a constant energy 
barrier and one varying with film pressure. 
Depending upon the equation of state de- 
scribing the surface film and the form of the 
assumed barrier function, a number of different 
relations between spreading force and time are 
obtained. We shall herein derive two such rela- 
tions, employing an equation of state applicable 
to films of soluble alcohols and two simple poten- 
tial functions of the Langmuir type. The extent 
to which either of these relations can be fitted 
to the experimental results will serve to indicate 
the possible correctness of assumptions involved. 


CASE I 


A constant energy barrier and a surface film 
obeying the two-dimensional Amagat equation 
are assumed. 

Solutions of many soluble substances, includ- 
ing the fatty acids and alcohols, form surface 
films obeying the equation!® ) 


F=N(FJ+L4T), (2) 


TABLE II. Energy barrier to adsorption of alcohols in 
surface film (at 20°C). 











Conc. of 
solution Fy 4 r 
(Wt.%)  (dynescm=!) (sec.X108) exp(A/kT) (cal. mole!) 
Isoamyl alcohol 
0.122 4.3 12.4 27.6 1930 
0.203 6.1 8.0 28.2 1940 
0.325 9.3 5.7 28.4 1950 
0.406 10.4 5.3 31.6 2010 
n-Hexyl alcohol 
0.050 5.4 24 10.4 1360 
0.100 9.0 9.5 8.9 1270 
0.150 11.4 6.0 9.3 1300 
0.200 13.2 4.9 10.3 1350 
n-Heptyl alcohol 
0.020 6.1 36 3.30 690 
0.040 10.0 17 3.31 700 
0.060 12.6 10 3.38 710 








10N. K. Adam, The Physics and Chemistry of Surfaces 
(Oxford University Press, London, 1941), p. 118. 


where N is the number of molecules of solute 
per cm? of surface and J and L are constants. 

Ward and Tordai® have shown that diffusion 
theory leads to the following expression for the 
rate of arrival of diffusing solute molecules at 
one cm? of surface of a solution, 


dM/dt=Dno/(xDt)', (3) 


where mp» is the bulk concentration of solute and 
D is the diffusion constant. 

We have assumed, however, that the actual 
adsorption of molecules in the surface film is 
opposed by a constant barrier of height \. 
Therefore, the rate of adsorption of solute is 
given by 

dN/dt=0M/dt exp(—A/RT). (4) 
Hence, 


N =2no(Dt/x)* exp(—A/RT). (5) 


Combining Eq. (5) and Eg. (2), we obtain the 
desired relation between spreading pressure and 
time: 


F/FJ+LkRT =2no(Dt/r)' exp(—A/kT), (6) 
which may be written in the form, 
exp(A/RT) = 20 FJ+LkT)(Dt/x)3/F. (7) 


Equation (3) was originally obtained by neglect- 
ing any backward diffusion of solute molecules 
from the surface, and accordingly can only be cor- 
rect at times short of equilibrium. In order to test 
Eq. (7), therefore, we shall consider the times 
required for the surface pressure to go only half 
of its way to equilibrium. The work of Addison*® 
on the rate of change of surface tension of solu- 
tions of isoamyl alcohol, n-hexyl alcohol, and 
n-heptyl alcohol furnishes enough experimental 
data to permit calculation of the approximate 
values of the energy barriers by Eq. (7). 

The required values of J and L may be ob- 
tained from the data given by Addison on the 
equilibrium surface tension and the rate of 
change of surface tension with concentration for 
solutions of the above mentioned alcohols. The 
diffusion constants (at 20°C) were taken from 
the curve of Ward and Tordai.* Table I gives 
the numerical values of the various constants of 
Eq. (7) so obtained, at 20°C. 

Table II gives the spreading force-time data 
for the three alcohols as interpolated from 
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TIME-DEPENDENCE OF SURFACE TENSION 


Addison’s experimental results, and the values 
of the energy barriers as calculated from these 
data and those of Table I by insertion in Eq. (7). 
In Table II, Fi denotes the spreading force of 
the surface film when the surface tension has 
proceeded one-half of its way toward equilibrium, 
and # is the time required for the spreading force 
to reach this value. Fy and & are used for F and ¢ 
in Eq. (7) for the reason previously cited. 

It is evident that with each alcohol the calcu- 
lated energy barrier values are essentially con- 
stant within the range of accuracy of the spread- 
ing force measurements. However, there is a 
consistent suggestion in the trend of calculated 
values that the barrier may increase slightly with 
increasing film pressure, or, in other words, that 
B of Eq. (1) should have a small positive value 
rather than the value of zero assumed in this 
case. It is also found that the height of the 
barrier decreases markedly with increase in 
molecular weight of the alcohol, a result which 
may be related to the fact that the lateral ad- 
hesion of the molecules in a surface film increases 
rapidly with increasing hydrocarbon chain length. 


CASE II 


An energy barrier of the Langmuir type, 
and a surface film obeying the two-dimensional 
Amagat equation are assumed. 

Expressed formally, these assumptions are that 


\=(A+BF)kT* 
and 
F=N(FJ+LkT).** 


Therefore, 


dN /dt=LkT/(FI+LkT)?(dF/dt). 


Substituting these values and Eq. (3) in Eq. (4), 
we obtain 


exp(BF)dF/(FJ+LkT)? 
=n(D/zt)* exp(—A)/LkTdt. (9) 


Integration of Eq. (9) leads to the desired 
telation between spreading pressure and time, 


* See Eq. (1). 
** See Eq. (2). 


namely 


exp(A) —LkT exp(A+BF)/FJ+LkT 
+BLkRT/J{ Ei. B(FJ+LkT)/J] 
— Ei. BLkT/J]} exp(A —BLkT/J) 


=2nJ(Dt/x)* (10) 


where Ei is the exponential integral function." 

When the coefficient B has the value of zero, 
the conditions are those of Case I, and Eq. (10) 
reduces to Eq. (7). 

For accurate evaluation of the constants A 
and B, the use of Eq. (10) is somewhat cumber- 
some, and, in addition, requires more precise 
experimental data on the relation between spread- 
ing pressure and time than are now available for 
non-electrolyte solutions. However, in the case 
of isoamyl alcohol, the appropriate constants 
from Table I and the experimental data from 
Table II have been employed in Eq. (10) to 
obtain approximate values of A and B. These 
are found to be 3.2 and 0.040, respectively. 
Therefore, for the adsorption of isoamyl alcohol 
in the surface of its aqueous solution at 20°C, 
the apparent energy barrier is given by 


d= (3.2+0.040F)kT 
or, in calories per mole at 20°C, 


Acal./mole = 1860+-23F. (11) 


This result conforms with the indicated slow 
increase of the energy barrier value with in- 
creasing spreading pressure observed in Case I. 


DISCUSSION OF RESULTS 


The results of calculations presented in Table 
II and under Case II show that the concept of 
a potential energy barrier, of reasonable magni- 
tude, hindering the adsorption of solute molecules 
in the surface is sufficient to describe the slow 
attainment of surface tension equilibrium in 
solutions of some organic non-electrolytes. For 
more accurate determinations of the energy 
barriers and their dependence on spreading pres- 
sure, more precise experimental data on the 
relation between surface tension and time would 


11 Numerical values of the exponential integral function 
are given by E. Jahnke and F. Emde, Tables of Functions 
(Dover Publications, New York, 1943), p. 6 et. seg.; and 
by A. N. Lowan, Tables of Sine, Cosine, and Exponential 


Integrals (Works Progress Administration, New York, 


1940), Vols. I and IT. 
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be desirable. Also, experiments on this relation 
at different temperatures would be helpful in 
evaluating the barrier. 

In considering the possible origin of a potential 
barrier to adsorption of a solute molecule, there 
immediately appears the work required to com- 
press the surface film sufficiently to make space 
for the molecule entering the film from solution. 
This energy requirement would be roughly equal 
to Fa, where a is the effective area occupied by 
the entering solute molecule. Such a potential 
barrier would, of course, be directly propor- 
tional to the spreading pressure and would 
increase with increasing molecular weight of 
solute. We have found, for the fatty alcohols 
discussed above, only a very slight dependence 
of \ on F and a smaller value of \ for the higher 
alcohols. The value of Fa will also be far too 
small, in general, to account for the A values 
calculated. Taking isoamyl alcohol as an ex- 





L. E. JACOBS, J. R. PLATT, AND G. W. SCHAEFFER 


ample, and using a molecular cross-sectional area 
calculated from its density at 20°C, we find Fa 
to have a value of about 57F cal. per mole. 
This is of the same order of magnitude as the 
spreading pressure-dependent portion of the po- 
tential barrier calculated for isoamyl alcohol and 
given in Eq. (11). It thus appears that, at most, 
the work of compressing the surface film can 
only account for the small portion of \ which is 
dependent on spreading pressure. 

Although in carrying out the present calcu- 
lations, Fy and # values have been employed in 
Eq. (7) and Eg. (10), it should be pointed out 
that the use of other corresponding F and } 
values up to the point where F is at least 80 
percent of its equilibrium value yields essentially 
the same values of \. Therefore, it appears that 
back diffusion from the surface probably leads 
to no important errors in the calculations, within 
the range of F values used. 
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Search for Near Ultraviolet Bands in Borazole* 


Lots E. Jacoss, J. R. PLATT 
Physics Department, University of Chicago, Chicago, Illinois 
AND G. W. SCHAEFFER 
Chemistry Department, University of Chicago, Chicago, Illinois 
(Received October 1, 1947) 


Borazole, BsNs3He, shows no absorption bands in the region 2000A-—4600A, at pressures of 
about 2 atmospheres in a 4-cm cell. Since bands analogous to those of benzene are expected on 
theoretical grounds; it is concluded that the analogues of the benzene 2600A and 1830A bands 
are the borazole 1995A and 1720A bands, with the analogues of the benzene 2080A bands 





perhaps being hidden in between. 


N attempt has been made to find bands in 
triborine-triamine or borazole, B3N3He, 
corresponding to the 2600A bands in benzene— 
that is, having the same general appearance as 
the latter bands and lying, like these, in the 
quartz ultraviolet. Pressures of about 2 atmos- 
pheres, higher than any previously used in this 
experiment, were obtained but with a 4-cm path 
no bands were found in the region 2000A—4600A. 
This experiment was suggested by Dr. E. H. 
Hutten. He attempted earlier to detect these 


* This work was assisted by the Office of Naval Research 
under Task Order IX of.Contract N6ori-20 with the Uni- 
versity of Chicago. 


bands in this laboratory, but only enough bora- 

*zole was available at that time to obtain vapor 
pressures of about 20 cm in a 4-cm cell. His 
work was the result of the belief that the ring 
structure and electronic structure of borazole 
are very similar to those of benzene (see Bibli- 
ography in reference 1). 

In support of this belief was later work from 
this laboratory on borazole solutions! in which 
bands at 1995A and a peak at 1720A seemed to 
be similar to the benzene bands at 2080A and 


1J. R. Platt, H. B. Klevens, and G. W. Schaeffer, J. 
Chem. Phys. 15, 598 (1947). 
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the strong N->V transition peak at 1830A, 
respectively. Also, Dr. R. S. Mulliken has shown 
that on theoretical grounds “‘a one-to-one corre- 
spondence and close resemblance between the 
electron configuration and over-all states of 
borazole and benzene are to be expected”’ (foot- 
note 12 of reference 1). In these solution experi- 
ments the bands similar to the 2600A benzene 
bands were not found even when concentrations 
20 times that necessary for detecting these 
bands in benzene were used. 

For the present work, 2.4;X10-* moles of 
borazole were heated to 77°C, or until all the 
material appeared to be in the vapor state, in a 
quartz cell of 4-cm length and approximately 
35.4-cc volume, giving a concentration of 6.87 
X10-* m/l. The borazole was prepared by the 
method indicated in reference 1. Plates taken on 
the small University of Michigan vacuum fluo- 
rite spectrograph! showed no bands in the 2000— 
4600A region. The benzene 2600A bands ob- 
served in a matching cell were still distinctly 
visible when the side arm was at a temperature 
of —68°C, corresponding to a pressure of about 
0.042 mm. It was estimated that they would be 
barely visible at pressures about five times lower, 
corresponding to a concentration of 6.8107’ 
m/l. This would indicate that if bands like the 
2600A benzene bands occur in borazole in the 
2000-4600A region, their intensity must be less 
than that of benzene by the factor of at least 10°. 
Therefore, one might conclude that no bands 
exist in borazole corresponding to the 2600A 
benzene bands, or that if they do exist, they 
have shifted to shorter wave-lengths and the 
upper state has become so mixed with other 
states that the bands are unidentifiable as such. 

Another explanation which seems more plaus- 
ible in view of the theoretical reasons for ex- 
pecting such bands, is that the 1995A borazole 
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bands should not be identified, as was done in a 
former paper,’ with the 2080A benzene bands, 
but are instead themselves the analogs of the 
2600A bands. On this view, it would be the 
analogues of the 2080A benzene bands which are 
concealed, perhaps being hidden by the strong 
N-V peak at 1720A. This interpretation would 
reduce the spread of the states of the upper 
electronic configuration from its value of about 
15,000 cm between the 'Bz, and '£j, states 
benzene to about 5000 cm~ between the '/ » 
and !E’ states of borazole, but it would preserve 
the sequence of corresponding states in the 
compounds. Also the comparative weakness of 
the 1995A bands, which was a difficulty in the 
former interpretation making them analogous 
to the rather strong 2080A benzene bands, now 
becomes quite reasonable, as they are several 
times as strong as the 2600A benzene bands. 
Their strength and their diffuseness may then 
be the result of their closeness to the strong 
N-V transition. 

The whiteness of boron nitride crystals, which 
was commented on in the former paper, also 
becomes somewhat more understandable on this 
interpretation. The resonance interaction which 
produces the spread in energies among the states 
of the upper electron configuration now appears 
to be only one-third as strong in a B;N; ring as 
it is in a Cs ring. Consequently it is no surprise 
that in a lattice of BsN; hexagons the resonance 
interactions do not produce even the 20,000- or 
30,000-cm= spread necessary to give a visible 
color, even though in the corresponding carbon 
lattice of graphite the states are spread through- 
out 60,000 cm or more and the compound is 
black. 

We are indebted to Mr. Riley Schaeffer for 
the preparation of the borazole used in the 
experiment. 


/ 
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Molecular Orbital Treatment of the Ultraviolet Spectra of Benzene and Borazole} 


C. C. J. RoorHAan* AnD R. S. MULLIKEN 
The University of Chicago, Chicago, Illinois 
(Received October 1, 1947) 


The ground state and the first excited states of benzene and borazole are treated according 
to the semi-empirical molecular orbital method. Comparing the observed spectra with the 
computed energies, a spectroscopic value of —2.6 ev for the bond integral 6 is obtained in 
the case of benzene. In the case of borazole this procedure includes one more parameter, which 
takes care of the difference in electronegativity of boron and nitrogen. Taking a reasonable 
value for this parameter, a value of —2.3 ev for the bond integral between boron and nitrogen, 
8*, is obtained in fair agreement with theoretical expectations. 





HE benzene molecule has been treated 
according to the molecular orbital method 
by Mayer and Sklar.'! These authors made a com- 
pletely theoretical computation. For the Hamil- 
tonian, they took for the field of the core the 
field of six carbon atoms minus six 27 electrons, 
then fed these six electrons into this field and 
calculated the energies of the ground state and 
the first excited states. The only empirical value 
they used was 1.40A for the C—C distance. 

In the present paper a semi-empirical method 
is used. The procedure is equally simple for 
benzene (CsH¢) and borazole (B;N3Hg), the only 
difference being the occurrence of one more 
parameter in the case of borazole. 

Neglecting interactions of non-neighbor ring 
atoms, we can write the secular equations for 
the six approximate molecular orbitals, derivable 
as linear combinations of 2pz atomic orbitals of 
the ring atoms, in the form: 








A i000 1 
Paot:sd & @ 
Pon 2.6.2 
“a ar ae a =0 (benzene), 
600141 
, 2 @& © 1 -s 
- 0 O 0 1 
i ._ 2 aa ae 
= ae ae 0 O 
0 0 a 0 =0 (borazole), 
oh fee ee Gee ee 
1 OO oO O ; 2 








+ This work was assisted by the Office of Naval Research 
under Task Order IX of Contract N6ori-20 with. the Uni- 
versity of Chicago. 





where 
a—E : a’—E at a’—E 
iit yt — S*E yt — S*E 


7 f Selleddr, «x f Jo'He'dr, 
al" = f ut bu'dr 


in J Jcllvetdr, y*= f Det pnd, 


S= | Vo'verdr 5* = J ¥ordn'dr. 

The y’s are atomic 27 orbitals, the subscripts 
referring to the kind of atom, the superscripts to 
the location in the ring. H and H* are the one- 
electron Hartree-field Hamiltonian operators for 
benzene and borazole, respectively. Substituting 


a—-E=-X, 
y—Sa=8, 
y*—S*a=6*, 
a’—a=9', 


o!'-a-s s”’, 


we get 
—X 5’—-X 6” —X 
a, ‘=—____., Al =——. 
B—SX B*—S*X B*¥—S*X 


* At present at the Catholic University of America, 
Washington, D. C. - 
1M. G. Mayer and A. L. Sklar, J. Chem. Phys. 6, 649 
(1938). See also V. Griffing, J. Chem. Phys. 15, 421 (1947). 
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SPECTRA OF BENZENE AND BORAZOLE 


The a’s, y’s, and 6’s are known to be negative 
quantities. The a’s are electronegativity param- 
eters, and can be taken proportional (with 
reversed sign) to the absolute electronegativities.’ 
In Pauling’s electronegativity scale,’ the electro- 
negativity increases by approximately a constant 
amount for successive members in the second 
row of the periodic table. We therefore assume 
—§' =" =5(<0). 

By a suitable similarity transformation‘ the 
secular equations can be reduced to the forms: 


coohmrwl lo 








oo 


Pe ee 








0 
0 
0 
, 0 
A 
1 


oo 


A” 


(The equation for benzene could be reduced 
completely by a different transformation based 
on the higher symmetry of benzene; however, 
in the given form the similarity between benzene 
and borazole comes out more clearly.) The 
molecular orbitals corresponding to the roots 
A=-—2 and A=~-—1 (twice), in benzene, or the 
corresponding orbitals in borazole, are each 
occupied in the ground state by two electrons. 
The average height above the ground state of 
the first group of excited states'® (three singlets 
and three triplets), corresponding in benzene to 
the excitation of one electron from one of. the 
A=-—1 orbitals to one of the A=-+1 orbitals, 
is given by the difference X,—X, of the two 
roots of 


A*—1=0, ie., of X?—(p—SX)?=0; 


a S. Mulliken, J. Chem. Phys. 2, 782 (1934); 3, 573 


(1 

*L. Pauling, The Nature of the Chemical Bond, (Cornell 
University Press, Ithaca, New York, 1940), p. 65. 
_ ‘ These transformations will be treated in a future pub- 
lication. 

® See footnote 12 in J. R. Platt, H. B. Klevens and G, W. 
Schaeffer, J. Chem. Phys. 15, 598 (1947), 


namely 


for benzene; and for borazole, of 
A'A" —1=0, i.e. of X2—62—(6*—S*X)?=0; 


namely 
2* 
X_—X),= ————(1+ (1 —S**)d?)}, 
— S*2 
where d= 6/8*(>0). 
The overlap integrals S and S* can easily be 
evaluated if we use Slater 27 orbitals. The 
results are:*7 


2 1 
S= (: +x +ot+—xt es 
5 15 


x 


for like atoms, with x=Z’R/n*, where Z' =ef- 
fective nuclear charge (3.25 for carbon), n* =ef- 
fective quantum number (2 for 2 orbitals), and 
R=internuclear distance, in atomic units. 


(ab)? )> 1 1 
s=| —{ coshn+~ sinh 
3(a+b)) 1»? n 


1\1+x 1 
+3(1-—) (coshns -— sinh) fe 
9g? x* nx 


for unlike atoms, with 
a=Z'/n*, b=Z"/n*, 
x=3(a+)b)R, 
a—b 


n= n* =2. 


a+b 


Z' =2.60 for boron, 
Z” =3.90 for nitrogen, 


The actual internuclear distances in benzene and 
borazole are 1.39A and 1.44A, respectively. With 
these values we get S=0.25, S*=0.22, and for 
the mean energies from the ground state to the 
first excited states: 


E.— E\ = —2.13~, 
E,—E;= —2.108*(1+0.952d?) 3. 


benzene: 
borazole: 


As was mentioned above, corresponding to the 
mean energy Ee, a group of six electronic states 


6 Handbuch der Physik, Vol. XXIV;, p. 643. 
7A. L. Sklar, J. Chem. Phys. 7, 991 (1939). 
8S. H. Bauer, J. Am. Chem. Soc. 59, 1804 (1937). 
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benzene borazole 
: » i 
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Ra -———A0 
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Oe} 


1A y———0.0 


Calculated 
(Mayer and Sklar) 


—_—_—§0.0 —_———0.0 
Observed 


Fic. 1. Calculated and observed electronic energy levels, 
for vertical excitation (that is, all internuclear distances 
remaining fixed during the excitation process). 


(three triplet and three singlet states) is expected 
theoretically for each of the two molecules.! 

In Fig. 1 we give the results of the calculations 
of Mayer and Sklar on benzene, and the energy 
levels for benzene and borazole based on ob- 
served spectra.*" Since the present calculations 
apply to vertical excitation, that is, to electron 
jumps without change in nuclear configuration, 
the frequencies of the intensity maxima of the 
band systems in the observed spectra correspond 
to our energy levels. (For those spectra which 
were observed in liquid or solution, we have 
estimated and corrected for the resulting red 
shift.) The energies above the ground state are 
given in ev. The group theory notation is essen- 
tially that of Placzek and Mulliken, slightly 
modified by Sponer and Teller.” 

For transitions from the ground state (A, in 
benzene, 'A;’ in borazole) to the six excited 
states, group: theory gives the following selection 
rules; corresponding states are listed side by side: 


* For benzene: H. Sponer, G. Nordheim, A. L. Sklar and 
E, Teller, J. Chem. Phys. 7, 207 (1939), where numerous 
references on experimental data are given; H. Sponer, J. 
Chem. Phys. 8, 05 (1940). 

10 For borazole: J. R. Platt, H. B. Klevens and G. W. 
Schaeffer, J. Chem. Phys. 15, 598 (1947). 

"For borazole: L. Jacobs, J. R. Platt, and G. W. 
Schaeffer, J. Chem. Phys. 16, 116 (1948). 

2 R. S. Mulliken, Phys. Rev. 43, 279 (1933); H. Sponer 
_ E. Teller, Rev. Mod. Phys. 13, 75 (1941): see Table 


MULLIKEN 


Benzene Borazole 
To'E,- '£’ (strongly) allowed 
1B,, 1A,’ electronically forbidden, but 
(weakly) allowed in combina- 
tion with vibrations of species 
E,*, E’ (or Bag, A2"") 
1Bo,, electronically forbidden, but 
(weakly) allowed in combina- 
tion with vibrations of species 
ZE,*, 2 
5Bo, electronically and spin for- 
bidden 
3E.- %E’ — spin forbidden 
3B,, 3A,’ electronically and spin for- 
bidden 


The observed benzene spectrum shows a 
strong band system at 7.0 ev, two weaker ones 
at 6.2 and 4.9 ev, and a very weak one at 3.8 
ev (the energies given are those for peak absorp- 
tion in each case). In obtaining these values, 
besides the liquid or solution correction already 
mentioned, we subtracted one vibrational quan- 
tum of the E,+ or E’ vibration which made the 
transitions from the ground state to the B or A’ 
levels allowed. 

The 7.0-ev band system is the allowed transi- 
tion ;“ the 6.2- and 4.9-ev band systems are the 
transitions to the 1B,,, 4 and 1B», 154 levels, re- 
spectively. The identification of the 4.9-ev band 
system has been thoroughly established by a 
careful vibrational structure analysis.*%* The 
3.8-ev band system!!* must be interpreted as 
the transition to the *B,,, 1° (or the *B2, !7) level, 


13R. S. Mulliken, J. Chem. Phys. 7, 984 (1939). 

4 The present interpretations of the 6.2- and 4.9-ev 
systems were first suggested as a possibility by Mayer 
and Sklar (reference i The definite interpretations were 
first given by G. Nordheim, H. Sponer, and E. Teller, 
J. Chem. Phys. 8, 455 (1940), as modified by R. S. Mul- 
em ‘yg footnote 9 in the paper of Nordheim, Sponer, and 
Teller). 

1 A, L. Sklar, J. Chem. Phys. 5, 669 (1937). 

158 The vibrational analysis shows that the 4.9-ev transi- 
tion is electronically forbidden, and vibrations of sym- 
metry E,* which make the transition allowed, can be 
indicated. This proves that the electronic level is either 
1B,, or 1By,. A transition to a 1Bi, level should also be 
made allowed by a vibration of symmetry B2,; the presence 
of such a vibration. in the band would establish the sym- 
metry of the electronic level as 'B,, for certain, its absence 
would make 1B2, probable. A more detailed vibrational 
analysis of the 4.9-ev band might confirm or disprove the 
labeling 'Bo,,. 

16 $ N. Lewis and M. Kasha, J. Am. Chem. Soc. 67, 994 
1945). 

17 From the numerical values of the integrals computed 

by Mayer and Sklar:the *B1, level should be the lowest, the 
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since the low intensity in absorption and the long 
phosphorescence lifetime!® indicates a doubly 
forbidden transition, and its appearance in 
phosphorescence shows that it must be from the 
lowest triplet level.'®!® Since no other band 
systems have been found in the visible or near 
ultraviolet, and the transition to the *E,~ state 
is expected to be about 10* times stronger than 
the observed singlet-triplet transition, the *Z,~ 
level must be at least as high as 4.7 ev; the 
transition to this level could then be masked by 
the 4.9-ev band system, which is expected to be 
10° times stronger.?° Since a larger shift is im- 
probable, we shall in estimating the locations of 
the unobserved levels assume 4.7 ev for the 
‘E,~ level. For the remaining unobserved *Be, 
(or *B,.) level we shall assume 5.6 ev, since 
Mayer’s and Sklar’s calculations give equal 
separations among the three triplet levels. The 
fact thus probably established, that the *B2, (or 
3B,~) level must be above the lowest singlet level, 
means that the calculations of Mayer and Sklar 
not only give numerical values of some inter- 
action integrals in disagreement with experi- 
ment, but even with the wrong sign. However, 
as these authors pointed out, any orbital type of 
approximate calculation tends to give triplet 
states somewhat too low compared to singlet 
states, 

The observed spectrum of borazole, if we 
accept the interpretation given in reference 11, 
shows the allowed transition at 7.7 ev, and the 
equivalent of the 4.9-ev transition in benzene at 
6.5 ev. To complete the picture of the energy 
levels, we assume the splittings of the six levels 
to be proportional to the splittings in benzene. 


*B», level the highest, among the triplet states. However, 
if the true values of the integrals should be somewhat 
different, the *Bs, might be the lowest, the *B,, the highest. 
In any event, the theoretical prediction that the *£,~ 
should lie in the middle is independent of the numerical 
values of the integrals. 

84.0 sec. (private communication from M. Kasha). 
This agrees roughly as to order of magnitude with the 
lifetime of about 0.1 or 0.2 sec. calculated from the ob- 
served absorption. (The absorption f value, estimated from 
Fig. 13 on p. 1000 of reference 16, is about 1.2 X 10-8.) 
(1944), N. Lewis and M. Kasha, J. Am. Chem. Soc. 66, 2100 

*° The transition to *B,, is about 10-5 times as strong as 
that to 1Be,; since the transition to 1B, is about 10-* times 
as strong (see I. R. Platt and H. B. Klevens, Chem. Rev. 
41, 301 (1947)) as that to 1Z,-, we expect the transition 
to*E,~ to be about 10* times stronger than that to *Bi, and 
thus about 10? times weaker than the 1B, transition. 


SPECTRA OF BENZENE AND BORAZOLE 


















benzene borazole 

1F! 7.7 

14 ,/— — —7.2 

ee 7.0 34 ,/— — —6.9 

14 «,———_6.5 

1Bi,—————6.2 3k’ ———6.4 

34 ,/———5.9 
3Bo,— — —5.6 
1Bs,——_—+4.9 
sEL-— — —4,7 
3By,———3.8 

14 1,———0..0 14 ,/————0.0 


Fic. 2. Observed and (dashed lines) estimated electronic 
energy levels, for vertical excitation. In benzene it is fairly 
likely that the labels of the two *B levels should be inter- 
changed, and similarly for the two 4A levels of borazole. 
Similar reversals may also be needed for the corresponding 
singlet levels.'* 


The final results for both molecules are shown in 
Fig. 2. : 

In the formulas which we derived for the 
energy difference E2,—E, between the first 
excited states and the ground state in terms of 
the bond integral 8, E» refers to the average of 
the energies of the six levels we have just been 
discussing (in averaging, each E level must be 
taken with double weight since these levels are 
orbitally doubly degenerate; corresponding sing- 
let and triplet levels are, however, given equal 
weight, a procedure whose correctness can be 
seen by referring to ref. 1). The resulting aver- 
ages are 5.5 ev for benzene and 6.8 ev for bora- 
zole. (These averages do not depend on a change 
in labels of the B or A levels.'®* 7.) 

The formula then gives for benzene, if we 
designate by spec the spectroscopically deter- 
mined 8, Bspec = —5.5/2.13 = — 2.6 ev. From spec- 
troscopic data on ethylene, the value Bspec= 
— 3.0 ev is obtained,” again for a 2b” C—C bond 
although for a somewhat smaller internuclear 
distance. From resonance energy data another 8 
value van be obtained for benzene:™ By.= 
—3.1 ev. 

For borazole, we find 


Bepec* = —6.8/2.1(1.95)§= —2.3 ev 
if we assume d=1.0, 


2 R. S. Mulliken and C. A. Rieke, Rev. Mod. Phys. 14, 
259 (A) (1942); and full papers to be published shortly, in 
which the significance of bond integrals from spectroscopic 
and other data will be discussed. 

































































































































































































































































































or 


Bspec* = — 6.8/2.1(3.14)= —1.8 ev 
if we assume d=1.5. 


This gives for 6 the respective values —2.3 and 
—2.7 ev. The first choice gives the closest re- 
semblance between benzene and borazole and 
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gives at the same time what is probably a 
reasonable value for 6.? 

It can be seen from the formulas for E2,— 
that the main reason for the higher frequency of 
the borazole spectrum as compared with benzene 
is the appearance of the parameter d in the 
formula for borazole. 
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The Effect of Certain Column Variables on Thermal Diffusion 


H. G. Drickamer, V. J. O’Brien, J. C. BRESEE, AND C. E. OCKERT 
Noyes Laboratory of Chemistry, Urbana, Illinois 
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Tests have been made on two thermal diffusion columns with different annular spacing. The 


results indicate that the dependence of the separation on temperature and annular spacing 
predicted by Furry, Jones, and Onsager holds. The dependence on pressure is slightly different 
than was predicted. The possibility that the index of repulsive force may depend on the pressure 


is pointed out. 








NOMENCLATURE 


C,=specific heat at constant volume 
D1:=coefficient of ordinary diffusion 
Dr=coefficient of thermal diffusion 

g=gravitational constant 
L=column length 
p=pressure in atmospheres 
q=separation factor (X17X22/X27X \) 
T =absolute temperature, °K 
T =arithmetic average temperature, °K 
AT=temperature difference between hot and cold 
walls 
V1, ¥g=convection velocity of species 1, 2 
X,, X2:=mole fraction of species 1, 2 in feed 
X,7, X27 =mole fraction of species 1, 2 in top reservoir 
X,*, X28 = mole fraction of species 1, 2 in b.t.m. reservoir 
w=one half of annular spacing 
a=thermal diffusion ratio= (D7/Dj2) 
6, e=empirical constants 
n = viscosity 
\=thermal conductivity 
p=density 
v=index of repulsive force in F= —(Kj2/r’) 


INTRODUCTION 


HERMAL diffusion—the transport of mole- 
cules caused by a temperature gradient, 
was predicted theoretically by Enskog! in 1911 


1D. Enskog, Physik. Zeits. 12, 56, 533 (1911). 






and independently by Chapman? in 1916. Its ex- 
istence was confirmed experimentally by Chap- 
man and: Dootson* in 1917. With the develop- 
ment of the thermal diffusion column with hot 
wire by Clusius and Dickel‘ and the refinement to 
a concentric tube column by Brewer and 
Bramley,’ it became a practical process for 
isotope separation. 

Furry, Jones, and Onsager‘ and later Furry and 
Jones’? have presented an extensive analysis of 
the theory of the column. This paper reviews 
their results briefly and compares the predicted 
effect of certain column variables on the separa- 
tion with experiment. 

THEORY 

Furry, Jones, and Onsager® use the defining 

equation : 


X1X2(%1 —Ve) 
= —Dy. gradX,:+Dr grad In7, (1a) 


=a S. Chapman, Phil. Trans. 216A, 279 (1916). 
( 3S. Chapman and F. W. Dootson, Phil. Mag. 38, 248 
1917). 

4K. Clusius and G. Dickel, Naturwiss. 26, 546 (1938). 

935) K. Brewer and A. Bramley, Phys. Rev. 55, 590(A) 
(1939). 

6 W. H. Furry, R. C. Jones, and L. Onsager, Phys. Rev. 
55, 1083 (1939). 

7W. H. Furry and R..C. Jones, Rev. Mod. Phys. 18, 
219 (1946). 















=S = fA = PA a he 


th 









-F, 
y of 
ene 

the 


1948 


S ex- 
hap- 
elop- 
1 hot 
nt to 

and 
; for 


y and 
sis of 
views 
licted 
para- 


fining 


(1a) 


38, 248 


938). 
590(A) 


'S. Rev. 


ys. 18, 





or 


X1iXe(% —Ve) 
=D,,(—gradX1+X1X20 grad InT), (1’) 


for the thermal diffusion coefficient, D7, and 
thermal diffusion ratio, a. Applying the Navier 
Stokes Law for viscous flow to a simplified 
hydrodynamical model of the column, a general 
differential equation is obtained expressing the 
flux of one component in terms of the molecular 
properties of the gas and the constants of the 
column. 

To solve this equation it is necessary to know 
the temperature dependence of the molecular 
properties and, therefore, to select a satisfactory 
molecular model. One of the simplest models 
giving approximately satisfactory results is the 
point source of repulsive force model, 


F=(—Ki2/r’). (2) 


For this model Chapman and Cowling® give 
the following relationships for the gas coefficients : 


n=Q(T, mM, v) =O Tt3/20-1) | (3a) 


T (3+) /2(e—1) 








Dyp=Y(T, m,X, ¥, p) =y’ (3b) 
h=7(v)nCs, (3c) 
a= (vy, m, x)(v—5/v—1), (3d) 

Dr =X \X2aD. (3e) 


It should be pointed out that more accurate 
molecular models indicate that a depends on the 
temperature. 

The general equation cannot be solved, but if 
one ignores vertical diffusion the equation sim- 
plifies and can be solved using (3a)—-(3d). A 
correction for vertical diffusion can then be 
appended to the solution. The solution is par- 
ticularly simple if one considers that the tempera- 
ture dependence of Diz, 7, and X is that of 
Maxwellian molecules (v=5). For many gases 
this simplification adds little inaccuracy. 

One then obtains for the separation factor 





8S. Chapman and T. G. Cowling, The Mathematical 
Theory of Non-Uniform Gases (The Cambridge University 
Press, New York, 1939). 
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q=(X17X28/X27X1*) the relationship 
63LaDnf(AT/T) 


Ing= 





2w‘pg 





1 
x| | (4) 
1+ 5670(DinT /w®gpAT)? 


The factor in the brackets is the correction for 
vertical diffusion. 

The quantity f(A7/T) is essentially equal to 
one for (AT/T) £0.50. 

Solving Eq. (4) for the thermal diffusion 
coefficient one obtains 


2 X:1X w'pg Ing 


;=— 


63 Ly 








DynT 4 
— | (5) 


Xx [ a 5670( 
pew'AT 


It is the purpose of this paper to consider the 
effect of certain column variables on the separa- 
tion and to compare predicted values with 
experiment. 

Assuming the Maxwellian temperature de- 
pendence of Dis, and n, a independent of temper- 
ature, and the ideal gas law, one obtains for (4) 


CT 1 
= — | (4’) 
w'PL1+(CoT/ Py) 





Ing 


where C, and C2 are independent of 7, p, and W. 
Equation (4’) can be maximized with respect to 
p, T, or W at constant values of the other two 
variables. Figure 1 shows Ing plotted against 
pressure, with temperature and VW as parameters. 
For this case the following values were assumed 
for the different variables: 


L=250 cm 
Dy»=0.1X 10T* 
7=0.3X10°T 
p=0.488(7T/P) 
.a=0,.012 


Different values for these functions would merely 
shift the curves without changing the qualitative 
relations discussed below. 

From this diagram it is convenient to study 
variations of Ing with (1) p and T at constant W, 
(2) p and W at constant 7, (3) T and W at 
constant p. 
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Fic. 1. Effect of temperature, pressure, and annular 
spacing on separation factor gq. 


The first category corresponds to the problem 
of establishing optimum operating conditions in 
an existing column. The second or third would be 
involved in design considerations where the tem- 
perature or pressure were fixed within narrow 
limits by external circumstances. 

At constant W one can always improve the 
separation by decreasing the temperature and 
pressure approximately in such a manner that 
(C2T°) /(W*p*) always equals one. 

At constant temperature the separation can 
always be improved by increasing p and de- 
creasing W, subject to a similar limitation. 

Likewise, at constant pressure the separation 
can be improved by decreasing 7 and W as 
above. 

The possibility of improving the separation by 
adjusting the above variables is subject to further 
limitations, as pointed out by Furry and Jones,’ 
because of remixing due to asymmetry, remixing 
due to turbulence, and unduly long separation 
times. 
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The possibility of asymmetry in construction, 
which would cause remixing currents, increases 
with decreasing W and therefore limits the feasi- 
bility of moving indefinitely in this direction. 

The Reynolds number, which is a measure of 
turbulence, is defined by Onsager and Watson® as 
proportional to (W*p2AT/T*). This limits the 
possibility of large pressure increases except at 
very small values of W, very high temperature, or 
very low values of AT. These factors are subject 
to further limitations discussed below. 

The relaxation time is roughly proportional to 
(TL) /(W7p5(AT)*). This imposes a limitation 
on the feasibility of decreasing W, p, or AT 
indefinitely, or of operating at very high tempera- 
ture. It also limits the practical length of a 
column. These problems, including the practi- 
cality of multi-stage and continuous operation 
are discussed in detail by Furry and Jones.’ 

The following sections present the results of 
experiments performed on two single-stage batch 
columns connected to reservoirs on each end, to 
test certain of the relationships described above. 





TABLE I. Constant pressure data p= 1.118 atmos. 














Column I 
Run Av. Percent CO: 
No, Temp.°K Feed Top B.t.m. Ing 
1 382 54.41 57.72 48.38 0.376 
1 382 54.41 57.59 48.85 0.352 
2 362 54.22 56.91 49.11 0.314 
2 362 54,22 55.98 48.52 0.299 
3 334 58.34 57.83 50.86 0.281 
3 334 58.34 57.44 50.66 0.273 
4 379 57.07 59.48 49.40 0.408 
4 379 57.07 59.65 49.56 0.409 
5 356 56.89 58.10 49.73 0.337 
5 356 56.85 56.90 49.60 0.295 
6 311 55.47 56.38 51.54 0.195 
6 311 55.47 55.77 50.94 0,194 
7 334 57.57 59.76 53.95 0.237 
7 334 57.57 59.66 53.79 0.237 
8 357 53.50 56.67 49.38 0.293 
8 357 53.50 56.41 49.09 0,294 
9 398 52.68 57.45 48.10 0.376 
9 398 52.68 57.20 48.96 0.376 
Column II 





10 405 58.22 77.73 36.29 1.813 
10 405 58.22 79.44 35.06 1.968 


11 379 55.29 74.29 35.55 1.576 
11 376 55.30 73.51 36.46 1.304 
12 355 55.83 68.94 37.64 1.317 
12 353 55.84 68.82 37.17 1.320 








9L. Onsager and W. W. Watson, Phys. Rev. 56, 474 


(1939). 
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EXPERIMENTAL EQUIPMENT AND PROCEDURE 


Experiments were performed on two concentric 
tube columns each 244 cm long and connected to 
200-cc reservoirs at each end. Column I had an 
annular spacing of 0.456 cm (W=0.228) and 
column JI an annular spacing of 0.308 cm 
(W=0.154 cm). 

The outside wall was cooled by circulating 
cooling water through a jacket. The water was 
maintained at 16.5+1°C and never gained over 
1°C in passing through the column. The inner 
wall was heated by radiant energy from a -hot 
wire passed up the inside tube and held in place 
by ceramic spaghetti and spacers. The inside tube 
was sealed off and evacuated to eliminate con- 
vection currents. 

The inside wall temperatures were measured by 
means of thermocouples inserted in the inside 
tube from each end and held in place by spring 
steel rings. The whole thermocouple, except for 
the head which was protected by the ring, was 
insulated to prevent exposure to radiation. 

On the average, the top and bottom tempera- 


TABLE II. Constant temperature data T =380°K. 








Column II 
Run Pressure Percent CO: 
No. (Atmos.) Feed Top B.t.m. Ing 





1.116 56.15 74.61 37.03 1.610 
0.836 55.02 72.80 36.92 1.520 
0.845 55.02 72.42 37.34 1.483 
0.625 51.02 64.98 38.95 1.068 
0.631 51.02 64.48 38.24 1.076 
0.397 50.07 56.48 43.47 0.523 
0.390 50.07 55.68 42.36 0.536 
0.998 53.21 70.95 34.81 1.521 
0.682 51.47 65.50 36.97 1,175 
0.686 51.47 65.24 36.13 1,199 
1.410 55.16 67.23 40.20 1.112 
1.395 55.16 67.43 39.51 1.154 
0.510 51.30 62.39 42.84 0.795 
0.510 51.30 62.79 42.37 0.831 
0.744 53.63 67.34 34.64 1.359 
0.744 53.63 66.69 34.71 1.326 
1.593 57.44 67.87 42.96 1.032 
1.590 57.44 67.66 43.08 0.994 
1.248 56.48 71.03 40.08 1.299 
1.247 56.48 71.70 40.12 1.330 


HH SSCS WMA AU EI 


12 0.312 55.14 56.24 48.75 0.301 
12 0.309 55.14 56.21 49.13 0.285 
13 1.805 49.58 62.97 41.99 0.854 
13 1.805 49.58 63.01 42.40 0.839 
14 1.153 53.49 68.97 36.01 1.374 
14 1.143 53.49 68.91 35.14 1.409 
15 1.082 53.49 69.59 34.67 1.462 
15 1.082 53.49 69.51 34.96 1.445 


16 0.924 53.49 69.66 34.21 1.485 
17 0.433 51.08 58.56 43.31 0.615 
17 0.439 51.08 5840 41.67 0,676 
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Fig. 2. Comparison of experimental and predicted sepa- 
ration factor as a function of temperature and annular 
spacing. 


ture differed by 6°C with a maximum difference 
of 12°C. It is felt that this difference represents 
errors in thermocouple readings due to exposure 
to radiation rather than actual differences in wall 


temperature. 


The gases used were carbon dioxide obtained 
from the Ohio Chemical Company and 99.9 
percent propane obtained from the Matheson 
Company. The method of purification and analy- 
sis followed closely that described by Wall and 
Leaf,!° making use of a Zeiss portable interferome- 
ter for analyses, with purified propane being the 
reference gas. 

The two columns were operated over a range of 
average temperature of 311°K-405°K at 1.118 
atmospheres pressure. The average temperature 
was considered to be the arithmetic average of 
the hot and cold wall temperatures. Then column 
II was operated at pressures from 0.309-1.805 
atmospheres with a constant average tempera- 
ture of 380°K. 


10 F. T. Wall and B. Leaf, J. Phys. Chem. 46, 820 (1942). 













































The feed analysis was maintained substantially 
constant at 55 percent carbon dioxide through- 
out. The experimental results are presented in 
Tables I and II. 

Runs were continued until samples indicated 
no further change in analyses. In general, the 
first sample was not taken until at least four 
times the approximate relaxation time, as dis- 
cussed by Furry and Jones,’ had elapsed. 


COMPARISON OF THEORY AND EXPERIMENT 


The first constant pressure runs on column | 
indicated that the best value of v for the point 
source of repulsive force model was 6.15. Using 
this value and the equations of Chapman and 
Cowling,® Chapters 9 and 10, to predict a and D, 
and the data of Trautz and Kurz" for the 
viscosity of carbon dioxide-propane mixtures, 
values of Ing were calculated from Eq. (4). In 
Fig, 2 these values are compared with the experi- 
mental values of Ing taken from Table II. The 
agreement obtained for the two columns indi- 


| | | 
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Fic. 3. Comparison of experimental and predicted effect of 
pressure on the separation factor. 


"1M. Trautz and F. Kurz, Ann. d. Physik 9, 981 (1931). 
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cates definitely that the dependence of Ing on W 
is correctly predicted by Eq. (4). Also, the equa- 
tion appears to predict the correct dependence of 
Ing on the temperature over the temperature 
range covered. It should be pointed out that since 
the cold wall temperature was fixed, it was not 
possible to vary the average temperature, 7’, 
independently of AT. | 

Seventeen runs were made at different pres- 
sures using column IJ. Since these data were the 
last obtained, after operating and analytical 
techniques had been perfected, it is felt that these 
runs are the most reliable ones. 

Evaluating the constants in Eq. (4) from data 
in the literature or in the case of a and Dy: from 
theory (using the 6.15 point source model) one 
obtains an equation for Ing 


2.607 1 
‘i | (6) 
P? L1+(0.495/P*) 





Ing 


Figure 3 compares experimental values of Ing 
with those predicted from Eq. (6). It can be seen 
that the general shape of the two curves is the 
same, but that the agreement is quantitative only 
in the pressure region around 1.2 atmospheres 
where the value of »=6.15 was first established 
from experiment. Further, it can be seen that no 
adjustment of the non-pressure dependent vari- 
ables will bring the curves into agreement. 

Possible reasons for the discrepancy include: 
(1) Some consistent experimental error. (2) The 
lack of spherical symmetry of the molecules, 
whereas the equations are derived for spherically 
symmetrical molecules. (3) The use of the inverse 
power model and of the Maxwellian model for 
Dy, 9, and d. (4) The method of solution of the 
differential equation for diffusion used by Furry, 
Jones, and Onsager, first ignoring vertical diffu- 
sion and then appending a correction for it. 
(5) The possibility that v and, therefore, a may 
depend on the pressure. Considering these in 
order, no source of consistent error in the experi- 
mental work was discovered. 

The second possibility can only be tested by 
the use of spherically symmetrical molecules, 
such as the rare gases. There appears, however, to 
be no a priori reason why the lack of spherical 
symmetry should affect the pressure dependence 
at these low pressures, 
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Fic. 4. Comparison of theoretical values of Dr with those 
predicted from empirical Eq. (8). . 


Third, the inverse power model is known to be 
somewhat crude. However, the more refined 
models so far studied do not predict any different 
pressure dependence for the gas coefficients. 

The fourth possibility is not so easily disposed 
of. The fact that Eq. (4) predicts the correct 
temperature and W dependence for Ing indicates 
that its form is correct. If one considers an 
empirical equation of the form: 


a - 
ne PL + (8/P)S 


values of 6 and ¢ can be obtained from experi- 
mental values of Ing, and Eq. (7) becomes 


_ 2.607 1 
~ Pp? L14-(0.67/P4-33) J 


Equation (5) for Dr, calculated from experi- 
mental data, is then modified 


LD sv lense, = 1.761 
0.67 
pt-33 








Ing (7’) 
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Fic. 5. Comparison of values of Dr predicted from 
Eq. (5) with values. obtained from Eq. (13) when V is 
taken to be a function of pressure. 


Using the relations of Chapman and Cowling,® 
one obtains the equation 


Dr=(4.60X10-)/p, (9) 


with the 6.15 inverse power model. Figure 4 
compares values of Dr predicted from Eq. (9) 
with those calculated from Eq. (8). Since the 6.15 
model was used in obtaining a for Eq. (7’) the 
agreement is to be expected. 

The fifth possibility has no direct theoretical 
basis. If the inverse power model were really 
adequate, v would be independent of temperature 
and pressure. However, in comparing viscosity 
data with theory it is found that » must be con- 
sidered a function of temperature to obtain 
agreement between theory and experiment. 
There is evidence from viscosity investigations at 
elevated pressures that vy may increase with pres- 
sure. Chapman and Cowling® point out that the 
solid elastic spherical model (v=) gives a 
better fit for nitrogen at high pressures than at 
atmospheric pressure where v=8.8. The data of 
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Comings, Egly, and Mayland” indicate that for 
numerous gases the temperature dependence of 
the viscosity decreases with increasing pressure 
even at relatively low pressures in the neighbor- 
hood of 10-20 atmospheres. 
The temperature dependence of the viscosity is 

given by 

n2/m = (T2/T1)%, (10) 
where 

S=3+(2/y+1). (11) 


So a decreasing value of S corresponds to an 
increasing value of v and, thus, to a “hardening”’ 
of the molecular model with increasing pressure. 
If » were only slightly pressure dependent, the 
value of S would not vary noticeably with pres- 
sure at ordinary pressures, but a, which is very 
sensitive to changes in v, might vary appreciably. 

If one substitutes Eq. (3d) for a in Eq. (4), 
considering ¢(v) as a constant to be evaluated at 
a pressure where the theoretical and experi- 
mental values of Ing coincide, one obtains as an 
approximate equation for v over the pressure 


12E, W. Comings, R. S. Egly and J. C. Mayland, Uni- 
versity of Illinois Exper. Stat. Bull. 354 (1944). 
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range covered : 


5 —0.206p9-388 


(12) 





v= ° 
1—0.206p°-*33 


Using this value in the theoretical equations 
one obtains approximately 
3.98 X 1074 
Dso————-. (13) 
ps 
The values of Dr predicted from this equation 
agree closely with those calculated from experi- 
ment by Eq. (5), which becomes 


Dr=1.761p[1+ (0.495/p4) ] 
X (Ing)exper.X 10%. (14) 


This agreement is shown in Fig. 5. 

Considerably more experimental work would 
be necessary before one could put much reliance 
in Eq. (12) for » and use it to extrapolate to 
higher pressures or to other systems. If such a 
relationship is really general it would provide an 
incentive for operating thermal diffusion columns 
at slightly elevated pressure. 

The authors wish to acknowledge the en- 
couragement and helpful suggestions of Dr. F. T. 
Wall of the Noyes Chemical Laboratory. 
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The pyrolysis of toluene and the xylenes was investigated by a new technique, which made 


it possible to study these reactions for degree of decomposition ranging from 0.01 percent up 
to a few percent. It was shown that these reactions are homogeneous, unimolecular reactions. 
It was found that the weakest bond in these molecules is the C—H bond in the methyl group. 
The bond energy of that C—H bond was found to be 77.5 kcal. for toluene and meta-xylene, 
75 kcal. for para-xylene and 74 kcal. for ortho-xylene. These data indicate that the resonance 
energy of the benzyl radical is 24.5 kcal. and that the hyper-conjugation in para-xylene 


RECENT investigation of the pyrolysis 

of benzyl iodide by C. Horrex' and the 
author revealed the great stability of the benzyl 
radical. This work gave rise to the suspicion 
that the resonance energy of the benzyl radical 


1 This will shortly be published elsewhere. 


decreases the C—H bond energy in the methyl group of this molecule by 2.5-3 kcal. 





is much higher than was previously thought and 
the present study was started in the hope of 
further clarifying this point. It was decided to 
investigate the pyrolysis of toluene, and in 
particular the nature of the first step of its 
decomposition. 
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The pyrolysis of toluene and the xylenes has 
been investigated in the past ten years by several 
authors, mainly in connection with the coking 
process, and most of the research was devoted 
to studying the influence of various catalysts. 
It was found by Tilicheev? that toluene decom- 
poses 13 times faster than benzene, and the 
energy of activation of this decomposition was 
found to be 70 kcal. It is interesting to compare 
this value with the value of 77.5 kcal. obtained 
in the present work. 

Jost and Miiffling? measured the rate of 
pyrolysis of a number of hydrocarbons, among 
them toluene, following the course of the reaction 
by the increase of pressure. They observed that 
Hg(CHs3)2 had no effect on the rate of decompo- 
sition of toluene at 690°C, which indicated that 
the decomposition is not a chain reaction. All 
these investigations suffered from the common 
difficulty that the decomposition was carried out 
under conditions in which it was impossible to 
derive the true mechanism of the reaction, 
owing to the complicated mixture of products 
obtained. 

The first successful attempt to elucidate the 
mechanism of the reaction was made by Hein 
and Mésee.* These authors pyrolyzed toluene in 
the presence of mercury vapor, using a technique 

? M. D. Tilicheev, J. App. Chem. U.S.S.R. 12, 741 (1939). 

°]. W. Jost and L. v. Miiffling, Zeits. f. Elektrochemie 


47, 766 (1941). 
‘F, Hein and H. J. Mésee, Berichte 76B, 430 (1943). 














similar to that developed by Rice. They suc- 
ceeded in isolating mercury-di-benzyl as a 
product of the reaction, and followed the course 
of the decomposition by the determination of 
the quantities of the mercury compound formed. 
The rate of formation of this compound seemed 
to depend on a number of uncontrollable factors 
which prevented the determination of an energy 
of activation. Yet the authors concluded from 
the formation of mercury-di-benzyl that the first 
step in the decomposition of toluene is the 
rupture of the C—H bond, 


C.sH;-CH;—-C.H;-CH2+H (1) 


which is followed by a reaction between benzyl 
radicals and mercury atoms leading to Hg(CHe 
-CsHs)2. However, it would seem that the evi- 
dence for this conclusion is not definite, for it is 
possible to explain the formation of the benzyl 
radical by a mechanism involving the rupture of 
the C—C bond as a first step: 


CsH;-CH;—-C,H;+ CH; (12) 


followed by the reactions: 
C,H; ° CH;+CH;—-C,H; . CH.+CHs,. 
C.H; ° CH;+ C,H s—CeH, ° CH2+ CoH. 


A stronger argument, put forward by Hein and 
Mésee, in favor of (1) is their observation of 
hydrogen among the products of the decompo- 
sition. There is, however, some ambiguity about 


(13) 
(14) 









































TABLE I, 
K +1072 
Run Toluene pestialty K-10-2 sec.for % % 
No. pyrolysed T°C sec.“ 845°C He CHs 
70 once at 835°C 836° 1.50 2.18 59 41 
71. ~=‘twice at 835°C 845° 1.90 1.90 60 40 
72 once at 845°C 845° 2.02 2.02 60 40 
73 twice at 845°C 848° 2.02 1.78 59 41 


74 three times at 847°C 852° 2.12 1.76 61 39 








the origin of this hydrogen, since Hein and 
Mésee worked in the region of 1000°C, when 
according to the present study complete cracking 
of toluene molecules takes place. 

In the present investigation the rate of py- 
rolysis of toluene was observed by measuring the 
quantity and composition of the gaseous mixture 
formed. The experimental conditions were so 
chosen that the decomposition ranged from 
0.01 percent to slightly over 1 percent. It was 
expected that in these conditions the reaction 
mechanism would be simplified. The experi- 
mental results confirmed this by showing that 
no other gases than hydrogen and methane were 
produced. 


APPARATUS AND TECHNIQUE OF MEASUREMENT 


- The pyrolysis was investigated by the flow 
method in the temperature range of 680°C-850°C 
and under a pressure of 2-15 mm Hg. The 
apparatus used is shown in Fig. 1. The hydro- 
carbon under investigation was kept in a small 
vessel, V, fitted with a capillary K;. This vessel 
was weighed before and after every experiment, 
the loss of weight giving the quantity of hydro- 
carbon used. During runs V was kept in a Dewar 
vessel containing water maintained at a temper- 
ature constant to 0.1°C. By opening the tap 71, 
the vapor of the hydrocarbon was led into a 


silica reaction vessel R heated by an electrical © 


furnace. The temperature of the furnace was 
kept constant by a constant voltage transformer 
and measured by a double chromel P-Alumel 
thermocouple. The variation of temperature 
along the reaction vessel was +2°C. The pressure 
before and after the reaction vessel was measured 
on two manometers M, and M;z by means of a 
cathetometer. 

The vapor leaving the furnace passed through 
a heated tube and heated capillary Kez into a 
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small trap U cooled by mixture of salt and ice 
to —10°C, and then through a tap T? into a big 
trap H cooled by acetone +COz mixture. In 
the trap U the less volatile products such as 
di-benzyl, di-phenyl and others were condensed. 
The undecomposed hydrocarbon was condensed 
completely in the trap H. The trap H was 
connected to a system of two mercury pumps 
P; and P2, capable of pumping all the gases 
into the storage vessel S. The pumps were 
chosen so as to maintain a vacuum of 10-* mm 
of Hg in the trap H, while the pressure in the 
storage vessel S could rise up to 12 mm of 
mercury. The trap E placed between the pump 
P, and the storage vessel S was cooled with 
liquid air to trap any ethane or ethylene pro- 
duced in pyrolysis. The McLeod gauge G was 
used to measure the pressure of the gases pro- 
duced. The volume of the storage part between 
the taps 73, T,, and Ts was exactly determined 
(at about 1.5 1.). The taps 74, Ts and 77 were 
used for evacuation of appropriate parts of the 
apparatus. 7; led to an analytical device and 
Ts was used to let air into the Trap H. By 
changing the capillaries K; and Kz and choosing 
the desired temperature of water bath, it was 





TABLE II. Toluene. 











P Time 
No. toluene of 
of mm contact % K-10 % % 
Remarks run T°C Hg sec. decomp. sec.-! Hz CHa 
94 738 0.310 0.0098 0.032 
87 740 0.310 0.0114 0.037 
88 741 0.310 0.0132 0.043 
89 742 0.296 0.0141 0.048 
90 742 0.291 0.0155 0.053 60 40 
91 745 0.280 0.0163 0.058 
93 747 0.310 0.0238 0.077 


61 39 


Pressure decreased 82 773 


Time of contact 83 773 0.825 0.140 0.170 61 39 
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increased 
Surface increased 86 776 1 0.262 0.0344 0.131 60 40 
15 times 
95 779 0.294 0.0740 0.252 
78 792 0.256 0.0730 0.286 60 40 
Pressure decreased 79 796 0.286 0.107 0.375 61 39 
Time of cont. 81 803 0.905 0.437 0.485 
increased 


Time of cont. 


increased. Pres-} 80 803 13.7 0.738 0.382 0.518 59 4! 
sure increased 
97 819 6.3 0.284 0.244 0.860 
98 821 6.3 .0.284 0.261 0.920 
71 845 7.6 0.306 0.581 1.90 60 40 
72 845 8.1 0.338 0.690 2.02 60 40 
73 848 59 0.246 0.500 2.02 
Pressure decreased 75 851 2.8 0.299 0.615 2.06 59 41 
74 852 6.3 0.246 0.520 2.12 
Surface increased 84 863 8.0 0.270 1.14 4.21 62 38 
15 times : 
Surface increased 85 864 8.0 0.270 1.25 462 61 39 
15 times 
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possible to vary both the pressure of the hydro- 
carbon in the reaction vessel and the time of 
contact. 

The quantity of H.+CH, produced was given 
by the pressure in the storage part with liquid 
air on the trap E. Changing the liquid air to 
the acetone +CO, bath and measuring any 
consequent rise in pressure, it was possible to 
find the quantity of any C,-hydrocarbons 
present. ! 

Numerous experiments showed that this ar- 
rangement gives exact results even where the 
quantities of gases produced are as low as 
5.10-* cc N.T.P./minute, so that this technique 
made it possible to follow the reaction down to as 
little as 0.01 percent decomposition. The Hz to 
CH, ratio was estimated by combusting the Hg, 
in a known quantity of H2+CH, mixture, by 
passage through a CuO cell. This method was 
tested by total combustion in oxygen and esti- 
mation of the CO: formed. 


MATERIAL USED FOR PYROLYSIS 


The first series of experiments in which toluene 
samples of different origins were used, showed 
that in spite of careful purification by shaking 
with sulphuric acid, distillation and crystalliza- 
tion, some impurities remained, the pyrolysis of 
which gave rise to considerable side reactions. 
These impurities could be removed by partial 
pyrolysis and the repeated partial pyrolysis of 
toluene produced samples of toluene, which 
pyrolyse at a steady rate when partial pyrolysis 
is repeated once or several times. This is seen 
from Table I. The final rate of pyrolysis we 
take to be that of toluene itself. Accordingly, 
the further work was carried out with toluene 
purified by twice repeated partial pyrolysis at 


- 845°C, then distilled over sodium, using an 


efficient column, and collecting a fraction which 
passed within a range of 0.3°C. 


EXPERIMENTAL RESULTS OF THE PYROLYSIS OF 
TOLUENE 


The pyrolysis of toluene produced the two 
gases hydrogen and methane in the constant 
proportion 60 percent of He to 40 percent of 
CH, (see Table II). No ethane or ethylene was 
observed. In addition to the gaseous product, 
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dibenzyl was isolated in quantities corresponding 
approximately to one mol. of dibenzyl per one 
mol. of gas produced. Dibenzyl was identified by 
its melting point and mixed melting point. The 
percentage of decomposition was calculated 
under the assumption that one mol. of gas pro- 
duced corresponds to one mol. of primarily 
decomposed toluene. By packing the reaction 
vessel with silica wool which increased the 
surface about 15 times it was shown that the 
reaction is a homogeneous gas reaction (see runs 
86, 84 and 85 in Table II). Reproducibility was 
very satisfactory, as can be seen comparing runs 
89 and 90, and 71 and 72 (Table II). By changing 
the pressure of toluene and the time of contact, 
it was shown that the reaction is of the first 
order. The plot of logk against 1/T (shown on 
Fig. 2) can be satisfactorily drawn as a straight 
line, with a slope corresponding to an activation 
energy of 77.5 kcal. and a frequency factor of 
2.10", as"represented by the equation k= 2.10" 
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-e-77-5/RT, A further discussion of the grounds 
on which these constants were determined is to 
follow later. 


MECHANISM OF PYROLYSIS OF TOLUENE 


The evidence given above shows that the 
reaction is a homogeneous gas reaction of the 
first order. The fact that the frequency factor 
has the theoretical value of 10", as deduced by 
Polanyi and Wigner® suggests that we deal in 
this case with a true unimolecular reaction. 

The appearance of hydrogen as one of the 
main products forces us to assume that the 
first step of the pyrolysis is the splitting of the 
toluene molecule into a benzyl radical and a 
hydrogen atom: 


CsH;-CH;-CsH;-CH2+H (1) 


and the body of experimental evidence suggests 
that this step is followed by a sequence of quick 
reactions: 


C.sH;-CH;+H—-C,H;:CH2+H, (2) 


and 
CsH;-CH;+H—-C,.H,.+CH; (3a) 
CsH;-CH;+CH;—CH,+C,H;:-CH2 (3b) 
or 
CsH;-CH;+H—-C,H;+ CH, (4a) 
CsH;-CH3;+CsH;-CsHe+CeHsCHe (4b) 


it being assumed that the benzyl radicals are 
finally removed by dimerization to di-benzyl, 


2C.H; ° CH,—-C,H; . CH, . CH, : Ce.Hs. 


This mechanism accounts for the observed 
products (H2, CH,, dibenzyl and benzene) and 
gives the correct proportion of one mol of 
dibenzyl per one mol of gas produced. 

We can show that the recombination of 
hydrogen atoms to He is very unlikely under our 
experimental conditions, and this justifies its 
omission from the suggested mechanism. We 
do not know the energy of activation of the 
reactions (3a) or (4a), but we may expect these 
energies of activation to be smaller than that of 
(4b), which was estimated by Taylor and Smith® 

5 M. Polanyi and E. Wigner, Zeits. f. Phys. Chem. A139, 
439 (1928). See also S. Glasstone, M. Laidler and H. 
Eyring, ‘“The Theory of Rate Processes’’ (1941). 


6H. S. Taylor and J. O. Smith, J. Chem. Phys. 8, 543 
(1940). 
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at 5.6 kcal. On this basis we calculate that in 
(3a) and (4a) at least one in ten collisions should 
lead to reaction. Since the chance of three body 
collisions should be smaller than this by several 
powers of ten, the recombination of H atoms in 
the gas phase should be negligible. The recombi- 
nation at the wall should in these conditions be 
also of little importance, especially since at these 
high temperatures the accommodation coefficient 
is probably extremely small.’ 

The absence of recombination of H atoms at 
the wall was confirmed by showing that the 
packing of the reaction vessel caused no change 
in the rate of reaction, nor any increase in the 
proportion of H» produced. 

The proposed mechanism leads to the result 
that the reaction (1) is the rate determining 
step and therefore the observed energy of activa- 
tion corresponds to the energy of dissociation of 
the C—H bond in the methyl group of toluene. 
This conclusion makes use of the usual assump- 
tion that there is no energy of activation involved 
in the recombination of a bond from a: free 
radical and a free atom. While this is not the 
place to re-examine this assumption, it may be 
stated that the theoretical reasons in its favor 
seem to be very strong. 

Accordingly, we give for the C—H _ bond 
energy in the methyl group of toluene the value 
of 77.5 kcals. The difference between this value 
and the bond energy of C—H in methane which 
may be estimated at 102 kcals.® is to be ac- 
counted for by the resonance energy of the 
benzyl radical, which thus appears to have the 
value of 24.5 kcals. 


PYROLYSIS OF META-XYLENE 


The meta-xylene used was of German origin, 
produced by the ‘Gesellschaft fiir Verwehrtung 
von Kohlenteer.” It was purified by repeated 
crystallization at —60°C, which was found to 
be sufficient. 

The pyrolysis of meta-xylene produced hydro- 
gen and methane in the same proportion as that 
of toluene (60 percent He, 40 percent CH,). No 


7F. Paneth and K. Herzfeld, Zeits. f. Elektrochemie 37, 
577 (1931); F. Paneth, K. Hofeditz and A. Wunsch, J. 
Chem. Soc. (1935) p. 372. 

8 G. B. Kistiakowsky and E. B. Van Artsdalen, J. Chem. 
Phys. 12, 469 (1944). D. P. Stevenson, ibid. 10, 291 (1942). 
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TABLE III. Meta-xylene. 








P 





meta- Time 
No. xylene of 
of mm contact % K-10 % %G% 
Remarks run T°C Hg_ sec. decomp. sec.~1 He CHa 
121 740 3.0 0.358 0.037 0.104 
122 741 3.0 0.360 0.040 0.111 
123 742 3.0 0.360 0.044 0.121 60 40 
101 760 2.3 0.386 0.094 0.242 
102 765 3.5 0.352 0.086 0.242 
105 767 2.8 0.316 0.086 0.272 60 40 
Surface increased 112 778 56 0.326 0.120 0.368 62 38 
15 times 
Pressure increased } 
Time of cont. 106 802 5.9 0.840 0.825 0.98 57 43 
increased 
Time of cont. 
increased 109 812 5.5 0.891 1.29 1.45 59 41 
Pressure increased 
Time of cont. 113 815 2.5 0.910 1.37 1.50 59 41 
increased 
128 822 2.8 0.336 0.80 2.38 60 40 
129 825 2.8 0.336 0.88 2.62 60 40 
Time of contact and {108 828 5.5 0.908 2.52 2.78 
pressure increased \|108 829 5.5 0.876 2.60 2.96 60 40 
103 850 3.0 0.311 1.51 4.86 60 40 
Pressure increased 114 851 58 0.320 1.60 500 59 41 
Surface increased 110 867 4.4 0.400 3.33 8.3° 61 39 
15 times {111 867 3.7 0410 3.54 8.75 61 39 








ethane or ethylene was observed. The droplets 
of liquid collected in the U trap were identified 
by their boiling point as 3, 3’ di-methyl-di-benzyl. 

The percentage of decomposition was calcu- 
lated as before by assuming that one mol of 
each gas produced corresponds to one mol of 
meta-xylene primarily decomposed. By changing 
the time of contact and pressure of xylene the 
reaction was proved to be of the first order, and 
packing experiments confirmed it to be a homo- 
geneous gas reaction. The results are summarized 
in Table III, and the plot of logK against 1/T 
is given in Fig. 2. 

The graph of logK against 1/7 demonstrates 
strikingly that the pyrolysis of meta-xylene is 
represented by a straight line parallel to that of 


toluene but displaced towards higher values of 


logK by logarithm of 2. Accordingly, the activa- 
tion energy is taken once more at 77.5 kcals. 
while the temperature independent factor is 
chosen as 4.0-10%—that is twice the value for 
toluene. A more precise justification of these 
figures will be given later. 

The interpretation which we place on these 
facts is that the mechanism of pyrolysis of meta- 
xylene is analogous to that of toluene. The 
recurrence of the same value for the energy of 
activation indicates the same C —H bond energy; 
while the twofold value of the frequency factor 
indicates a twice increased probability of reac- 
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tion, owing to the presence of two independent 
methyl groups. 

The recurrence of the same composition of gas 
in the products of the pyrolysis of meta-xylene 
as was observed for toluene is in itself a powerful 
argument for assuming the same mechanism of 
decomposition in both cases. 


PYROLYSIS OF PARA-XYLENE 


The para-xylene used was manufactured by 
Chance and Hunt. After purification by repeated 
crystallization the substance had a sharp melting 
point at 12.8°C, and gave steady results in 
repeated pyrolysis showing that it was pure 
enough for our purposes. 

Table IV gives the experimental results, and 
the plot of logK against 1/T is to be found in 
Fig. 2. Hydrogen and methane are once more the 
only gaseous products and no ethane.or ethylene 
was observed. The gas produced is slightly richer 
in hydrogen than the gas obtained by pyrolysis 
of toluene and meta-xylene, the composition 
being 65 percent of hydrogen and 35 percent of 
methane. The packing experiments proved that 
the reaction is a homogeneous gas reaction. The 
rate of decomposition is at all temperatures 
three times greater than that of meta-xylene. 

From these results we conclude that in the 
pyrolysis of para-xylene the first step is again 
the rupture of a C—H bond of a CHs group. 
The main result clearly consists in the fact that 
para-xylene is proved to decompose three times 
faster than meta-xylene. As a further discussion 
will presently show, this can be best accounted 
for by a reduction of the activation energy 
(which is the energy of the C—H bond) to 75 


TABLE IV. Para-xylene. 











P. 
para- Time 
xylene of 
Run * mm contact % Ki@® % %&% 
Remarks No. T°C Hg_ sec.“ decomp. sec.-! He CHa 
134 745 3.2 0.419 0.158 0.378 
135 745 3.2 0.419 0.159 0.381 67 33 
130 746 3.3 0.425 0.168 0.395 
136 747 3.3 0.428 0.178 O415 67 33 
138 781 2.6 0.239 0.344 1.44 67 33 
Surface increased 133 782 6.2 0.715 0365 135 64 36 
15 times 
137 783 2.0 0.238 0.372 1.56 
139 866 3.7 0.344 7.73 224 
140 866 3.8 0.354 7.73 218 65 35 
Surface increased 131 860 4.5 0410 10.0 25.0 63 37 
15 times 
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kcal; the temperature independent factor being 
retained at the 4-10". 

The slight change in the proportion of hydro- 
gen to methane is probably the result of the 
decrease in the C—H bond energy, as this would 
increase the probability of the reaction between 
the hydrogen atom and para-xylene molecules 
which leads to the production of hydrogen 
molecules. (Reaction 2.) 

The pyrolysis of para-xylene differs from that 
of toluene and of meta-xylene in the following 
particular. The last two give rise to dimers of 
the free radicals primarily produced. The radical 
produced in the pyrolysis of para-xylene gives 
no dimer but only a high polymer which accumu- 
lates in trap U and trap H. 

The problems arising from this polymerization 
are discussed elsewhere.® It appears that the 
primarily produced para-xylyl radicals dispro- 
portionate in the gas phase according to the 
following equation : 


ig i. a. al 
+} Jo] |+ 

| | | | 

CH; CH:- CHs CH;: 


and the biradical (or its quinonoid form 


CH=—€__ >=CH:) 


polymerizes rapidly on condensation at a cooled 
surface. 


PYROLYSIS OF ORTHO-XYLENE 


The ortho-xylene used was a B.D.H. product, 
which was additionally purified by crystallization 
at —30°C, followed by a careful distillation. 
Since this product was not pure enough, it was 
further purified by partial pyrolysis. The neces- 
sity of this purification is seen from Table V. 
(Compare runs 153 and 150, 151 and 152, also 
146 and 147 with 148 and 149.) The concordant 
results obtained with pyrolyzed ortho-xylene, 
proved that they represent the pyrolysis of 
ortho-xylene itself. 


9M. Szwarc, “Discussion on labile molecules,” Trans. 
Faraday Soc., in press, 
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The main features of the pyrolysis of ortho- 
xylene are the same as those of the previously 
investigated hydrocarbons. The principal gaseous 
products are hydrogen and methane. The rate 
of their formation is somewhat (about 70 percent) 
greater than that observed for para-xylene, the 
percentage of hydrogen being smaller (50 percent 
Hz 50 percent CH,). In this case the formation 
of a small quantity of Cz hydrocarbon was 
observed ; about 10 percent of H2+CH,. 

It seems reasonable to explain the formation 
of Cz hydrocarbons from ortho-xylene by the 
presence of two adjoining methyl groups. This 
may also account for the low H2/CH, ratio if 
the steric interaction of the two CH; groups is 
assumed to decrease the Ph—C bond energy and 
consequently increase the probability of the 
reaction leading to the formation of CHg. 

These considerations leave unaffected the 
conclusion that the main process in the pyrolysis 
of ortho-xylene is again the rupture of the C—H 
bond. The slightly higher rate of decomposition 
of ortho-xylene compared to that of the para 
isomer may be taken to indicate that the C—H 
bond is weaker in the former. 

The values of logK are plotted against 1/T in 
Fig. 2, and are represented as lying on a straight 
line corresponding to an energy of activation of 
74 kcal., indicating the decrease in C—H bond 
energy of ortho-xylene as compared to meta- 
xylene. 


FURTHER DISCUSSION OF THE RESULTS 


We shall now present the more detailed con- 
sideration on the ground of which the activation 
energies and the temperature independent factors 


TABLE V. Ortho-xylene. 











ortho- Time % decomp 
No. xylene of decomp. intoC: 
of mm contact into K-10? hydro % % 
Remarks run 7°C Hg sec. H2t+CHs sec.-1 carb. He CHs 
156 730 3.5 0.370 0.140 0.379 0.010 
P 157 734 3.6 0.385 0.185 0.480 0.014 
Partially pee 155 777 3.2 0.838 0.745 2.20 0.028 
at 836° 154 788 3.0 0.340 1.06 3.12 0.045 46 54 
150 815 34 0.330 2.26 6.85 0.22 52 48 
. 151 815 3.4 0.330 2.33 7.07 0.25 
rte  — 152 815 3.5 0.339 2.30 6.90 020 52 48 
a 
Not lyzed 153 815 3.5 0.339 3.24 9.55 2.72 42 58 
Part | 148 836 4.0 0372 483 13.0 0.35 62 48 
oY ER e Ee SS 
. .00 18.4 y 
Not pyrolyzed 147 $37 35 0326 620 190 24 
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quoted above were determined. Using the method 
of least mean squares we obtain for the assumed 
linear dependence of logK on 1/T the parameters 
of the equation 


E 
i wr 


logK = 
2.3RT 


shown in Table VI. The standard deviations 
listed in the second column show that the 
accuracy is not sufficient to distinguish definitely 
between the activation energies of the reaction 
proceeding at different speeds. The situation 
can, however, be clarified further by assuming 
that the temperature independent factor is the 
same for the decomposition of one methyl group 
in each of the four investigated hydrocarbons. 
We note that all the results can be represented 
on this. basis within the experimental errors. 
Such a representation limits the possible values 
of the temperature independent factor for the 
decomposition of one methyl group to vo=2.0-10" 
—2.5-10"%. We can make use of the choice of 
any particular value within this range for a more 
accurate determination of the differences between 
the activation energies of the four pyrolytic 
reactions. 

Table VII shows the calculated values for 
energy of activation on the basis »=2.0-10"%, 
the following equations were used: 


E=2.3-RT(13.3—logK) for toluene. 
E=2.3-RT(13.6—logK) for xylenes. 


Calculations are based on experiments at low 
temperature. 

From these results we conclude that the C—H 
bond energy is the same in toluene and meta- 
xylene, while it is smaller by 2.3-2.5 kcal. in 
para-xylene and by 3.5 kcal. in ortho-xylene. 
We may recall the observations of Dobryanskii 
and Saprykin" on the relative stability of xylenes 


TABLE VI. 








Mean deviation of 





Hydrocarbon E kcal. E=4E kcal. y 
toluene 77.3 +1.3 2.10" 
meta-xylene 77.1 +1.9 4.10% 
Para-xylene 76.2 +1.5 5.10% 
ortho-xylene 74.8 +1.1 5.1018 












” A. F. Dobryanskii and F. Y. Saprykin, J. Gen. Chem. 
(USS, R.) 9, 1313 (1939). 


BOND ENERGY IN TOLUENE AND XYLENES 











TaBeE VII. 
No. of run “C E SE 
Toluene 
87 740 77.9 +0.4 
88 741 77.9 +0.4 
89 742 77.5 0.0 
90 742 77.5 0.0 
91 745 77.0 —0.4 
92 749 77.1 —0.3 
93 747 76.9 —0.5 
94 738 78.3 +0.8 
average 77.5+0.14 
Meta-xylene 
101 760 77.0 —0.2 
102 765 77.1 —0.1 
105 767 77.4 +0.2 
121 740 77.4 +0.2 
122 741 77.4 +0.2 
123 742 77.0 +0.2 
average 77.2+0.18 
Para-xylene 
130 746 74.9 +0.1 
134 745 74.9 +0.1 
135 745 74.9 +0.1 
136 747 75.0 +0.2 
137 783 74.5 - —0.3 
138 781 74.6 —0.2 
average 74.8+0.20 
Ortho-xylene 
156 730 73.9 —0.1 
157 734 73.7 —0.3 
155 777 74.9 +0.9 


average 74.0+0.55 








by pyrolytic isomerization. They found that 
meta-xylene was the most stable and did not 
isomerize, para-xylene isomerized only slightly 
to meta, while ortho-xylene isomerized very 
easily giving mainly meta- and partially para- 
xylene. These facts agree well with the sequence 
of activation energies derived above. 

The weakening of the C—H bond in the 
methyl groups of para- and ortho-xylene as 
compared with toluene and meta-xylene was to 
be expected on the assumption of hyperconjuga- 
tion occurring for example in para-xylene as 
follows: 

i re i 
H- CH i se ~*~ -H 








































136 J. 


The effect of hyperconjugation weakening in the 
C—H bond in C2H¢ as compared with that in 
CH, has been estimated at 7 kcals." A detailed 
discussion of the disparity between this value 
and that obtained for para-xylene seems pre- 
mature in view of the inaccuracy of the results 
obtained up to date. But we would like to 
emphasize that in the case of the weakening of 
C—H bonds in para-xylene as compared with 
meta-xylene we can see no other possible expla- 
nation than hyperconjugation, whereas in the 


11 E. C. Baughan, M. G. Evans, and M. Polanyi, Trans. 
Faraday Soc. 37, 377 (1941). 
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case of the weakening of the C—H bond in 
ethane as compared with methane, the proximity 
of the two interacting methyl groups offers 
additional possibilities—for example, steric ef- 
fects or dipole interactions—of an explanation. 
The author wishes to express his sincerest 
thanks to Professor M. Polanyi for his constant 
encouragement and valuable help throughout 
the work and the preparation of this paper. His 
thanks are also due to Manchester University 
for the I.C.I. Fellowship which enabled him to 
carry out this work, and to the Anglo-Iranian 
Oil Company for a grant. 
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The vacuum ultraviolet spectra of methyl and ethyl cyanide have been determined in the 


region 1000A to 1800A. Three members of a Rydberg series were found for methyl cyanide 
starting at 1292A, the limit of the series corresponding to an ionization potential of 11.96 ev. 
Vibrational band changes which have been found to accompany the first two members have 
been assigned to fundamental frequencies of the molecule by comparison with Raman data. 
The electronic transition, expressed by the Rydberg series, was ascribed to excitation of a + 
electron in the triple bond. At higher pressures, a continuous absorption appeared which was 
followed only to 1600A. No structure could be observed in the spectrum of the ethyl cyanide, 
although it paralleled that of methyl cyanide in having a stronger absorption starting at 
1342A and a weaker one appearing at higher pressures and spreading to longer wave-lengths. 


INTRODUCTION 


HE following communication contains a 
study of the absorption spectra of methyl 
and ethyl cyanides in the vacuum ultraviolet. 
Previously, Herzberg and Scheibe! studied the 
spectrum of methyl cyanide down to about 
1540A, finding a weak band extending from 
1820A to 1600A with a maximum at about 
1670A. Below 1600A, there was an indication of 
the beginning of a stronger absorption which 
would extend to shorter wave-lengths. They also 
* Taken from a dissertation submitted to the faculty of 
the Graduate School of New York University in partial 
fulfillment of the requirements for the degree of Doctor of 
Philosophy. 
** Present address: Polytechnic Institute of Brooklyn, 
Polymer Research Institute, Brooklyn, New York. 


1G, Herzberg and G. Scheibe, Zeits. f. physik. Chemie 
B7, 390 (1930). 





observed the same indications in the spectra of 
the methyl halides. Price? extended this work on 
the halides to shorter wave-lengths and found 
discrete Rydberg series, assigning the excitation 
to the pz electrons on the halide. From this it 
was felt that if the work on methyl cyanide were 
also extended to shorter wave-lengths, a Rydberg 
series might be identified. 


EXPERIMENTAL 


The far ultraviolet spectra were photographed 
on a Lyman continuum as a background, using 
hydrogen as the inert gas, by means of a vacuum 
spectrograph at nearly normal incidence. A lightly 
ruled glass grating was used, of focal length 1 
meter, with 30,000 lines per inch, ruled by R. W. 


2 W. C. Price, J. Chem. Phys. 4, 539 (1936). 
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ALKYL CYANIDE SPECTRA 








Fic. 1. Spectrum of 








methyl cyanide. 








Wood. The dispersion was about 8.43A per mm. 
Eastman Kodak III-O, ultraviolet sensitized 
plates and Ilford Q1 plates were used. The 
Eastman plates were made on x thin glass so 
that they could be bent to fit the Rowland 
circle. Three exposures could be made on one 
plate, each exposure covering 500A. Discrete 
bands were measured from microphotometer 
traces using the strong Si and O lines in the 
source as standards. ‘ihe spectrum of Mel was 
taken as a test plate for the region 1900-1300A 
since Price® has published an illustration of the 
sharp bands found in that region. The absorption 
spectra of the compounds were taken with the 
vapor filling the entire spectrograph, thus pro- 
ducing a total absorbing path of 2 meters. The 
concentration of the vapor was controlled by the 
temperature of an acetone bath cooled by dry 
ice, surrounding the distilling flask leading into 
the spectrograph. 


METHYL CYANIDE 


The material supplied by the Eastman Kodak 
Company was dried for several days over calcium 
chloride and then distilled in air in an all glass 
still. The main fraction came over from 81 to 
81.5°C. The methyl cyanide was immediately 
put into a distilling apparatus connected to the 
spectrograph and redistilled in vacuum before 
being vaporized into the spectrograph. 

Typical spectra of methyl cyanide are repro- 
duced in Fig. 1 where they are enlarged about 
three times. This compound has a set of strong, 
discrete bands which start abruptly at 1292A at 
the lower vapor pressures studied. As the pres- 
sure is increased these discrete bands widen, and 
a continuous absorption appears to spread from 
the short wave-length end. The continuum is 
believed to be caused by the absorption of 
methane and ethane in this region, starting at 
1300A and 1350A, respectively.* 4 As the pressure 
is increased, the absorption no longer starts 


ED 


(934) F. Duncan and J. P. Howe, J. Chem. Phys. 2, 851 


‘W. C. Price, Phys. Rev. 47, 444 (1935). 








abruptly at 1290A, but a continuous absorption 
occurs which travels to the longer wave-lengths. 
At the highest temperature used, —25°C, this 
continuous absorption extends to 1600A. The 
weak band from 1820A to 1600A reported by 
Herzberg and Scheibe was not observed, probably 
because the vapor pressure was never high 
enough. 

The spectrum of methyl cyanide (see Figs. 1 
and 2) is not an ideal one for measurement and 
interpretation. The bands are not very sharp, and 
it is impossible to find indications of structure 
which might facilitate in the analysis. In this 
respect the spectrum of methyl cyanide resembles 
the shorter wave-length spectrum of HCN as 
reported by Price.** The overlapping continua 
extending from the shorter wave-lengths are a 
serious difficulty, obscuring the shorter members 
of the discrete banded systems. Nevertheless, it 
seems very probable that a Rydberg series is 
represented by the strongest bands 1, 7, and 10. 
These bands can be represented by the formula 


vo" = 96,426 —[R/(n+0.40)?]. 


The constants were chosen to fit the first member, 
n=2, since the measurement of the member 
corresponding to »=3 is made difficult by inter- 
fering O III emission lines, and the band corre- 
sponding to =4 is quite broad. The agreement 
between observed and calculated values is shown 
in Table I. The convergence limit, 96,426 cm- 
=1035A, is obscured in the C—C and C—H 
continua ; the corresponding ionization potential 
is 11.96 volts. 

The strong bands 1 and 7 are accompanied on 











TABLE I. 
Observed cm= Calculated cm=1 
n=2 77,374 77,374 
n=3 86,953 86,933 
n=4 90,853 90,757 








5 W. C. Price, Phys. Rev. 46, 529 (1934). 
6 W. C. Price, and A. D. Walsh, Trans. Faraday Soc. 41, 
381 (1945). 
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Fic. 2. Microphotometer tracing of methyl cyanide spectrum. 


the short wave-length side by weaker bands 
which represent transitions to vibrating elec- 
tronic states. The differences of the main bands 
of the 1292A group are listed in Table II. 

There seems little doubt that the strongest 
vibrational band associated with the Rydberg 
transition of longest wave-length, with a fre- 
quency difference of 2001 cm from the O-O 
band, originates in the excitation of the C=N 
valence vibration (v2),*** of frequency 2249 in 
the normal state, as determined by the Raman 
spectrum.’ The diminution in frequency in the 
excited state from its value in the normal state, 
indicating a diminution of bonding strength, is a 
common occurrence, being analogous, for ex- 
ample, to the change in the symmetrical C=C 
vibration in ethylene from its value of 1623 cm! 
in the ground state to 1370 cm in its first 











TABLE II. 
Differences from 
Band cm~1 0-0 
1,00 77,374 
2 78,187 813 
3 78,342 968 
4 78,613 1239 
4a 78,934 1560 
5 79,375 2001 
6 80,098 2724 








*** Notation in G. Herzberg, Infrared and Raman Spectra 
of eet Molecules (D. Van Nostrand Company, Inc., 
New York, 1945), p. 333. 

(19 A W. Reitz and R. Skrabal, Monats. F. Chem. 70, 398 


Rydberg excited state.* Similarly, the band dis- 
placed by 2724 cm from the O-O band is 
probably the symmetrical C—H vibration of 
frequency 2942 cm-! (y;) in the ground state. 
The band displaced 1239 cm- from the O-0 
band may represent the symmetrical deformation 
frequency of the methyl group, 1376 cm in the 
ground state (v3). The ground state has a totally 
symmetrical ““C—C”’ frequency v4, at either 916 
cm or 1124 cm~",*** more likely the former, and 
it therefore seems probable that the frequency 
interval 813 cm or 968 cm-, more likely the 
former, represents this vibration in the excited 
state. Herzberg assigns the second of the fre- 
quencies just mentioned in the normal state to 
the degenerate rocking frequency of the methy! 
group (v7), and it is possible that the 968 cm” 
band represents this frequency in the excited 
state, although the intensity seems unusually 
large for a degenerate vibration. The very weak 
band 4a (Fig. 2), with a frequency separation of 
1560 cm— from the OO band, may be 2y,, and 
a weak band at about 4000 cm- from the O-0 
band may be 2». 

The vibrational transitions associated with the 
second Rydberg member at 1150A are less 
certain because of difficulties of measurements 
imposed by interfering emission lines in the 
source. A strong band displaced 2009 cm™ 
(band 9) from the O-—O band doubtless repre- 


8 W. C. Price and W.T. Tutte, Proc. Roy. Soc. London 
A174, 207 (1940). 
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sents the C=N vibration (v2), and the frequency 
interval of 1107 cm (band 8) may represent 
the methyl deformation vibration (v3). 


ETHYL CYANIDE 


No discrete bands could be found for ethyl 
cyanide. A diffuse absorption starts at 1342A 
and extends to shorter wave-lengths. This ab- 
sorption occurs at low pressures similar to those 
at which the Rydberg type of absorption was 
predominant in methyl cyanide. At higher pres- 
sures the continuous absorption extends to longer 
wave-lengths similar to the continuous absorp- 
tion of methyl cyanide. 


DISCUSSION 


The discrete Rydberg band series starting at 
1292A in the case of methyl cyanide is assigned 
to an excitation of one of the four m electrons 
making up the C=N bond. The fact that the 
C=N vibration is most strongly excited lends 
support to this assumption. This excited elec- 
tron, however, can set up methyl and C—C 
vibrations since its excited orbital will probably 
extend, to some extent, into the rest of the 
molecule. As in the case of methyl iodide (where 
the excited electron has been proved to be the pr 
electron on the iodine),? symmetrical methyl 
vibrations in another part of the molecule accom- 
pany the early Rydberg members. Since only 
the O—+1 vibrations are excited strongly in each 
Rydberg transition, we may conclude that the 
relative positions of the atoms in the molecule 
do not change very much by excitation. 

It is interesting to note that one of the Rydberg 
series found by Price and Walsh® for methyl- 
acetylene, 


vo" = 91,240 —[R/(n+0.48)*], 


has a constant in the term value denominator 
(0.48) which is similar to that found for methyl 
cyanide. Price and Walsh showed that it was a + 
electron from the triple bond of methylacetylene 
which was excited. The term values observed 
for methylacetylene and methyl cyanide are 
compared in Table III. 

Although the term value for the ground state 
in methylacetylene is 91,240 cm as compared 
with 96,430 cm- for methyl cyanide, the term 
values for the corresponding higher states of -the 


ALKYL CYANIDE SPECTRA 








TABLE III. 








Methyl cyanide (obs.) Methylacetylene (obs.) 





n=2 19,052 17,643 
n=3 9,473 9,003 
n=4 5,573 5,488 








mw electron of the triple bond are very similar. 
The lowest members are the ones most affected 
by the immediate configuration of the surround- 
ing molecule. Since the two molecules are iso- 
steric, it is interesting to note the similarity in 
term values and one may probably conclude that 
a m electron of the triple bond is also excited in 
methyl cyanide. 

The ionization potential of HCN is given by 
Sponer as 14.8 ev,® as compared with 11.96 ev 
found for methyl cyanide in these experiments. 
This reduction is similar to that found by Price 
and Walsh® for methylacetylene as compared 
with acetylene, and can be accounted for in the 
same way, by an inductive effect and hyper- 
conjugation. The inductive effect transfers nega- 
tive charge into the triple bond from the 
polarizable methyl group. In the case of hyper- 
conjugation, the three C—H bonds behave some- 
what as a triple bond and resonance occurs as in 
ordinary conjugation, resulting in putting more 
charge into the triple bond. This conjugation 
can be represented as: 


+ _ 
H;=C—C=N=-H;=C=C=N 


or 
+ 
H “ 
HC=Ce=N . 
H 


Hurdis and Smyth” explained the increase in 
the dipole moment of methyl cyanide (3.94 
Debyes) over that of hydrogen cyanide (2.93 
Debyes) on the basis of the hyperconjugation 
and inductive effect in methyl cyanide. Both 
effects would put extra negative charge into the 
C=N and, therefore, the ionization potential 
would be lowered. 

There seems to be some difficulty reconciling 

; ; : : 
1935)" we Molekulspektren (Julius Springer, Berlin, 

10 FE. C. Hurdis and C. P. Smyth, J. Am. Chem. Soc. 65, 
89 (1943) ; N. B. Hannay and C. P. Smyth, paper presented 


at the American Chemical Society meeting in Atlantic City 
on April 9, 1946. 


140 i I. 


the value of the ionization potential of methyl 
cyanide as deduced in these experiments, 11.96 
ev, with the value of about 13.8 ev deduced for 
cyanogen by Price and Walsh® from the diffi- 
cultly measurable absorption spectrum of that 
compound. Methyl cyanide bears the same rela- 
tionship to cyanogen as methylacetylene to 
diacetylene. The ionization potential of methy]l- 
acetylene, deduced from a well established Ryd- 
berg series by Price and Walsh, is 11.36 ev, and 
that of diacetylene, similarly obtained, is 10.74 
ev. This is the relationship to be expected since 
resonance in the normal conjugated system 
represented by diacetylene will lower the ioniza- 
tion potential more than in the hyperconjugated 
system, methylacetylene. Similarly, the ioniza- 
tion potential of cyanogen is expected to be 
lower than the 11.96 ev of methyl cyanide, which 
in turn would be lower than that of hydrogen 
cyanide, 14.8 ev. 

The first Rydberg type band in cyanogen is 
at 1320A as compared with 1292A in methyl 
cyanide, and this suggests the possibility that 
the first ionization potential of cyanogen should 
be lower than that of methy! cyanide. Similarly, 
the first Rydberg band of diacetylene is at 
1630A (I.P.=10.74 ev) while the first Rydberg 
band of methylacetylene occurs at 1540A (I.P. 
= 11.26 ev). If Price and Walsh’s value for the 
ionization potential of cyanogen is taken, then 


A. CUTLER 


in the light of this experiment, the resonance 


-effects among cyanogen, methyl cyanide, and 


hydrogen cyanide are different from those among 
diacetylene, methylacetylene, and acetylene. 
The beginning of the absorption of ethy! 
cyanide at 1342A may represent the limit of the 
C—C absorption or it may correspond to the 
first Rydberg member of methyl cyanide. Its 
shift to longer wave-lengths would indicate a 
lowering in ionization potential of a few tenths 
of an electron volt. Since the actual members of 
a Rydberg series could not be measured, this 
evaluation is very rough. The lowering of the 
first ionization potential of ethyl cyanide com- 
pared to methyl cyanide may be ascribed to the 
increase in induction of negative charge into the 
triple bond from the ethyl group as compared 
to the methyl group. The dipole moments are 
4.00 Debyes and 3.94 Debyes respectively.” The 
lack of structure in the absorption of ethyl 
cyanide may be due to the relatively stronger 
C—H and C—C continua (more C—H and C—C 
bonds) or low frequency vibrations may be 
present giving the appearance of a quasi-con- 


‘tinuum. 


In conclusion, it is a pleasure to acknowledge 
my indebtedness to Dr. William West for his 
supervision and encouragement during the course 
of this work and for his helpful suggestions on 
the preparation of this paper. 
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The quantum correction has been expressed as a deviation 7 from InQ., where Q, is the 


energy, specific heat, etc., caused by quantiza- 


tion can be expressed in terms of n and certain thermodynamic quantities. 7 itself is expressed 
in terms of thermodynamic quantities. The equations obtained give for solid neon the same 
results, within the limits of error, as the usual treatment of a Debye solid with the aid of the 
Griineisen equation. To apply the equation to liquid neon, various thermodynamic quantities 
have been obtained from known values for argon by use of the law of corresponding states. 
The difference in volume and specific heat of liquid neon from the values expected from the 
law of corresponding states have been computed and found to agree reasonably well with the 
observations. However, the melting point, and the AE and AS of fusion are found experi- 
mentally to be much closer to the classical values than one would expect from the theory. 


No explanation of this discrepancy is apparent. 





T has frequently been noted that the properties 

of neon are somewhat out of line with those of 
the other condensed rare gases. These deviations 
have been ascribed to quantum effects in neon,? 
which, it is believed, would otherwise obey the 
law of corresponding states. Quantum effects in 
the solid are easily taken care of with the help of 
the Debye theory ; the liquid presents somewhat 
more of a problem. 

The quantum corrections for a liquid have re- 
cently been given in a general form by Mayer and 
Band.4 Their results, however, involve some 
multiple integrals containing distribution func- 
tions and second derivatives of the interatomic 
potential function. They are, therefore, difficult 
to evaluate, and it seems desirable to derive some 
rougher equations which will be more easily 
handled. It will be the purpose of this paper to 
express the magnitude of the quantum effects in 
terms of thermodynamic quantities which have 
either been measured or may be approximated 
from the law of corresponding states, and apply 
the results to neon. 





1This document is based on work performed under 
Contract No. W-35-058-eng-71 for the Atomic Energy 
Project at the Clinton Laboratories. 

On leave 1946-47 from the University of North 
Carolina, Chapel Hill, North Carolina. 

*See K. S. Pitzer, J. Chem. Phys. 7, 583 (1939). 

‘J. E. Mayerand W. Band, J. Chem. Phys. 15, 141 (1947). 
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I. EQUATIONS FOR A ONE-PHASE ISOTHERMAL 
AND ISOBARIC SYSTEM 


We shall express the quantum deviation by 
means of a quantity » defined by the equation® 


InQ=InQ.+7. (1) 


Here Q is the actual partition function, and Q, is 
the classical partition function of the system, at 
the given temperature T and molal volume V. Q, 
is found by taking the limiting form of Q for high 
temperatures, and inserting the actual value of T 
in this limiting form. In calculating Q only the 
thermal part of the energy, Z;, and the zero-point 
energy, E,, are considered. The Helmholtz free 
energy is therefore given by 


A=-—RT InQ+E,, (2) 


where E, is the potential energy (excluding any 
potential energy involved in E; and £,) ; that is, 
E, is the value the energy would have with all the 
atoms at rest in their positions of equilibrium.® 
From Eq. (1) this may be written 


A=-—RT InQ.—RTn+E£, (3) 
=A.—RTn. 


5 If » is very small we may write Eq.,(1) in the form 
Q=Q.(1+7), which brings out the significance of 7. How- 
ever, to actually define 7 by this equation would be less 
convenient for purposes of computation than to use Eq. (1). 

6 For a solid the positions of equilibrium form a regular 
lattice, but this is, of course, not true for a liquid. 
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In this equation we assume either that E, is the 
same for the classical and quantum case if the 
volume remains fixed, or, better, that any change 
is included in 7. Since E, is the energy after 
subtracting out all thermal and zero-point energy, 
it would appear that any change in it which could 
be caused by quantization could be referred to a 
change in the average coordination number 
caused by quantization. Such a change could 
occur in the liquid, but may be included in 7. 
This, then, eventually becomes merged with the 
assumptions regarding 7 which are later made, in 
order to estimate it for the liquid. 

The Gibbs free energy F is equal to the 
Helmholtz free energy if the pressure p=0, as is 
essentially the case for condensed systems around 
atmospheric pressure. It is to be noted, however, 
that, in Eq. (3), A. is to be evaluated at the same 
molal volume as A, not the same pressure. If A is 
evaluated at zero pressure, A, must be evaluated 
at that pressure, 6p, at which the system would 
have the same volume if the classical statistics 
held. 


Since the pressure is given by 
p= —(0A/dV)r; (4) 


and since 6p is the pressure the system would 
have if Q. were the correct partition function, we 


may write 
5p = —(0A./dV)r. (5) 


The true pressure being zero, we obtain from 
Eqs. (3), (4), and (5) 


6p = —RT(0n/0V)r. (6) 


If we let A.,o be the classical value of A at zero 
pressure, we have’ 


op 
| A=Aqot [ (aA ./ap) ndp 


bp 
=A. +f (0A ./dV)7r(dV/dp)r, dp. 
0 


Using Eq. (5), and integrating, assuming 
(9V/dp)7,- to be constant, 


Ac=Ac,o— (8p/2)%(8V/0p) 1, 
Aside from the assumption that (0V/0p)7,. is 


7 We use, in general, a subscript ¢ to designate the clas- 
sical value of any quantity. 
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constant, this equation is exact. If 7 and hence 6p 
is a small quantity, it shows that A.-—Az< o is a 
small quantity of the second order. Such quanti- 
ties may, in general, be neglected in the case of 
the liquid, but may not be entirely negligible for 
the solid. The form of the dependence of A,.—A. 4 
on 6p is that which is to be expected from the fact 
that for p=0 the condition of equilibrium is that 
A be a minimum. 

There will be a first-order effect on the volume 
of the system. Let V, be the volume the system 
would have at zero pressure if the classical sta- 
tistics held, and let V be its actual volume. A 
classical system at volume V would have the 
pressure 6p. To go to zero pressure its volume 
must change by —(0V/0p)7,.6p. We therefore 
have 


V—V.=(0V/0p)7, 6p 
= —RT(dn/dV)r(0V/dp)7,c. (7) 


There will also be a first-order effect on the 
energy. We have, in general, for the sum of the 
thermal energy and zero-point energy 


E.+E,=RT?(d InQ/dT)y. (8) 
At a given volume V, therefore, 


E,+E,—E1,.=RT{(a InQ/aT)y 
— (0 InQ./dT)y ]=RT*(dn/dT)y, (9) 


where £,,, is the classical thermal energy at the 
volume V. As noted above, quantization causes 
no change in £,, since any effect is taken care of 
by 7. However, at volume V the pressure of a 
classical system would, as we have seen, be 4p. 
Going back to zero pressure the energy would 
change by —(dE/dV)7,.(0V/dp)7,-6p. Further, 
at zero pressure we may write (dE/dV)r,- 
=T(0S/0V)7,-=T(0p/dT)yv,¢. Therefore, for fixed 
zero pressure we may write (using Eq. (6)) 


E—E,=RT*(dn/0T)y —RT2(dn/8V) 17 
X (0p/8T)y, (8V/dp) r, 
= RT*[(dn/8T)y 
+(8n/AV)1(AV/AT) y,c]. (10) 


Provided the differential coefficients are constant, 
both Eqs. (7) and (10) are exact. 

The quantum contribution to the specific 
heat, 65C,, can be obtained by differentiating 
Eq. (10). We note that (0/8T),=(0/0T)v 
+(AV/dT),(8/AV)r. This expression is to be 


~, — Af ~*~ ~~ ma — — Ae 
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applied to the true quantized system rather than 
a hypothetical classical one, so the (8V/dT) >, 
which appears here (without subscript c) is the 
actual, rather than a classical value. However, 
the operator, (0/07), or (0/0T)y and (0/dV)r, 
applied to a classical expression gives a classical 
expression. We therefore obtain 


6Cp=2RT[(0n/dT)y + (0n/dV)7(0V/0T) pc] 
+RT*[(0°n/dT?) y+ (0°n/dTAV) 
X {(0V/8T) p+ (dV/AT) pc} 
+ (8°n/V?)7(8V/dT) ,(0V/AT) pc 
+(dn/9V)7(d°V/dT) p<]. (11) 


In the case of liquid neon, (@V/0T), and 
(0V/AT) p,- will be nearly enough alike so it will 
not be necessary to distinguish between them. 


II. THE EVALUATION OF THE QUANTUM 
CORRECTION, n 


The evaluation of the quantum correction for 
the solid is very easy, if we assume that it obeys 
the Debye specific heat law. In this case we can 
write the expression for InQ in the form® 


InQ=In[ (7/0) *e ]—(3/40)(9/T)? 
+ (3/6720)(@/T)*—-++, (12) 


where 9 is the Debye characteristic temperature. 
If we break off Eq. (12) at the second term and 
compare with Eq. (1) we see that this gives 


n= —(3/40)e10.-?= —0.14640,-. (13) 


In many respects the liquid resembles the 
solid, and we shall assume that Eq. (13) holds for 
the liquid as well as the solid. It is true that there 
is a translational, as well as a vibrational, con- 
tribution to the partition function of the liquid, 
and there may be some complications arising 
from the possibility of a change in coordination 
number in the liquid, as indicated in the discus- 
sion following Eq. (3). However, the quantum 
correction will, in general, be expected to de- 
crease as the partition function increases, and in 
view of the reasonably close similarity between 
liquid and solid, it would appear that Eq. (13) 
cannot be too far incorrect for the liquid. 

We may now refer back to Eq. (2) to get 7 in 
terms of thermodynamic quantities. Since A.-— Ep 
=E.—TS.—E,=E1<.—TS., where Ez. is the 


5 J. E. Mayer and M. G. Mayer, Statistical Mechanics 
(John Wiley and Sons, Inc., New York, 1940), p. 255. 
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classical thermal energy, we may write from 
Eq. (13) 


n= —0.146 exp(2E,,-/3RT—2S./3R). (14) 


Since (0E,,./8T)y =Cy,- and (0S./d8T)y =C,,./T, 
we have 


(0n/8T) y = —(2/3)nE,,-/RT?. 
We also have 


(0n/0 V) 7 = (2/3)n[ (OE: -/8V)7/RT 
—(0S./0V)r/R]. (16) 


We note that 7(0S./0V)7=(0E./8V)r+p. Sup- 
stituting this into Eq. (16), and remembering 
that E.=E,-+E, (the zero-point energy does 
not occur in the classical case) we get 


(dn/d V) , 
= —(2/3)n(dE,/dV)/RT —(2/3)np/RT. 


We will, in general, apply this with p=0. 

To apply Eq. (11) we also need the second 
derivatives. In the final expressions for the 
second derivatives we set p=0. 


(0°n/dT*) y = (4/9) nE:,.?/R°T' 
“7 (2/3)nCy, e/ RT? 
+(4/3)nE.,-/RT*, 
d*n/dVOT = (4/9)nE1,-(dE,/dV)/R?T® 
— (2/3)n(OE,, -/8V)r/RT?, 
(0°n/d V?) 7 = (4/9)n(dE,/dV)?/R*T° 
— (2/3)n(@?E,/dV?)/RT 
—(2/3)n(dp/dV) 7, -/RT. 


(15) 


(17) 


(18) 


(19) 


(20) 


Equation (19) can be somewhat transformed. 
Using the relation at zero pressure, 


(0E./0V)7=(0E.,-/8V)r+dE,/dV 
=T(0S./0V)r=T(dp/dT)v,., 
we have 
0°n/d VOT = (4/9) nE,, .(dE,/dV)/R?T? 
+(2/3)n(dE,/dV)/RT? 
—(2/3)n(0p/dT)y,-/RT. (21) 


Ill. APPLICATION TO SOLID NEON 


Application of the above results to a Debye 
solid should yield nothing essentially different 
from application of the Griineisen equation,® 


p=T(E.+ E,) —dE,/dV, (22) 


9 See O. K. Rice, J. Chem. Phys. 12, 289 (1944), 
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where 
r=—0-d0/dV. (23) 


However,-by way of illustration, it will be of some 
interest to apply Eqs. (7) and (11) to solid neon. 
In this case, we shall not use Eq. (15) or (17), but 
will go directly to the relation which follows from 
Eq. (12), 


= — (3/40) (O/T)?+ (3/6720)(O/T)4. (24) 
From this we obtain 


(dn/d V)r=T'(6/40) (O/T)? 
— (12/6720)(@/T)*] 
&—2Fn[1—0.0060(0/T)2]. (25) 


We may also write 
(0V/dp)7,-= —(9V/0T)»,-/(AP/AT)v,c. (26) 


(0p/dT)y,- is readily obtained by differentiating 
Eq. (22), taken for the classical case. Since T and 
dE,/dV are independent of T this gives 


(0p/8T)v,e= Cy, I =3RP. (27) 
(8V/8T)»,- may be written in the form 
(8V/8T) ».-= —I“T—(8 In@/a InT),,-. (28) 
Using Eqs. (25) to (28) in Eq. (7), we obtain 


T'(V— V.) = (2/3)(0 In@/d InT) », 
xX [1—0.0060(0/T)?]. (29) 


If we have a series of substances which obey 
the law of corresponding states, the value of 
(d In@/d InT),,- will be the same for all of them.’ 
In the case of the rare gases xenon may be con- 
sidered as practically classical at its melting 
point. The melting points are very near to being 
corresponding states; hence we may take the 
value of (@1n@/d InT), given by Murphy and 
Rice for xenon at its melting point" as the 
(0 In@/d In7),,. for neon at its melting point. 
If we set 9, equal to the value of 0 at V., then 
it is seen that the left-hand side of Eq. (29) will 
be —(@—9.)/Ow, where O,, may be taken as the 
mean of © and ©,. © at the melting point of Ne 
(or rather the temperature of the corresponding 


10 See G. W. Murphy and O. K. Rice, J. Chem. Phys. 14, 
518 (1946). 

11 We will use the calculations based on Murphy and 
Rice’s “Curve II.” These calculations do not fit the experi- 
mental data on neon and xenon perfectly, but they are 
reasonably good and they do exactly fit the law of corre- 
sponding states as extended by Murphy and Rice. 
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point, based on the melting point of xenon, which 
differs by only a few tenths of a degree) may be 
found by a slight extrapolation in Table V of 
Murphy and Rice, while ©, may be estimated 
from the results with xenon. We find 


—(O-—0.)/Ow 
= — (59.2 —65.8)/62.5=0.106. (30) 


Evaluating the right-hand side of Eq. (29), using 
@=59.2 and T = 24.9 to calculate yn, we obtain 


(2/3) (0 In@/d InT) », nl 1—0.0060(0/T)?] 
= 0.667 X 0.385 X 0.410 X0.965 =0.101. (31) 


This is an excellent check, especially since 
(8V/dp)7,- will actually vary with volume, but is 
evaluated at V, (this is effectively what is done, 
since the right-hand side of Eq. (28) is evaluated 
at V. by the above procedure). 

The value of T for argon has been taken as 0.18 
mole cc. For neon it can be calculated by 
multiplying by the inverse ratio of the critical 
volumes,!? according to the law of corresponding 
states; it therefore is 0.325 mole cc. This gives 
V—V.=0.33 cc/mole from Eq. (30). The actual 
volume of solid neon given by Clusius'* is 13.98 
cc/mole. If we multiply the volume of solid 
krypton" at its melting point by the ratio of the 
critical volumes we get 13.42 cc/mole for V.. 
Similarly, from argon’ we get 13.59 cc/mole for 
V.. The calculated value of V— V, seems to agree 
about as well as one might expect. 

We have also evaluated Eq. (11) for the solid. 
This requires some further differentiation, in- 
cluding differentiation of Eq. (28), which brings 
in the quantity (0? 1nQ/d InTdT),,., which has 
been estimated as —0.00405 X (289.8/44.8) from 
Table II, Murphy and Rice (the ratio in paren- 
theses is the ratio of critical temperatures of 
xenon and neon, taken from Guggenheim). We 
also have to evaluate (9V/dT), (non-classical) 


. which is done with the aid of Table V of Murphy 


and Rice. We have obtained for 6C, at the 
melting point —2.0 cal. mole deg.—!. Compari- 
son of Tables II and V of Murphy and Rice shows 
that this is about right. 


2 E. A. Guggenheim, J. Chem. Phys. 13, 253 (1945). 
13 K. Clusius, Zeits. f. physik. Chemie B31, 459 (1936). 
14K, Clusius and K. Weigand, Zeits. f. physik. Chemie 
B46, 1 (1940). 
15 Q, K. Rice, J. Chem. Phys. 14, 321 (1946). 
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IV. APPLICATION TO LIQUID NEON 


The checks obtained in Section III merely 
confirm the consistency of the present calculation 
with that of Murphy and Rice. The comparison 
of experimental and calculated V— V, may have 
some significance, but again the calculation may 
be handled equally well by either method. For 
the case of the liquid, however, only the equa- 
tions of the present paper are available. 

To calculate V— V., or better (V—V.)/V., we 
clearly need, from Eqs. (7), (17), and (26), values 
of V., dE,/dV, (8V/dT)>»,-, and (dp/dT)y,.. To 
get 6C, we need also Ez,., Cv, c, (0?V/dT?) pc, and 
d’?E,/dV*. These quantities are not known for 
neon; they are known or may be readily calcu- 
lated for argon from related quantities.!® Since 
the quantities appearing in the expressions for 
(V-—V.)/V- and 6C, are, except for 7, in a 
combination which is independent of the sub- 
stance if the law of corresponding states holds, 
the values for argon should be sufficiently good. 
The values for the second derivatives show a 
rather sudden change at a molar volume equal to 
the molar volume at the normal melting point of 
argon. However, the volumes we are interested in 
are equal to or greater than the normal melting 
volume, so no ambiguity occurs here. The second 
derivatives can be obtained by differences from 
the first derivatives; the values are necessarily 
rough, but they fortunately do not contribute 
greatly to 6C,. d?E,/dV? is sufficiently small, and 
the terms containing it are sufficiently unim- 
portant, so that they can be neglected. 

We have yet to consider 7 for neon. This is to 
be obtained from Eq. (13). Q, can be obtained by 
comparison with argon. We have previously 
shown that Q, can be expressed approximately in 
the form!é 


Q.= Qo? =[(2e)/*T/O 


+ (2emkT /h2)¥22-1°’ ~ ae) 8. (32) 


On going from one substance to another which 
obeys the law of corresponding states, the two 
terms in the bracket transform in the same way. 
(For the transformation of © see Murphy and 
Rice.” To convert © from argon to neon we use 
fe (argon) = 1.530 and fg (neon) = 1.405.) Qo has 
been calculated for argon in connection with the 





*O. K. Rice, J. Chem. Phys. 14, 324 (1946); see 
especially Tables V and VI. : 
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TABLE I. Thermodynamic quantities for argon; 7 for argon 
and neon. 








V=28.02 cm’ mole 
T (fusion) = 83.78°K 
dE,/dV =50.5 cal. cm=* 
E./3RT =0.7961 
C,/3R=0.856 
(8V/dT),=0.1178 cm? mole deg. 
(0p/8T),=0.556 cal. cm~* deg.“ 
(@?V /dT?),=0.00165 cm* mole deg.~? 
Qo (argon) = 3.681 
—7 (argon) =0.01080 
—7 (neon) =0.1030 








calculations of the article cited, whence Q» for 
neon and 7 for neon can be obtained as indicated. 
The various thermodynamic quantities required 
and 7 are listed in Table I. 

This table gives the values of the various 
quantities at the melting temperature and 
melting volume of argon, or, in the case of 7, at 
the corresponding point for neon. Neon at its 
melting point (24.55°K), or at its corresponding 
temperature (24.9°, which is the temperature we 
have actually used), is slightly expanded beyond 
its corresponding volume. The errors which will 
be caused by this discrepancy are not sufficient to 
make it worth while to perform further calcula- 
tions. It is easy to find from reference 16 how the 
various quantities involved change with volume 
at constant pressure. This should give the order 
of magnitude of the errors involved, and it is found 
to be of the order of 10 percent in (V—V.)/V. 
and 15 percent in 6C,. However, these are 
probably exaggerations since the most affected 
quantity is 7, and it will be affected considerably 
less at constant temperature than at constant 
pressure. 

Use of Table I enables us to estimate 
(V—V.)/V. as 0.026. Figure 2 of Guggenheim’s 
article! enables one to estimate (V—V.)/V. by 
comparing the ratio of density to critical density 
for neon and other gases. It appears that neon is 
about 4 to 5 percent expanded over what one 
might expect from the law of corresponding 
states. However, just the points near the melting 
points near the melting point for neon are a little 
out of line, so 4 percent expansion is probably the 
better estimate. Pitzer states that liquid neon has 
a density about 3 percent less than that of the 


“perfect liquid.’"” So it appears that the theo- 


17 See p. 588 of reference 3. 
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retical estimate of (V— V.)/V. cannot be far out 
of line. 

5C, is calculated as —1.42 cal. per mole. The 
specific heat of liquid neon is given" as 8.64 near 
the melting point, that of argon'* as 10.05 or'® 
10.5, that of krypton” as 10.6, and that of xenon” 
as 10.7. The classical values should all be the 
same by the law of corresponding states. It is, 
therefore, obvious that there is reasonably good 
agreement between the calculated and experi- 
mental values of 6C,. 


V. TWO-PHASE SYSTEMS 


We will now consider the effect of the quantum 
correction on the equilibrium between the solid 
and the liquid. We shall first consider the change 
in the temperature of melting at zero pressure. 
At zero pressure the Gibbs free energy F is equal 
to the Helmholtz free energy A. We can conse- 
quently use Eq. (3), substituting F for A and F, 
for A,.. F. really corresponds to A.o, since we 
should compare classical and quantum F’s at the 
same pressure. The difference between A. and 
A. is, as we have shown following Eq. (6), a 
second-order quantity. It is not entirely negli- 
gible for the solid, but is not important and may 
be neglected for our purposes. 

The condition for equilibrium between two 
phases may be written 


AF=AF,.—RTAn=0. (33) 


If there were no quantization equilibrium would 
occur at a temperature JT, with AF,.=0, while 
Eq. (33) holds at a temperature 7 where AF,=0 
+ (dAF,/dT),(T—T.) = —AS.(T—T.). Substi- 
tuting this into Eq. (33), we obtain 


T—T,= —RTAn/AS.. (34) 


Using the values of 7 for solid and liquid given in 
Sections III and IV, and taking for AS, the 
entropy of melting, 3.35, of argon, we obtain 
T—T,= —4.5°. The observed value of T—T, is 
about 24.55 —24.9= —0.35°. The discrepancy is 
very large, and quite surprising, in view of the 
apparently good results obtained with solid and 
liquid separately. 

18 A. Eucken and H. Hauck, Zeits. f. physik. Chemie 134, 
161 (1928). 

19 K, Clusius, A. Kruis, and F. Konnertz, Ann. d. Physik 
33, 642 (1938). 


* K. Clusius and L. Riccoboni, Zeits. f. physik. Chemie 
B38, 81 (1937), 
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We can analyze the situation further by con- 
sidering AE—AE, and AS—AS, for the fusion 
process. Since the melting temperature is so 
nearly a corresponding temperature with respect 
to the melting temperature of xenon, say, we can 
neglect any effect of change in temperature on 
these quantities. This means that we can calcu- 
late AE—AE, by applying Eq. (10) directly to 
both solid and liquid at 24.9°K. For the solid we 
apply the equations of Section III. For the liquid 
we apply the equations of Section II and the data 
of Section IV, and convert from argon to neon. 
We obtain E— E,=54.3 cal. mole for solid neon 
and E—£,.=18.3 for liquid neon. Thus, 


AE—AE,.= —36.0 cal. mole, calc. 


From the heat of fusion of krypton” we get for 
the corresponding heat of fusion of neon 83.6 cal. 
mole-!, while the actual heat of fusion is about 
80.1 cal. mole. Thus, 


AE—AE,= —3.5 cal. mole, obs. 


A similar discrepancy occurs for the entropy. 
We can write 


AS—AS, 
= (AE—AE,—AA+AA,)/T 
= (AE—AE,.)/T+RAn 
= —1.45+0.61 = —0.84 cal. mole deg.—. 


The actual value of AS—AS, is around —0.1 
cal. mole deg.—. 

The results with AE — AE, and AS—AS, appear 
to indicate that the quantum correction for the 
liquid is nearly as great as for the solid. This 
would seem very strange in itself, but it is still 
more puzzling in the light of the results of 
Section IV. The calculation for 6C, may indeed 
indicate that Eq. (13) gives something of an 
underestimate of 7; if it were increased by per- 
haps 40 percent a better value of C, would result. 
On the other hand, it must also be noted that the 
derivatives of » in Eq. (11) involve the quantity 
dE,/dV, which is somewhat indirectly calculated 
from the data. It does not appear likely, however, 
that this can cause appreciable error. We have” 


dE,,/dV = (@E/aV)r—(dE,/aV)r. 
(0E,/aV)r is certainly negative ; therefore dE,/dV 


*1 Zero-point energy is not included because argon is 
essentially classical. 
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should be slightly larger than (@E/0V)r which at 
zero pressure is equal to 7(dp/dT)y. Now the 
value of dE,/dV used is only slightly larger than 
T(0p/dT)y; it seems likely that dE,/dV should 
be increased rather than decreased, and this 
would actually increase | 5C;|, so that |n| would 
have to be increased less. The use of values at the 
normal melting temperatures and volumes gives 
a slight overestimate for |6C,|, which might 
mean that |7| should be increased by a little 
more than 40 percent, but this effect is not very 
great. 

Another quantity which occurs in the deriva- 
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tives of 7 which is not entirely directly determined 
experimentally is E,.. We have also used the 
calculated rather than the experimentally de- 
termined value of C,,.. However, it does not 
appear that any possible change in these quanti- 
ties could clear up the discrepancy. 

On the other hand, it appears likely that the 
discrepancy is too great to be accounted for on 
the basis of deviations from the law of corre- 
sponding states. 

We are, therefore, forced to leave the question 
of the quantum correction for neon in a some- 
what unsatisfactory state. 
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Bryce L, CRAWFORD, JR. AND WILLIAM Horwitz** 
School of Chemistry, University of Minnesota, Minneapolis, Minnesota 


(Received November 5, 1947) 


A redetermination of the depolarization factors of the Raman lines of nickel carbonyl shows 


the presence of only two polarized lines instead of three, as previously reported by other 
investigators. The data obtained not only eliminate the only evidence against the tetrahedral 


LTHOUGH the molecular structure of nickel 
carbonyl has been investigated several 
times, no clear-cut decision as to its configuration 
has been arrived at. The present status of the 
problem is well summarized by Herzberg,! who 
lists the structure of this molecule as tetrahedral, 
but adds in a footnote, ‘It must be mentioned 
that there is a serious discrepancy between this 
structure and the observed Raman spectrum in 
that the three strongest Raman lines are polar- 
ized whereas on the tetrahedral model there are 
only two totally symmetrical vibrations.” 
Raman spectra of nickel carbonyl were ob- 
tained by Dadieu and Schneider® and by Ander- 





_ *A portion of the thesis submitted by William Horwitz 
in partial fulfillment of the requirements for the degree of 
Doctor of Philosophy at the University of Minnesota, 
August 1947, 
_** Present address: U. S. Food and Drug Administra- 
tion, Minneapolis 1, Minnesota. 

1G. Herzberg, Infrared and Raman Spectra of Polyatomic 
Molecules (D. Van Nostrand Company, Inc., New York, 
1945), p. 356. 

* A. Dadieu and F. A. Schneider, Akad. Wiss. Wein 68, 
191 (1931). 


structure for this molecule but supply additional evidence in its favor. 





son,*? but the only conclusion that could be 
reached from the evidence of these investigations 
was that the intact carbon monoxide group was 
present in nickel carbony] rather than the organic 
carbonyl group. A more complete investigation of 
the Raman spectrum was performed by Duncan 
and Murray.‘ The data were interpreted on the 
basis of a planar square structure, a structure 
which these authors incorrectly ascribed to 
Pauling.® On the basis of the directed theory of 
valence, Pauling predicted that the complexes of 
bivalent nickel should have a planar configuration, 
but on the basis of the same theory complexes of 
neutral nickel should havea tetrahedral structure. 
Not only was the theoretical basis for the planar 
structure as used by Duncan and Murray incor- 
rect, but, as Wilson® pointed out shortly after- 
wards, the assignment was invalid since they had 


8 J. S. Anderson, Nature 130, 1002 (1932). 

4A. B. F. Duncan and J. W. Murray, J. Chem. Phys. 2, 
636 (1934). 

5 L. Pauling, J. Am. Chem. Soc. 53, 1367 (1931). 

6 E. B. Wilson, Jr., J. Chem. Phys. 3, 59 (1935). 
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TABLE I. The Raman data of nickel carbonyl.* 








Frequency Relative Depolarization Assign- 
cem~! intensity factor ment** 





791 (10) 100 0.864 0.07 (18) _ 


W5 


38144 (8) 25 ~0.01 (17) “ 
46143 (4) 10 0,830.11 (5) i 
598 (2) 3 <04 (1) Que 
2039+5 (9) 47 0.730.06 (18) e 


212147 (4) 3 <0.4 (1) wetws 








* Number of observations in parentheses. 

** From Crawford and Cross (see reference 8). The assignment of 
2121 cm™~ has been changed from w2+ws. 
violated the selection rules for the Raman effect. 
Their conclusion that nickel carbonyl was planar 
was, therefore, also invalid. Actually their ex- 
perimental data do not agree with either the 
tetrahedral or the square model, since either of 
these structures should have only two polarized 
Raman lines while Duncan and Murray reported 
the presence of three. 

Brockway and Cross,’ using electron-diffraction 
methods, found evidence clearly favoring a tetra- 
hedral structure, but a surprisingly small value 
was found for the nickel-to-carbon distance which 
may tend to alter the tetrahedral structure. 
Crawford and Cross® performed an infra-red ab- 
sorption study and a quantitative vibrational 
analysis on this molecule and found that a 
reasonable assignment of frequencies could not be 
made on the basis of the square model; a satis- 
factory assignment could be obtained for the 
tetrahedral model only if the Raman polarization 
data of Duncan and Murray were in error. The 
purpose of this investigation was to re-examine 
the state of polarization of the Raman lines of 
nickel carbonyl to determine whether Crawford 
and Cross were justified in rejecting the polariza- 
tion data of Duncan and Murray which represent 
the only evidence against the tetrahedral struc- 
ture for nickel carbonyl. 


APPARATUS AND MATERIALS 


Nickel carbonyl? of commercial grade was 
purified by several vacuum distillations. The 


7L. O. Brockway and P. C. Cross, J. Chem. Phys. 3, 
828 (1935). 

8 B. L. Crawford, Jr. and P. C. Cross, J. Chem. Phys. 6, 
525 (1938). 
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Raman spectra and depolarization factors were 
obtained with the apparatus described in a previ- 
ous publication.!° 

Wave-length nieasurements were obtained by 
comparison with a bracketing iron reference spec- 
trum. Ordinary nails were used for the arc in 
order to supplement the iron spectrum with 
manganese lines in the 2000-cm™ region. Meas- 
urements were made with a Bausch and Lomb 
spectrum measuring microscope or with a Leeds 
and Northrup recording microphotometer. The 
wave numbers (vacuum) were obtained by a 
linear interpolation over a limited region or by 
the use of the Hartmann formula over the entire 
spectrum. The frequencies and standard devia- 
tions are listed in Table I. The relatively low 
precision of these observations is the result of the 
use of a large number of different methods of 
measuring and calculating the spectra in an 
attempt to determine the most accurate and con- 
venient method of frequency determinations. 

The exposure time required to obtain a density 
of one unit for the strongest line (79 cm) of 
nickel carbonyl was 2.1 minutes as compared 
with 30 minutes for the 459-cm™ line of carbon 
tetrachloride and 13 minutes for the 992-cm™ line 
of benzene. This line is one of the most intense 
Raman lines which has been reported, but despite 
this fortunate circumstance the depolarization 
factor is difficult to measure because it lies on the 
“wing’’ of the exciting line as shown in the 
microphotometer tracing of Fig. 1. When the 
“minimum adjacent background’”’ method for 
correcting for background which was developed 
for the lines of carbon tetrachloride’ was applied 
to this line, the impossible values for p of greater 
than 1 were obtained. The background for this 
line was obtained by extrapolating the estimated 
contour of the exciting line and taking the inter- 
section of this curve with the position of the peak 
as the background density. The ‘‘minimum adja- 
cent background” was used for the correction 
with the 381-cm— and 461-cm=! lines while no 
background correction was required for the 2039- 
cm line. 


® We are indebted to Mr. O. B. J. Fraser of the Inter- 
national Nickel Company, Inc. for this material. 

10 B. L. Crawford, Jr. and W. Horwitz, J. Chem. Phys. 
15, 268 (1947). 
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RESULTS 


The only Raman polarization data on this 
compound available hitherto are the qualitative 
results of Duncan and Murray,‘ who report the 
following: the 463-cm™ line, depolarized; the 
382-cm™ line, polarized ; the 82- and 2043-cm7! 
lines, partly polarized. The quantitative depolar- 
ization factors, as determined’ in this investiga- 
tion, are presented in Table I together with the 
standard deviations and the number of observa- 
tions. There is no question that the 381-cm~ line 
is highly polarized. On 16 of 17 observations there 
was no trace of a parallel component; on only one 
tracing was there any suggestion of a parallel 
component, and this was only a slight jog in the 
tracing, of the same magnitude as several other 
jogs of no significance (i.e., arising from grain 
fluctuations, spots, dust, etc.). Even if this 
parallel component is considered significant, the 
estimated degree of depolarization of this line is 
about 0.01, or as close to zero as it is possible to 
demonstrate experimentally. 

These data also show that the 79-cm~ line is 
depolarized. The fact that this line actually lies 
on the “wing’’ of the Rayleigh line with sub- 
stantially different background densities on its 
two sides may have led Duncan and Murray‘ to 
the conclusion that it was partially polarized. 
Even under high microphotometer magnification 
we have been unable to confirm Duncan and 
Murray’s classification of this line as unsym- 
metrical. An unsuccessful attempt was made to 
resolve this line into the two lines postulated by 
Crawford and Cross,° using the long focal length 
camera of the Steinheil spectrograph which is 
capable of resolving doublets with a spacing of 
10 cm—. 

The 461-cm- and the 2039-cm— lines were 
found to be depolarized and partially polarized, 
respectively, asreported by Duncan and Murray. 
The latter is unsymmetrical. It too is a postulated 
doublet which could not be resolved with the long 
focal length camera. 

A single measurement of the degree of depolar- 
ization of the weak lines at 598 cm=! and 2121 
cm™~ showed them to be at least partially polar- 
ized. No attempt was made to study the still 
weaker lines reported by Duncan and Murray,‘ 
since the continuous background was approaching 
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Fic. 1. Microphotometer tracing of the polarized Raman 
spectrum of nickel carbonyl (upper—perpendicular; lower 
—parallel). The Raman lines are from different plates. 







substantial levels with the exposure time re- 
quired. 

The relative intensities found in this study, as 
shown in Table I, are in good agreement with 
those reported by Duncan and Murray.‘ 












DISCUSSION 






The present investigation reveals the presence 
of only two strong polarized lines, which agrees 
with the selection rules for either the tetrahedral 
or the square configuration for nickel carbonyl. 
The partially polarized line reported by Duncan 
and Murray at 82 cm™ is actually depolarized. 
As pointed out by Crawford and Cross,® a 
polarized line in this region is unaccountable by 
either model. The fact that the 381-cm™ line is 
completely polarized is in itself very strong evi- 
dence in favor of the tetrahedral structure since 
complete polarization is expected theoretically 
only for molecules of the cubic point groups. For 
all non-cubic point groups, which includes the 
planar square, the degree of depolarization of the 
totally symmetric Raman lines is intermediate 
between 0 and 6/7. 
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The vibrational analysis of Crawford and 
Cross!® which rejected the polarization data of 
Duncan and Murray for the 82-cm~ line, is 
shown, on the basis of the data presented here, to 
be completely justified. As a result of polarization 
measurements on the two weak overtone and 
combination lines at 598 cm and 2121 cm™ it 
has been necessary to change the assignment of 
the latter from we+ws to ws-++ws in order to allow 
for a polarized combination. The original assign- 
ment of 598 cm—! as 2ws remains unchanged, since 
this overtone may have a polarized component. 


L. CRAWFORD AND W. HORWITZ 


The tetrahedral structure for nickel carbonyl 
may be considered as established on the basis of 
the following evidence: 


1. The theory of directed valence predicts a tetrahedral 
structure for quadricovalent neutral nickel complexes. 

2. The electron-diffraction study of Brockway and Cross’ 
clearly favors the tetrahedral structure. 

3. As shown by Crawford and Cross,* a satisfactory 
vibrational analysis of the molecule can only be made on 
the basis of the tetrahedral model. 

4, The present investigation supplies the Raman polar- 
ization data which indicate a tetrahedral structure and 
which justify the previous vibrational analysis. 
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Origin of the Potential Barrier Hindering 
Rotation in Ethane and Related 
Substances 


Epwin N. LASSETTRE AND LAURENCE B. DEAN, JR. 
Department of Chemistry, The Ohio State University, Columbus, Ohio 
December 9, 1947 


TTEMPTS have been made by Aston, Isserow, Szasz, 
and Kennedy! to systematize the calculation of poten- 
tial barriers hindering rotation around single bonds by 
assuming a repulsion between hydrogens which was found 
empirically to vary inversely as the fifth power of the dis- 
tance. Since the energy of interaction of two quadrupoles 
varies inversely as the fifth power of the distance? these 
results have suggested to us that the barriers arise, not 
from the repulsion of hydrogen atoms, but from the electro- 
static interaction of the charges of one bond with the charges 
of another bond. We have investigated the potential barrier 
for ethane on the basis of this model with results which are 
in rather satisfactory agreement with experiment and 
which suggest that the ideas are of general applicability. 
The results are briefly described below. 

In considering the magnitude of interaction between 
bonds the following specific hypotheses are employed. (i) A 
single-bond charge distribution consists of a unit positive 
charge at each nucleus and a pair of electrons. (ii) The 
electron distribution of a single bond is cylindrically sym- 
metric with respect to the internuclear axis. (iii) The charge 
distributions of two bonds, which have no atom in com- 
mon, do not significantly overlap and, hence, exchange 
interactions are negligible. (iv) The electrostatic potential 
due to a single-bond charge distribution is sufficiently well 
characterized by its dipole and quadrupole moments, calcu- 
lated relative to an origin midway between the nuclei. 
These hypotheses are in qualitative accord with the ac- 
cepted ideas concerning the distribution of molecular elec- 
trons. If, however, this theory is to account satisfactorily 
for potential barriers, it is essential that the bond quadru- 
pole moments have reasonable magnitudes. 

In the first calculations, an approximate quadrupole mo- 
ment for the negative charge of the C—H bond was calcu- 
lated by assuming that the electron pair forming the bond 
is distributed in a localized molecular orbital which is a 
linear combination of atomic orbitals for carbon and hydro- 
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gen. The molecular orbital used is 


o =aNjr, exp — (1.6271)($+ § cosd)+bN2 exp(—r2), 


where 7; and 72 are the distances, in atomic units, of an 
electron from the carbon and hydrogen nuclei, respectively, 
and @ is the angle between r; and the line of centers. N; 
and N? are normalizing constants. The factor ($+ 3 cos@) is 
Pauling’s tetrahedral orbital,’ and the radial factor is ob- 
tained from Slater’s functions and screening constants.‘ 
The constants a and } can be determined from the normal- 
izing condition and the bond dipole moment (0.4 Debye). 
Since, however, the sign of the dipole moment is uncertain, 
we have preferred to vary a/b over a considerable range 
which includes both signs for the dipole moment. This 
procedure also indicates how sensitive the calculated 
barrier is to alterations in bond ionic character since a/b 
can be considered a measure of ionic character. The results 
are given in Table I. In Table I the second column contains 


TABLE I. 








P 
pX10*18 atomic D-D D-Q Q-Q Barrier V (stag.) 
k ke kcal. 





a/b e.s.u. units cal. kcal. al. kcal. 

1.0 —2.30 3.16 1.58 —3.33 5.63 3.88 147.1 
1.5 —1.51 3.22 0.68 —2.22 5.84 4.30 96.4 
2.0 —0.98 3.22 0.29 —1.45 5.84 4.68 70.3 
4.0 0.03 3.13 0.00 0.05 5.52 5.57 38.3 
5.75 0.40 3.07 0.05 0.56 5.30 


5.91 33.1 








the calculated bond dipole moment and the third column 
the bond quadrupole moment. The “‘quadrupole moment,” 
P, is defined for a cylindrically symmetric distribution of 
charge density p by the formula 


P= f (r°—32%)pdr, 


where x is the coordinate perpendicular to the axis of sym- 
metry. This is the only parameter required to characterize 
the quadrupole field of such a distribution. The interaction 
energy of the three C—H bonds of one methyl group with 
the three of the other methyl group was calculated for both 
the staggered (D3z) and the eclipsed (Ds) configurations. 
The difference between these is given in column seven and 
the separate dipole-dipole, dipole-quadrupole, and quadru- 
pole-quadrupole energies are also given in columns four, 
five, and six, respectively. The interactions between C —C 
and C—H bonds need not be calculated since they remain 
unchanged as a methyl group is rotated around its sym- 
metry axis. 

The results may be briefly summarized as follows: 

1. The staggered configuration has lower energy for 
every value of a/b. 

2. The calculated potential barrier is larger than that ob- 
served (2750 cal.5), but is of the proper order of magnitude. 

3. For those values of a/b which give the correct magni- 
tude to the bond dipole moment (0.4D) the quadrupole- 
quadrupole interaction predominates, and hence an inverse 
fifth-power dependence on distance between bonds is 
confirmed. 

4. The barrier is not very sensitive to bond ionic 
character. 

5. In performing the calculations it was observed that, 
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although certain quadrupoles attract and others repel, the 
net interaction of methyl groups in ethane is one of repul- 
sion both for the staggered and the eclipsed configurations. 

6. The calculated potential barriers are small, being the 
difference between interactions in two configurations, but 
the total interaction energy is large. (See column eight.) 

Since a considerable weight of indirect evidence now 
exists in favor of the staggered configuration for ethane,’ ® 
the above are seen to be in as good agreement with experi- 
ment as could be expected in view of the rough approxima- 
tion to the negative charge distribution which has been 
employed. Accurate agreement between the calculated and 
the experimental barrier (2750 cal.) is obtained with 
quadrupole moment 2.03 (atomic units) if 4=0.4D and 
2.35 if w= —0.4D. The calculated interaction energies are 
given in Table II. It is interesting to note that both of these 
quadrupole moments are considerably smaller than those 
of Table I. This implies that the actual bonding distribu- 
tions overlap even less than is indicated by the rough mo- 
lecular orbital used here, hence is consistent with our initial 
hypothesis. It seems, therefore, that the potential barrier 
for ethane can be accounted for by means of a quadrupole 
moment of reasonable magnitude, and the staggered con- 
figuration is favored. Thus far the mutual polarization of 
bonding distributions has been ignored. It is obvious that 
such polarization occurs but it is difficult to see how such 
interactions can stabilize the staggered configuration. 

It is of interest that quadrupole moments obtained from 
observed potential barriers (which depend only on energy 
differences) may be used to calculate the force exerted by 
one bond on another which contains no common atom. The 
large energy of interaction of methyl groups in ethane 
(see Table II, column 7) indicates that such forces are 
appreciable. In this way a more quantitative treatment of 
the dependence of interatomic distances, hydrogenation 
energies, dissociation energies, etc., on bond environment, 
along the lines proposed by Conn, Kistiakowsky, and 
Smith,* becomes possible. 

Although we have found the same dependence on dis- 
tance for interacting bonds of small dipole moment as did 
Aston and collaborators, the present model also provides 
for dependence on the quadrupole moments of two inter- 
acting bonds as well as their orientation relative to each 
other and to the line joining them. Moreover, in a complex 
substance, such as tetramethyl methane, substantial inter- 
actions between C—H and C—C bonds are expected, and 
hence the similar dependence on distance may be fortuitous. 

The conception of interaction of bonding distributions 
as constituting the source of potential barriers is essentially 
similar to that proposed by Kistiakowsky, Lacher, and 
Ransom’ several years ago. In certain respects the present 
treatment may be considered as a quantitative extension 
of their ideas. 











TABLE II. 
P 
uw X10*18) =aatomic D-D D-O Q-0 Barrier V (stag.) 
e.8.u. units kcal. kcal. kcal. kcal. kcal. 
0.40 2.03 0.05 0.37 2.33 2.75 14.2 
—0.40 2.35 0.05 —0.43 3.13 2.75 29.0 
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bonding distributions is in progress along two lines: (a) A 
more accurate calculation of the electrostatic potential 
than is provided by dipole and quadrupole terms, (b) the 
empirical calculation of quadrupole moments from the ob- 
served potential barriers of compounds other than ethane, 
and the use of such moments in the prediction of other 
properties in order to check the accuracy of these values. 
The calculations will be fully reported when this investiga- 
tion is complete. 


1 J. G. Aston, S. Isserow, G. J. Szasz, and R. M. Kennedy, J. Chem. 
Phys. 12, 336 (1944). 

2H. Marganau, Rev. Mod. Phys. 11, 1 (1939). 

3L. Pauling, Nature of the Chemical Bond (Cornell University Press, 
Ithaca, New York, 1940), p. 86. 

4 J. C. Slater, Phys. Rev. 36, 57 (1930). . 

5G. B. Kistiakowsky, J. R. Lacher, and F. Stitt, J. Chem. Phys. 7, 
289 (1939). See also J. D. Kemp and K. S. Pitzer, J. Chem. Phys. 4, 
749 (1936). 

6 J. B. Conn, G. B. Kistiakowsky, and E. A. Smith, J. Am. Chem. 
Soc. 61, 1868 (1939). 

7G. B. Kistiakowsky, J. R. Lacher, and W. W. Ransom, J. Chem. 
Phys. 6, 900 (1938). 





Physical Basis of a Concept about the 
Mechanism of Viscosity of Liquids 


R. PARSHAD 


Research Laboratory of Electronics, Massachusetts Institute 
of Technology, Cambridge, Massachusetts 


December 15, 1947 


THEORY of viscosity of liquids should explain the 

following: (1) decrease of viscosity with temperature; 

(2a) increase of viscosity with pressure; (2b) the rate of in- 

crease of viscosity increasing with increase of pressure; 

(3) increased viscosity of associated liquids over that of 
normal liquids. 

The current theories of viscosity generally relate cohesion 
with viscosity. Cohesion may be taken as a measure of 
total attraction energy of molecules and is indicated in 
magnitude by latent heat of vaporization or inverse func- 
tion of compressibility. While the relation between cohesion 
and viscosity is generally true, there are significant de- 
partures from this correspondence. These originate in the 
difference between the mechanisms of viscosity and co- 
hesion. Viscosity arises because of the resistance to the dis- 
placement of a molecule at right angles to the intermolecu- 
lar line (taking an ideal case of two molecules), and is, in 
general, due to the resistance offered to change of con- 
figuration of the molecules, or the resistance against chang- 
ing of places of individual molecules in the same configura- 
tion, while the intermolecular distances remain unchanged. 
Cohesion arises due to resistance offered to a direct change 
of intermolecular distances, the configurations remaining 
the same. These differences have been discussed more fully 
elsewhere.! 

We assume that viscosity, or the resistance offered to 
change of position of molecules, is due to the existence 
of potential barriers around the molecules in the liquid. 
The steeper these potential barriers, the more the viscosity 
will be. This is why liquids of anisotropic molecules have, 
in general, higher viscosity than liquids of more sym- 
metrical molecules. Now we postulate that these potential 


A further investigation of the electrostatic interaction of 
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barriers are caused by the interaction of the two opposing 
intermolecular forces of attraction and repulsion acting 
individually. For high viscosity not only high attracting 
forces are needed but also high repulsive forces. By the 
interaction of these two opposite forces the molecules tend 
to get fixed up in space so to speak; this fact means high 
viscosity. The departure from complete correspondence 
between cohesion and viscosity arises because cohesion is 
a measure of difference of attractive and repulsive energies 
while viscosity depends individually upon both. Our con- 
cept at once explains the following pertinent facts. 

(1) Increase of viscosity with pressure. With decrease of 
volume under pressure, both the attractive and repulsive 
forces increase, leading to an increase of viscosity. More- 
over, as these forces are inverse functions of distance, their 
rate of increase will increase with the decrease of the inter- 
molecular distance. Further, since the repulsive forces are 
of much shorter range they will, after a stage, increase 
much faster than the attractive forces. This leads to a rate 
of increase of viscosity which increases fast with pressure. 
On the other hand, the rate of fall of compressibility with 
pressure will be rapid in the beginning and slower after- 
wards as it is related to the difference of attractive and 
repulsive energies, which increases in the beginning but 
decreases at still higher pressures. 

(2) Viscosity and association. It has been shown? that 
molecules of associated liquids do not only have great 
attractive forces but also have significant repulsive forces. 
This may be physically explained as being caused by the 
hydrogen bonds, in which the distance between two oxygen 
atoms of interacting molecules (OQ—H——O) decreases so 
much as to produce repulsive forces as well. This tendency 
is further increased due to the semipermanent nature of the 
hydrogen bond links, so that, due to the hindered relation 
of the molecules concerned, the intermolecular forces 
become more of a first order than second order. 

(3) Behavior of water. Remembering that compressibility 
decreases with breaking of associations,’ and that the as- 
sociations break under pressure and with increase of tem- 
perature, the odd behavior of water in its change of vis- 
cosity and compressibility with pressure and temperature 
will be explained. 

(4) Decrease of viscosity with temperature. This is the 
easiest thing to explain and also follows from the current 
theories. With the increase of temperature at constant 
volume, both the attractive and repulsive energies will 
increase because of the smoothing out of intermolecular 
distribution functions. Even if that is not so or if the effect 
is negligible, the increase of kT will increase the probability 
of jump of the potential barriers, leading to decreased 
viscosity. 

At constant pressure, the intermolecular distances in- 
crease, leading to a more rapid decrease of both kinds of 
forces and so leading to a faster decrease of viscosity. 

Further details of the theory and its quantitative aspects 
will be taken up in a later paper. 

1For example, Andrade’s Theory, Phil. Mag. 17, 494, 698 (1934); 
Eyring’s Theory, Glasstone, Laidler, and Eyring, Theory of Rate 
Processes (McGraw-Hill Book Company, Inc., New York, 1941). 


* Parshad, Doctorate Thesis, Punjab University, India, 1946. 
* Parshad, Ind. J. Phys. 19, 47 (1945). , 
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Small-Angle X-Ray Scattering from Metal 
Deposits Made by Evaporation 
BENJAMIN CARROLL 
The Newark Colleges, Rutgers University, Newark, New Jersey 
AND 
I. FANKUCHEN 


The Polytechnic Institute of Brooklyn, Brooklyn, New York 
December 2, 1947 


HE formation of dull or black films during the evapo- 
ration of a metal has been the subject of considerable 
conjecture. Numerous diffraction studies of “evaporated” 
metal films have been reported,'** but the small-angle 
scattering properties of these deposits have been over- 
looked. We have examined aluminum and copper deposits 
produced at 10-2 and 10-* mm by means of Debye- 
Scherrer and small-angle diagrams, and have found that 
the black deposits exhibited marked small-angle scattering 
of x-rays (Fig. 1) and bright deposits did not; however, the 
Debye-Scherrer diagrams were substantially alike regard- 
less of the wide variations in the reflectivity of the metal 
surface toward white light. 

Aluminum was evaporated from a tungsten wire and 
deposited on glass surfaces at various distances from the 
source. In the case of copper, a tantalum heater was used 
permitting the evaporation of relatively large quantities 
of metal. To minimize the possibility of compound forma- 


| I 


i F 


Fic. 1. Small-angle x-ray scattering of bright and dark metallic 
deposits. A. Bright aluminum foil made at 10-4 mm. B. ey. copper 
foil made at 10-4 mm. C. Aluminum black made at 10-? mm. D. Copper 
black made at 10-2? mm. E. Bright aluminum foil made at 10-? mm 
(when condensation is close to source). F. Bright copper foil made at 
10-2 mm (when condensation is close to source). 
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tion for the metal deposits formed at 10-? mm, the pressure 
in the evaporation chamber was regulated by allowing 
helium to leak into the system. The deposits were scra 

from the glass base and examined with filtered copper 
radiation. The distance between sample and photographic 
plate was 50 cm for the small-angle work and the method 
used to evaluate particle size was the same as described 
by Jellinek and Fankuchen.‘ The identical specimens were 


-then used in the Debye-Scherrer camera for the purpose 


of observing line-broadening effects. 

The Debye-Scherrer patterns of all samples showed no 
particular orientation effects and aside from general darken- 
ing of the background of the photograms for the black 
deposits, line broadening, due to small particle size was not 
apparent (see Fig. 2). The average particle size for alu- 
minum and copper black as calculated from the small-angle 
scattering was 110A and 90A, respectively; the bright 
deposits formed close to the heater during the same evapo- 
ration experiment (pressure 10-? mm) gave negligible 
scattering at small angles, indicating that the crystals were 
about 10 cm. 

Jellinek and Fankuchen® showed that for fairly homo- 
geneous specimens of y-alumina, crystallite sizes calculated 
from small-angle scattering and line-breath measurements 
of Debye-Scherrer patterns were in good agreement. If 
gross inhomogenities exist in a specimen as in the case 
of the aluminum and copper black, small-angle scattering 
will indicate the presence of the fine particles or crystallites 
only, whereas the Debye-Scherrer diagrams will reveal the 
existence of the coarser crystallites. 

The formation of black deposits by evaporation can be 
considered a general method of preparing colloidal metal 
particles, the colloidal condition apparently being caused 
by the arrival of clusters of atoms at the surface. Reorgani- 
zation and crystal growth of the film depends upon the 
surface mobility of the condensed atoms. Since the activa- 
tion energy for migration over the surface of a solid is 
greater for an atom in a cluster than for a single atom, on 
the surface of the sub-phase, we have an explanation for 
the invariable occurrence of black deposits when the mean 








Fic. 2, A. Aluminum black made at 10-2 mm. 
B, Aluminum bright foil made at 10~* mm. 
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free path inside the vacuum chamber was much smaller 
than the distance between the source and the condensing 
surface. The formation of bright deposits and relatively 
coarse crystallites close to the source when the pressure 
was 10-? mm were probably not only due to the compara- 
bility of the mean free path and the distance from the 
source to the condensing surface, but also to the somewhat 
higher temperature that prevailed near the source. 

1 Andrade, Trans. Faraday Soc. 31, 1157 (1935). 

2 Barrett, Structure of Metals (McGraw-Hill Book Company, Inc., 
New York, 1943), p. 440. 

3 Milligan and Focke, J. Phys. Chem. 45, 107 (1941). 

4W. C. Caldwell, J. App. Phys. 12, 779 (1941). 


5 Jellinek and Fankuchen, Ind. Eng. Chem. 37, 158 (1945); also see 
Jellinek, Solomon, and Fankuchen, ibid., 18, 172 (1946). 





A New Concept about the Mechanism of 
Passivity of Iron in Nitric Acid 


R. PARSHAD* AND L. C. VERMAN** 


Physical Laboratory, Council of Scientific and Industrial Research, 
Delhi, India 


December 19, 1947 


ANY theories! have been advanced from time to time 
to explain the passivity of iron in concentrated nitric 
acid. In the following we offer a new concept. 

In 1944 electrochemical evidence was obtained? showing 
the existence of a hydrogen or hydride film on aluminum, 
spaced between the well-known oxide film and the metal 
itself. In our view, the tenacity and strength of the oxide 
film is due to this hydrogen film. 

We were led to believe that the passivity of iron in con- 
centrated nitric acid might also be due to this intermediate 
hydrogen film, with the difference that, whereas in the case 
of aluminum the film could be produced in the atmosphere, 
on iron it would come from the acid itself, possibly as a re- 
sult of an initial metal-acid interaction. This hydrogen film 
adsorbed on iron would, in our view, attract an oxygen or 
oxygen-rich film (helped in this process by the attraction 
of iron itself for oxygen), and the latter oxygen-rich film 
would stop further oxidation. 

Experimentally, when a piece of iron is inserted in nitric 
acid, a potential due to this hydrogen film should be ex- 
pected before the known noble potential due to the oxygen- 
rich film is recorded. A platinum electrode was dipped in 
concentrated nitric acid and some time later an iron elec- 
trode was at once inserted in the acid. The potentials pro- 
duced right from the making of the cell were amplified by 
a d.c. amplifier and measured on a cathode-ray oscillo- 
graph. In acid of concentration 68.8 percent by weight at 
32.6°C, the potential of the cell rose to 1.25 volts almost 
instantaneously, and then in 0.05 second fell to a steady 
potential of about 0.2 volt..During its decrease, the poten- 
tial passed through a “‘plateau” where its rate of fall was 
the lowest. As concentration of acid decreased, the initial 
rise of potential occurred in two steps: first an almost in- 
stantaneous rise as before and then a continued slower rise 
through about 0.1 volt. The rate of fall of potential after it 
had reached the maximum also decreased. For example, in 
acid 45.6 percent by weight, at 32.6°C, the time taken to 
reach the steady low potential was one second. 

We have reasons to believe that the initial potential is 
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' due to the formation of the hydrogen film on iron and the 


continued rise of potential at lower concentrations indi- 
cates the continued but slower formation of the film. 

The following observation is a significant support to our 
theory. An iron electrode, which had been kept unused for 
some years, took about 10 seconds of chemical action before 
gaining passivity in acid, 41.4 percent by-weight, while the 
same iron, with its surface renewed by having it abraded 
with emery-cloth, took about 25 seconds. According to 
previous theories of passivity, this fact is due to the forma- 
tion of a thin invisible film of iron oxide or of adsorbed 
oxygen on the old surface of iron. The potential measure- 
ment of the cell showed, however, a higher potential when 
the electrode of old iron surface was used than that ob- 
tained for the use of the fresh surface. In the older theories, 
the relative magnitudes of the potentials should have been 
in the opposite order. In our view, the greater potential of 
the cell for the older surface of iron was due to its adsorp- 
tion of some hydrogen film from the atmosphere, which 
helped the iron to gain passivity in the dilute acid. 

Many other phenomena associated with passivity are 
easily explained by the proposed mechanism. For example, 
the corroding action of halides is due to the affinity of the 
halogens for hydrogen, resulting in destruction of the film 
and so of passivity. 

Details of the experimental data and application of the 
concept to behavior of stainless steels will be published 
elsewhere. 

* Now at Research Laboratory of Electronics, Massachusetts Insti- 
tute of Technology, Cambridge, Massachusettes. 
** Now at Indian Standards Institution, New Delhi, India. 

1 For a review see Evans, Metallic Corrosion, Passivity, and Protection 

(Edward Arnold and Company, London, 1946); Uhlig and Wulff, Am. 


Inst. Min. and Met. Eng. 135, 494 (1939). 
2 R. Parshad, Nature 154, i78 (1944). 





The Infra-Red Emission and Absorption of the 
Carbon Monoxide-Oxygen Flame* 


S. SILVERMAN, G. A. HORNBECK, AND R. C, HERMAN 


Applied Physics Laboratory, The Johns Hopkins University, 
Silver Spring, Maryland 


December 5, 1947 


ECENTLY developed combustion engines, such as 
rockets and jet-engines of various types, have led to 
a renewed interest in the infra-red emission of flames, ex- 
plosions, and combustion reactions in general. Numerous 
investigators, beginning with Julius and Paschen in the 
1890’s have studied the infra-red emission of different 
flames and hot gases.! The mission bands corresponding to 
the well-known absorption bands of the HO and CO; 
molecules were observed in the Bunsen flame in the region 
from about 1 to 124. Until recently? the only observation 
of radiation in the longer wave-length region, from about 
10 to 204, was made by Rubens and Aschkinass? in 1898, 
who found emission throughout that region from both the 
Bunsen flame and hot COs. Subsequently, investigators of 
the CO-oxygen flame‘ reported the strong emission bands 
of COe at 2.8 and 4.44 but failed to observe the fundamental 
at 14.95y. 
The CO-oxygen flame has been studied to determine, if 
possible, the vibrational excitation of the CO2 molecules 


formed in this reaction. The flame employed in these in- 
vestigations was obtained by burning a mixture of CO 
and O2 in an ordinary torch. The CO was commercial grade 
and contained as its chief impurities about 1.5 percent each 
of Hz and Ne, but no H:O emission was observed in the 
infra-red. The iron carbonyl was removed from the CO in 
the usual manner. The instrument used was a Perkin- 
Elmer infra-red spectrometer with NaCl or KBr optics. It 
was flushed continuously with dry Nz to minimize atmos- 
pheric absorption. The receiver was a bolometer using 
chopped radiation at 10 cycles/sec. with tuned amplifica- 
tion for drift-free recording purposes. 

The emission spectrum of the CO-oxygen flame shows 
in addition to the peaks previously observed a number of 
rather weak emission bands between about 10 and 20uy. 
The two strongest bands lie at 667 and 721 cm™, and agree 
to within about 1 cm with the transitions (0110-000) 
and (10°0-—01'0) of the CO» molecule in the ground elec- 
tronic state. There are other bands in this region as shown 
in Fig. 1. Relative to the 4.44 band the intensities of the 
bands are quite low and it is not surprising that some earlier 
workers were unable to detect them. We have also repeated 
the study of the emission of hot CO: gas streaming from a 
furnace at about 1400°K. The emission spectrum is com- 
pared to that of the flame in Fig. 1. 

To determine the absorption spectrum of the flame, 
radiation from a Globar was modulated and passed through 
it. In this way the emission from the flame itself is elimi- 
nated except for random intensity variations. Absorption 
bands are found at 2.8, 4.4, 14.1, and 14.95 of which only 
the two fundamentals are strong. 

Further work on the CO-oxygen and other flames is 
being conducted and a more complete report will be made 
in the future. We wish to express our thanks to Dr. E. K. 
Plyler of the National Bureau of Standards who has shown 
us the results of similar studies on the Bunsen flame? and 
who reports finding the long wave-length bands of COs. 
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Fic. 1. Curve A: Emission from CO-oxygen flame in 15, region. Slit 
width =0.3 mm. Curve B: Thermal emission from hot COs, uncorrected 
for background radiation from furnace, which is shown separately in 
curve C. Slit width =0.5 mm. All the curves were obtained with a KBr 
prism and are unsmoothed for noise variations; the maximum noise is 
about 2 percent of the peak deflection. 
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We are also indebted to Mr. B. W. Bullock for valuable 
help in experimental matters. In conclusion we express our 
appreciation to Professor D. M. Dennison and Dr. E. O. 
Salant for their interest and many valuable discussions. 


* The work described in this letter has been supported by the Bureau 
of Ordnance, U. S. Navy, under Contract NOrd-7386. 

1For an extensive bibliography, see Jean Lecomte, Le Spectre 
Infrarouge (Les Presses Universitaries de France, Paris, 1928), and 
A. G. Gaydon, Spectroscopy and Combustion Theory (Chapman Hall 
Limited, London, 1942). 

2 E. K. Plyler, J. Opt. Soc. Am. 37, 984A (1947). 

3H. Rubens and E. Aschkinass, Wied Ann. 64, 584 (1898). 

4W. E. Garner and C. H. Johnson, Phil. Mag. 3, 97 (1927); C. R. 
Bailey and A. H. Lih, Trans. Faraday Soc. 25, 29 (1929). 





A Cylindrical Model for the Small 
Soap Micelle* 


WILLIAM D. HARKINS 


George Herbert Jones Chemical Laboratory, University of Chicago, 
Chicago 37, Illinois 


December 16, 1947 


HE model presented here (Fig. 1) is based on exten- 

sive papers on x-ray relations by Mattoon, Stearns, 

and Harkins,! and on a great deal of other work by Corrin, 

Oppenheimer, Mittelmann, Herzfeld, Roginsky, and Lind. 

A quotation from a Cambridge University Symposium? 
(1947) is as follows: : 


There are different schools of thought on these matters, 
for example McBain postulates two types of micelle, 
‘ionic’ and ‘‘neutral,’’ Hartley only one, the ‘‘spherical’’ 
ionic micelle. The latter view is now generally accepted, 


and will be followed here. 


McBain’s views are expressed more accurately in his 
1942 review.* His ionic micelle is spherical (adopted by 
Hartley) and his neutral micelle identical with Stauff’s 
“Grossmizellen” or lamellar micelles. 

Soap micelles are: 






























































Fic. 1. Highly idealized cross section of a soap micelle. A. Without 
solubilized oil. B. With solubilized oil. The solubilization of an oil is 
found to increase the number of soap molecules in the micelle. Although 
a cylinder seems to represent most of the properties of the micelle as 
now known, it is imperfect in that it represents too large an interface 
at the side between hydrocarbon groups and water. Thus it seems not 
improbable that a model which lies between a cylinder and a cylindrical 
type of spheroid may be found to more perfectly represent the energy 
relationships. This would amount to a distortion of the side of the 
cylinder by means of which polar groups cover to some extent the 
non-polar hydrocarbon chains. 


2. Not spherical and exhibit a thickness, as measured 
by a new x-ray M-band,™ equal to the double length of 
the molecule. 

3. This thickness is increased by the solubilization of a 
non-polar oil, but not by alcohols or amines. 

4. When a salt is added the critical micelle concentration 
is decreased, the micelle size increased by increase in the 
diameter of the cylinder, and the principle of ionic strength 
and Debye-Hiickel theory becomes invalid (Corrin and 
Harkins‘). 

5. The number of soap molecules per micelle is increased 
by solubilization of an oil and by long chain alcohols and 
amines. 

By diffusion Hartley® found the cetyl pyridinium chloride 
micelle to be small and to contain about 70 molecules. Our 
x-ray work on different detergents indicates sizes of this 
order, but both smaller and considerably larger. 

Our earlier x-ray work appeared to favor the lamellar 
model of McBain and German investigators,*»* which 
consists of 5 or 6 large double layers of soap with inter- 
layers of water of definite thickness (7) as calculated from 
our data for potassium myristate. 

After a long and intensive study, in an endeavor to find 
a mechanism which would control the interlayer or 
“‘bound”’ water in this way, the writer concluded that there 
is no such mechanism, and that these water layers are 
non-existent, but instead each double layer of molecules is 
a micelle and moves freely as a large ‘“‘molecule’”’ through 
the solution, until a second critical concentration (ca. 7 
percent for potassium myristate) is reached. Above this 
concentration each double layer is a micelle, but the mi- 
celles now come closely enough together to restrict their 
motion (libration or rotation) so that above this concen- 
tration the x-ray J-band appears. This is related to the 
distance between the micelles. 

Calculations of the thickness (7) of the oil layer solu- 
bilized by the micelle, as made by ‘the writer, proved 
incorrect the conception of this process as related to the 
lamellar model. Those who developed this model assumed 
(McBain’s review, p. 124) that the increment (Ad;) in the 
long spacing due to the oil gives the thickness of the oil 
layer. The mean thickness (7) should equal the volume of 
oil in the micelles divided by the area of the large double 
soap layers, which is known, since each pair of molecules 
occupies 28 sq.A. In calculations for n-heptane, triptane, 
and ethyl benzene the mean thickness (r) was found at all 








Water layer thickness 
Percent soap rinA 


34 

31.7 
26.3 
22.5 
19.6 
17.2 
15.2 











concentrations of the oils to be only 0.4 as large as Ad;. For 
example, for 3.5 percent m-heptane in 25 percent potassium 
laurate r=5.5A but Ad;=13.6A or 2.5 times larger than 
it should be. It was concluded that Ad; is not the distance 
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between layers in a lamellar micelle, but the change in the 
intermicellar distance (d;). 

With respect to solubilization, Corrin found that Hart- 
ley’s spherical model gave much better agreement than this. 

The cylindrical model is superior in that its height is 
measured by the M- or micelle thickness band, and it has 
a diameter which may be varied and can be calculated 
approximately at high concentrations from the J-band, or 
at low concentrations from light scattering (Debye). 

In any model the energy relations make it essential to 
make the interface hydrocarbon chains-water as small as 
possible, which means that the cylinder may be somewhat 
distorted at the side to give coverage of the hydrocarbon 
chains by a few polar-ionic groups of the soap. However, 
the model is essentially a cylinder with these groups at 
the ends. 


* This investigation was carried out under the sponsorship of the 
Reconstruction Finance Corporation, Office of Rubber Reserve, in 
connection with the government’s synthetic rubber program. 

la R, W. Mattoon, R. S. Stearns, and W. D. Harkins, submitted for 
publication (September, 1947); > ibid. J. Chem. Phys. 15, 209 (1947); 
eW. D. Harkins, R. W. Mattoon, and R. Mittelmann, bid. 15, 763 
(1947). 

2A, E. Alexander, Colloid Science. A Symposium (Chemical Publish- 
ing Company, Inc., Brooklyn, 1947), p. 43. McBain now considers the 
“neutral’’ micelle to be the somewhat ionized “‘large lamellar micelle.” 

3a J. W. McBain, Advances in Colloid Science (Interscience Publishers, 
Inc., New York, 1942), edited by E. O. Kraemer, p. 124, for extensive 
references to work on the “‘large lamellar micelle’ see: > W. D. Harkins, 
R. W. Mattoon, and M. L. Corrin, J. Am. Chem. Soc. 68, 220 (1946), 
and ¢ J. Colloid Sci. 1, 105 (1946). 
as fl L. Corrin and W. D. Harkins, J. Am. Chem. Soc. 68, 679, 683 

47). 
as G. S. Hartley and D. M. Runnicles, Proc. Roy. Soc. A168, 420 

938). 





The Solubility of Sodium Soap in 
Organic Solvents 


A. BonpbI 
Shell Development Company, Emeryville, California 
November 18, 1947 


HE solubility of sodium palmitate and sodium stearate 

in organic solvents has been determined by several 
members of the Stanford Group.'~* These data seem rather 
difficult to correlate with each other or with other solubility 
data. The latter is usually most fruitful if the comparison 
is made on the basis of Hildebrand’s relation: 


AH;/1 1 


where AHy the heat of fusion (or transition) and Ty the 
melting (or transition) temperature.*® 

In view of the fact that the sodium soaps exhibit not 
only a melting point but a series of transition points in 
the solid state, the ideal solubility curve changes its slope 
at as many points. The transition points and the transition 
heats have been determined by Vold.* It is therefore 
possible to construct the ideal solubility curve of these 
soaps. Such a curve for sodium palmitate is shown— 
marked ‘‘ideal’’—on Fig. 1. The other curves on the same 
figure are the so-called 7; curves,’ giving the temperature 
at which the isotropic melt changes into an anisotropic 
liquid crystal. The proximity of several curves to the 
Raoult’s law curve is rather surprising, since one usually 
associates “ideality” with simpler physical behavior. It 





should be noted that the slope of all solutions in hydroxylic 
compounds deviates in the low temperature (<180°C) 
region but little from the slope of the ideal curve, while 
the slope of the other solutions, particularly in hydro- 
carbons, is very much steeper. This suggests that the 
excess free energy of mixing with the alcohols consists of a 
negligible heat of mixing and a strongly negative entropy 
of mixing. The solution of soap in hydrocarbons requires 
on the other hand both rather large (negative) enthalpy 
and entropy changes. One would hardly have expected 
that the solubility of stearic acid in similar solvents 
deviates more from ideality than that of soap, as is evident 
from the data on Fig. 2. This applies particularly to 
alcoholic solutions. 

In the high temperature region all solutions behave 
quite similarly in exhibiting solid solution phenomena, as 
is also apparent in visual observation of these systems. 
The effect of chemical structure of the solvent on the shape 
of the solubility curve is sometimes quite marked, as e.g., 
in the case of the cresols. The position of p-cresol solutions 
is given on the graph. The solubility curve of soap in 
o-cresol, on the other hand, is practically identical with that 
of cetyl alcohol and has, therefore, not been drawn in. In 
other cases there is little specificity as is exemplified by 
the solubility in cyclohexane and in low mole weight 
aromatic hydrocarbons, or cetane and white medicinal oil. 

The humps in the solubility curves which are observed 
in soap solutions in cyclic (low molecular weight) hydro- 
carbons, in water, and in glycerol are somewhat reminiscent 
of the shape of solubility curves of systems forming two 
immiscible liquid phases in a certain temperature range. 
While presentation of the soap solubility data in the 
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Fic. 1. Solubility of sodium palmitate in various organic solvents 


(from data of references 1, 2, 3) compared with ideal solubility, 
curve (1). 


(3) Diethylene glycol 
(4) Ethyl alcohol 
(S) Cetyl alcohol 


1) Ideal solubility 
2) p-cresol 


(6) n-butylamine 

(7) n-cetane 

(8) Nujol 

(9) Cyclohexane, sodium stearate in 
(10) Ethylbenzene, sodium stearate in 
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Fic. 2. Solubility of stearic acid in various organic solvents (from data 
of reference 8) compared with ideal solubility, curve (1). 
(1) Ideal solubility 
(2) Cyclohexane 
(3) n-butyl alcohol 
(4) Methyl alcohol 








manner shown in Fig. 1 simplifies the discussion of these 
complicated systems, it leaves many phenomena still 
unexplained. 

The writer was unable to find any solubility data of 
simpler compounds undergoing phase transition in the 
solid state, such as camphor, C2Cl¢, and the like, to com- 
pare with the ssqunay more thoroughly studied soap 
systems. 

The writer wishes to acknowledge discussions with 
Professor J. W. McBain and Dr. O. Redlich, which led to 
this work. 


1R. D. Vold, C. W. Leggett, Jr., and J. W. McBain, J. Phys. Chem, 
44, 1058 (1940). 

2C. W. Leggett, Jr., R. D. Vold, and J. W. McBain, J. Phys. Chem. 
46, 429 (1942). 

3G. Ly Smith, Ph.D. Thesis, Stanford, 1946. 

4j. Hildebrand, Solubility of Non-Electrolytes (Reinhold Pub- 
lishing gcc New York, 1936). 

5 A. Bondi, J. Phys. Chem. 51, I (1947). 

®R. D. Vold, J. 7. Chem. Soc. 63, 2915 (1941). 

7J. W. McBain, L. H. Lazarus, and A. V. Pitter, Zéits. physik. 
ey A147, 87 (1930). 

8 C. W. Hoerr and A. WwW. Ralston, J. Org. Chem. 9, 329 (1944). 





The Triplet State of Ethylene 


D. P. CRAIG 
University College, London, England 
December 9, 1947 


N dealing with the separation between the ground and 

triplet states of ethylehe, Laidler! has drawn attention 

to the remarkable agreement existing between the experi- 

mental value of 74 kcal. obtained by Lewis and Kasha,? and 

the calculated value of 72 kcal. found by Hartmann*® by 
the method of antisymmetrical molecular orbitals. 

It is not, however, indicated in Hartmann’s paper what 


values were used for the various integrals that occur, and 
in treating the same problem by the same approximation 
during the course of other work, the writer has obtained 
a considerably greater value for the singlet-triplet separa- 
tion. Using values for the integrals obtained by interpola- 
tion between the values given by Sklar and Lyddane,‘ the 
singlet-triplet separation comes out to be 95 kcal. or 
greater, the exact value depending on the assumptions 
made about the configuration of the triplet state. There- 
fore, there is probably not the good agreement between 
calculated and experimental values that Hartmann’s value 
suggests. Hartmann’s value for the excited singlet state, 
given in the same paper, appears to be in error due to the 
use of a wrong wave function, and the corrected value again 
proves to be considerably greater than the experimental. 

The values for the integrals used in these calculations are 
so closely in agreement with those given recently by 
Griffing’ that there is no need to quote them here. 

K. J. Laidler, J. Chem. Phys. 15, 712 (1947). 

2G. N. Lewis and M. Kasha, J. Am. Chem. Soc. 66, 2100 (1944), 
These authors commented that Hartmann’s value — improbably low. 

3H. Hartmann, Zeits. physik. Chemie B53, 96 (19 


43). 
4A. L. Sklar and R. H. Lyddane, J. Chem. Phys. 7, 374 (1939), 
5 V. Griffing, J. Chem. Phys. 15, 421 (1947). 





Erratum: Infra-Red Spectrum and Depolariza- 
tion Factors of the Raman Lines of Spiro- 
pentane and the Raman and Infra-Red Spectra 
of 1,1-Dimethylcyclopropane, 1-Methylcyclo- 
butene, and Methylenecyclobutane 


Forrest F. CLEVELAND 


Department of Physics, Illinois Institute of Technology, 
Chicago 16, Illinois 


AND 
M. J. Murray AND W. S. GALLAWAY 
Universal Oil Products Company, Chicago, Illinois 


[J. Chem. Phys. 15, 742 (1947)] 


UE to a faulty distillation and an oversight, the infra- 

red spectrum of methylenecyclobutane reported in 
Table IV of the above paper contains spurious bands at 
782, 990, 1051, and 1155 cm™!, which correspond to the 
prominent bands of spiropentane. If these are deleted, the 
remaining bands are those of methylenecyclobutane. (It 
‘should also be noted that the strongest band of spiro- 
pentane is nearly coincident with the 877-cm~ band of 
methylenecyclobutane.) The Raman data, however, were 
obtained with another sample which was pure enough not 
to show the strongest line of the spiropentane. 

This correction of the experimental data necessitates the 
following changes in the assignments for eight of the weak 
bands: the combination or overtone frequencies whose 
values are given in the last column of Table IV as 801, 
1509 or 1510, 1883, 2013, 2061 or 2063, 2118 or 2109, 2177 
or 2188, 2310, and 2416 should be deleted; they may be 
replaced by the combination frequencies 1167 —373 =794, 
354+1167 =1521, 728+1167 =1895, 954+-1070 = 2024, 895 
+1167 =2062, 954+1167=2121, 1029+1167 =2196, and 
1029+ 1395 =2424. Thus eight of the nine carbon skeleton 
frequencies and eighteen of the twenty-four C—H fre- 
quencies are now tentatively assigned. 
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An Equation of State for Gases at Extremely 
High Pressures and Temperatures from the 
Hydrodynamic Theory of Detonation 


STEWART PATERSON 


Imperial Chemical Industries Limited, Explosives Division, 
Stevenston, Scotland 


December 1, 1947 


HE thermal equation of state for any two-parameter 
substance is included in the general form: 


pvu=nRT+ pa(v, T), 


where p, v, JT, n, R represent pressure, specific volume, 
temperature, number of moles per g, and gas constant per 
mole, respectively. In attempting to determine the equa- 
tion of state for gases at high pressures and temperatures 
from the observed velocities of detonation (D) of explosive 
compounds, at various densities (A), P. Caldirola! has 
assumed that the “covolume” a(v, 7) can be regarded 
under such conditions as a function of v alone. This as- 
sumption, criticized by S. R. Brinkley,? finds support in a 
recent valuable paper by M. A. Cook* on the same topic. 
Cook’s arguments in favor of the temperature independence 
of a are partly of an a posteriori nature, deduced from the 
coincidence of the a(v) curves for explosives yielding widely 
different calculated detonation temperatures, and as such 
appear to carry considerable weight. The relative con- 
stancy of the theoretical second virial coefficient between, 
say, 2000 and 6000°K provides independent support for the 
same view, which is largely shared by the writer. 

Cook, however, also claims to have evaluated a directly 
from the experimental detonation velocities, not only on 
the assumption that a=a(v), but also by an independent 
and exact method which allows for possible variation 
with 7. If this were so, it would of course enable one to 
dispense with extraneous and a posteriori arguments re- 
garding the temperature dependence of a. However, it does 
not seem that Cook has in fact carried out such an “‘exact” 
solution, nor indeed that it is possible to do so on the basis 
of experimental velocities for one explosive alone. For since 
a(v) can evidently be determined to fit the observed veloci- 
ties for one explosive, these velocities alone do not allow us 
to discriminate between such a solution and any alternative 
solution, e.g., Kistiakowsky and Wilson’s‘ in which a 
depends also on T. 

Cook’s position is somewhat obscured by the large num- 
ber of interdependent equations which he writes down. The 


solution depends, in fact, upon Cook's Eqs. (2), (6), (9), 
and (16), together with any two of (10), (11), (12), and (13); 
for these last four are not independent, any two being 
immediately deducible from the remaining two. Cook ig- 
nores Eq. (9), which amounts to discarding the Chapman- 
Jougnet condition. Then, since W appears only in (2) and 
p: only in (6), we can get these equations also aside, leaving 
(16) and any two of (10), (11), (12), and (13) to be solved 
for T2, ve, a, B. 

These are apparently the equations from which it is 
proposed to deduce 72, v2, a, 8 by successive approximation. 
Cook evaluates a from (13) and then vz from (12), there- 
after calculating 6 from (11). Since (11) follows at once 
from (12) and (13), this process would clearly not be 


expected to define a unique solution. On the other hand, 


if Cook commenced these ‘‘successive approximations” 
with the values of 8 given by “‘method b”’ which as he says 
preceded ‘‘method a’’ in the development of the theory, it 
seems natural that method b should be confirmed. 

As indicated above the clearest proof that Cook has not 
accomplished an exact solution, and that such is not pos- 
sible on the basis of the experimental velocities for one 
explosive alone, is given by the fact that equations of state 
markedly different from his own are capable of*giving good 
agreement with observed velocities. Since Cook’s analysis 
is based solely on the D(A) relation, it is clear, therefore, 
that this relation alone cannot lead to a unique equation 
of state. If it does not provide proof that Kistiakowsky and 
Wilson’s equation, for example, is not correct, it evidently 
cannot prove that an alternative equation is correct. 

Regarding ‘method 6” itself, a further point may be 
made. It is stated by Caldirola, and seems implied by Cook, 
that one can “unequivocally”’ determine a(v) from the 
known D(A) relation. One of the five equations available 
to define v2, p2, T2, W, and a@ is however a differential 
equation in a, and the solution therefore would be expected, 
in general, to embrace not one only but an entire family of 
a(v) curves. That this is in fact the case may be confirmed 
by a step-by-step integration, as distinct from the method 
of successive approximation used by Cook, which should 
and does lead to the envelope of the family. Physical con- 
siderations certainly appear to favor this singular solution, 
but the point seems to deserve attention. 

1 P, Caldirola, J. Chem. Phys. 14, 738 (1946). 

2S. R. Brinkley, J. Chem. Phys. 15, 113 (1947). 


3M. A. Cook, J. Chem. Phys. 15, 518 (1947). 
4G. B. Kistiakowsky and E. B. Wilson, OSRD Report 69. 





